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PREFACE TO THE ENGLISH TRANSLATION 


Å SERIES of admirably clear and instructive courses of lectures, 
covering essentially almost the whole field of Theoretical Physics, 
were delivered by Professor Lorentz at the University of Leiden. 
Seven of these courses, bearing on subjects of fundamental 
importance, were, up to 1922, edited by Lorentz’s pupils and 
published in the Dutch Janguage. During Professor Lorentz’s 
visit in the United States in 1922 the plan was conceived to make 
these courses of lectures accessible to. English readers. 

After some delay necessitated by circumstances the present 
Volume I. is being issued. It contains the English version of 
two of these courses of lectures, viz. Aether Theories and Aether 
Models, edited in the Dutch language by Prof. H. Bremekamp, 
and Kinetical Problems, edited by Dr. E. D. Bruins and Dr. 
J. Reudler. To preserve the peculiar charm of Lorentz’s own 
style and exposition the translation of these lectures has been 
made as literal as was compatible with the nature of the English 
language. Only a very few changes or explanatory additions 
were made in the text. The latter are placed in square brackets. 

Volumes II. and IIJ., now in preparation, will contain the 
remaining five lecture courses, namely, Thermodynamics, Entropy 
and Probability, Theory of Radiation, Theory of Quanta, and The 
Principle of Relativity for Uniform Translations (Special Theory 
of Relativity). 

Recently, 1925, one more course of lectures, on Maxwells 
Theory, was brought out at Leiden, and arrangements will be 
made to include its English version into this publication. 


L. S. 
A. P. H. T. 
ROCHESTER, N.Y., 
July 1926. 
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AETHER THEORIES AND AETHER MODELS 
(1901-1902) 


I 
ABERRATION OF LIGHT 


THE first question which presents itself when we try to form an 
idea of the nature of the aether is that concerning its relations 
to ponderable matter. More especially we will consider the 
question whether a moving body, such as a planet, does drag 
the surrounding aether. The theory of the aberration of light 
may give us in this respect some information. 


1. Stoxrs’ Tneory: THE EARTH DRAGS THE SURROUNDING 
AETHER 


Stokes imagines the acther streaming in the neighbourhood 
of a planet and determines the direction in which according 
to this hypothesis a star should be seen in the following 
way. If we consider a wave-front near the planet, we can 
determine any of its successive positions by the Huygens con- 
struction, provided we take account of the velocity of the 
acther at each point of the original wave-front. Since the 
velocities at different points are different, the wave-front will 
be slightly tilted, and since the direction in which we see the 
light source is determined by the normal of the wave when it 
reaches our eye, this rotation will give us an explanation of the 
aberration. In order to find the amount of the rotation thus 
produced, we note first of all that the translational [orbital] 
velocity p of the earth is 10 times the light velocity V. Terms 
of the order (p/V)* will be neglected. 

Let the plane AB (Fig. 1) represent the wave-front at the 
instant ¢, and let our z-axis OZ be perpendicular to this plane at 
a point M. Then, if wp, vo Wo be the velocity components of 
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the aether at M, we can write for the velocity components at a 
point x, Y, Zə of AB chosen in the neighbourhood of M, 


Ou 
+-9, 


u =o +5, T dy 


Ov v 
V = Vo Tae Pay 
Ow 


Ow 
w=w ta; TE yY 


O 
M 





Z 
Fıc. 1. 


During the time di the point will be displaced to the new 
position 


+(u Hero y)de +( Ov )de 
T+ | Uo taz ay! sY en i 


(w+ ow | Ow ya 


The plane drawn through all the ‘aie thus obtained is parallel 
to the new wave-front. Its equation will be 


Z = Zg + (uy +5 w, Ej wy it 


This relation, in whose second term x and y stand for the old 
co-ordinates, will hold, provided the differences between these and 
the new co-ordinates are very small. The direction-cosines of 


this plane are 


Ow Ow 
nd, — Dy dt, 1, 
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so that the wave normal makes m the X-axis the angle ” 9 + di, 
and with the Y-axis the angle 5 gä The changes of de a 


first direction-cosines, which ns were zero, that is to say, 
the magnitudes 


Ow Ow 
= aut and — dy 
can be taken as a measure of the rotation of the wave-front 
about the y- and the z-axis. 
The wave-front propagates itself in the time dt over a distance 
Vdt, so that the rotation, per unit length, has the components 


dt, 


1 ðw 10w 
Væ Vy 
The total rotation undergone by the wave-front until it will have 
reached our eye will, of course, be expressed by an integral. 
Phe preceding reasoning gives us the contribution to that integral 
for a wave-front which is normal to the adopted Z-axis. The 
deviations from that orientation remain, however, so small that 
these results can be used for all the elements of the integral. 
Thus the components of the total rotation will be 


OU 7, a 
-y A 


Ox 


where the integrations are to be extended from a point at which. 
the earth’s motion is still imperceptible up to the observer's 
station P. 

Now, it is actually possible to account for aberration by 
means of this rotation of the wave-front, provided that the 
motion of the aether is assumed to be irrotational. There exists 
then a velocity potential ¢, such that 


Thus also Cw/dx=Cu/dz, Ow/ey=Ov/dz, and our expressions for 
the rotation components become 


Ou 
- pfd- a nd - - fide- - y 
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It will be readily seen that this result agrees with observation, 
and with the elementary theory as well, if it be assumed that 
the velocity of the aether at the observing station P is equal to 
the translational velocity of the earth. 


2. VELOCITY POTENTIAL IN AN INCOMPRESSIBLE AETHER 


Here, however, a serious doubt arises, whether the existence 
of a velocity potential is compatible with the requirement 
that the velocity should be equal all over the surface of the 
earth (viz. equal to the carth’s velocity in its annual motion), 
whereas the assumption of such a potential is indispensable for 
obtaining the correct value of the aberration. As a matter of 
fact, there is nothing uncommon about a motion for which a 
velocity potential exists. In a frictionless liquid in which there 
are, at a given instant, and therefore also at any later time, no 
vortices, no other motion, in fact, is possible. But, as was just 
mentioned, the assumption of a velocity potential for an incom- 
pressible aether cannot be reconciled with the other requirement 
that the velocity at the earth’s surface should be everywhere the 
same, in size as well as in direction. In fact, in an incompressible 
aether the motion is completely determined, if it be assumed 
that a velocity potential exists and that the normal velocity 
component of the aether at the surface of the earth is every- 
where equal to that of the translation velocity of the earth ; but 
the tangential components are then found to be different. If 
the origin is at the centre of the earth, the X-axis in the direction 
of the translational motion of the earth, and if the velocity of 
this motion is p, then all requirements are satisfied (and according 
to the theory of the Laplace equation this is the only possible 
solution) by the potential 


Cx 
$= pa’ 
so that 
ni IE ,_ ry Z 
u=0(2- 5 h = -30% w= -305 


where C is determined by the condition that the normal com-. 
ponent of the velocity at the surface of the earth should be 
equal to the component of p in the same direction. If R be the 


I ABERRATION OF LIGHT 7 


radius of the earth, and @ the angle between the direction r and 
the X-axis, then, since 0¢/or = — 2Cz/r’, the last condition gives 


-20 cos 6/R? = p cos 9, 
and therefore C= —}pR’. 


Now, in order to see that the velocity of the aether cannot be 
the same all over the surface, viz. equal to that of the earth, it 
is enough to notice that for points of the surface in a plane 
through the centre and perpendicular to OX the velocity along 
the X-axis becomes —3p, whereas at the intersection point of 
the surface with the X-axis we find the velocity p. In order to 
avoid this difficulty, one might perhaps take advantage of the 
circumstance that the velocity potential need not exist in the 
whole space around the earth, since we are concerned only with 
a limited region. This, however, would lead to very artificial 
and improbable concepts. 


3. PLANCK’S THEORY. (COMPRESSIBLE AETHER 


It was shown by Planck how these difficulties can be met by 
giving up the incompressibility of the aether and by assuming 
that the acther is subject to the earth’s attractive force. To 
investigate the details of the motion, it will be convenient to 
attribute to the whole universe a velocity —p, while consider- 
ing the earth to be at rest. We have then to investigate the 
disturbance produced by the fixed sphere in the otherwise 
uniform aether stream. Planck adjusts his hypotheses so that 
the velocity at the surface of the sphere should be small. In 
order to see how this requirement can be satisfied, it is enough 
to realise that a plane through the centre of the sphere perpen- 
dicular to the direction of the stream must be traversed in a 
given time by the same amount of aether as in the absence of 
the sphere. The quantity of aether which in the latter case has 
to stream through the circle cut out by the sphere from that 
plane, must find its way through the remaining part of the plane 
surrounding the circle, and especially through that nearest to 
the circle. This is, with a small velocity, only possible if the 
density in the neighbourhood of the sphere is large. It is 
necessary, therefore, to assume that the aether [to be thus 
condensed] is attracted by the earth. It must further be assumed 
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that the velocity of light in this very condensed aether is the 
same as in the aether of normal density. Planck assumes also 
that the aether while being cundensed behaves like a gas. 

We will now work out in detail these ideas. Let the Z-axis 
be chosen in the direction of motion of the aether at infinity. 
Let, as before, & be the velocity potential and, therefore, the 
velocity components u=0¢/Cx, v=0¢/¢y, w=0g/0z. The dis- 
tance of a point (x, y, x) from the centre of the earth will be 
denoted by r, and the radius of the earth by rọ Let p be the 
pressure, and %k the density of the aether, and let k/p= 
be assumed to be constant. Further, let V be the potential of 
the attractive force, per unit mass, so that the components of 
this force, again per unit mass, are -OV/ex, -OV/oy, -OV |. 
Then V=g/r, and the force, which is radial, is, per unit mass, 
-0V /0r =q]r?, and therefore at the surface, g/752?= -g. Thus q 
is determined, and we have 

2 
V=- 
r 

Let us now write down the equations of motion for a stationary 
state, in which the velocity and the density at every point have 
always the same values. First of all, the equation of continuity, 


le) Hall) + talk oh) =O. pe 2) 


Next, the equation of motion, 


uot toes wou = OV lop 

Ox ay” dz De kon’ 
where the left-hand member represents the acceleration along the 
X-axis, since, the state being stationary, Cu/%t=0. The right- 
hand member is the force per unit of mass in the same direction. 
Two similar equations will hold for the Y- and the Z-directions. 
Owing to the existence of a velocity potential an integral can be 
found at once. In fact, the equations of motion can be written 


bai ae) +58) +E) ea, 


cx | \Ox 02 dx kõr 
etc., whence 


[Pov +y HGE) + (+ =constant. . (2) 
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Let us now assume that the term containing the square of the 
velocity can be neglected. Then, 


J 2 + V =constant. 


This equation shows that, to the assumed degree of approxima- 
tion, the density distribution in the aether can be taken to be 
the same as in the state of rest. Thus we can write 


1 To 
-log p-g 5 =const., 
or also 


1 ro 
. log k -g „ sonst. 


The integration constant can be expressed in terms of the aether 
density k, at the surface of the earth, 


] ry 1 
> log k -g - E log ky - gro» 


instead of which we can write 


log nal, _ Sy 


where a = ugr. 

For r=» we find the limit value of k, 

ka = kge te, 

In the next place, Planck determines the velocity potential 
from the equation (1), where now k is a known function of the 
co-ordinates. What is required is a solution representing a 
motion which at large distances reduces to a stream along the 
Z-axis. Let vy be the velocity of this stream.* Then, at large 
distances, =z. This suggests the form ġ =zf(r). Then, 

Cp _ 22 df Op _y20f op _ 2 df -f, 


— So — 


dx rdr dy rros rdr 


h 2 2 Of | x2 af vz af 
(ee rdr r? dee y2 dr? 











* This was previously denoted by —p. In accordance with Planck’s 
notation p will now be used for the pressure. 
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Op _zdf yzdf ye df 


— —  —— 


dy? rd P dr r? dr? 
h _3z df _2 df 2 df 


02 rdr Pdr dr? 
and the equation (1) becomes 


df ptf) ak ‚az df ok yadf kf #4) =0. 


+2 2 dr 


(g ta ro y ru oN ae 


Substituting here the value of log k just found, we have 


ees ot „Iwazdf, iy yz af 2 (2 df 
Spar ae “Vr dr’ rr dr? (eg). = 


or, after some reductions, 
df = A a 
dr *\r "r2/dr pt = u 


The solution of this equation is 


fag- ag 


The constants a and b will be determined by considering the 
state of things at infinity and at the earth’s surface. For r= 
we must have ¢ =2y, and since ¢ =</(r), 


f(2)=y, 


and b-a=y. At the surface of the earth the aether cannot 
have any velocity perpendicular to that surface, 7.e. 0¢/0r=0, 
and since 


Fe taf’), 


or 


we have, at the earth’s surface, f(r) +7f'(r)=0. Substituting the 
value of f(r) and putting r=7,, we find 


-a +be~ ire +02 (5 + 1)e-em =Q. 
To `2To 


ae 2 
Thus a =b á + + 1 Je=ei so æ. g a A) 
o To 
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Now for the sliding of the aether along the earth’s surface. 
In order to find this, we have to calculate the Z-component of 
the velocity. From ¢=2f(r) follows 


nee), 


02 


and, at the surface of the earth, where f(r) +rf’(r)=0, 


Ob sas 
Ay =f(r) (1 - 3 = sin? Of (r) 
where @ is again the angle between r and the Z-axis. If w is 
the velocity in the tangential direction, then, since 0¢/0z = sin 6, 
we have 


=sin @. f(r), 
whence, by (3) and (4), 


2 
wo =sin 0 (eet a + Ue = -1) ton + 1 been = p u sin 0, 
Our aim now is to make the coefficient of sin @ small in com- 
parison with y. This will be the case if a/r, is large. For then, 
in virtue of (4), the constant a will be small as compared with b, 
and therefore, since b -a = y, b will be approximately equal to y. 

Planck has made some numerical estimates under different 
assumptions. Let us first take u=k/p and «=ugry as great 
as for air at 0° C. and for g the value which holds for 
ordinary ponderable matter. Then the consequences hold for 
the air actually surrounding the earth. Thus, a/r,=800, and 
the ratio of the density at the earth’s surface to that prevailing 


at a large distance is 
etlo = e800. 


The result of the computation is therefore, practically, that our 
atmosphere can well be kept by the earth, though at the same 
time the theory shows that some residual sliding is unavoidable, 
provided that friction be disregarded. Now, with reference to 
the aether we need not choose our assumptions so as to be 
driven to such an extravagant condensation. It can easily 
be calculated how far we must go with this.’ „The aberra- 
tion constant is known only up to 4 per “sent. In other 
words, the effect of a [sliding] velocity smaller than xov 
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of the velocity of the earth cannot be detected. It is enough, 
therefore, that the velocity of the aether at the earth’s surface 
should remain below this lmit. For this purpose it is only 
necessary to have a/r,>11. Even so we are still left with a 
condensation e™! [or about 60,000]. The natural question arises 
whether such a strong condensation could not be tested in an 
independent way. Notice that the condensation around the sun 
and similar bodies will be considerably greater, since the attrac- 
tion is proportional to the mass. The aether condensation due 
to the presence of the sun will also increase somewhat the con- 
densation e””* at the earth’s surface. A number of questions 
suggest themselves here, but must be left unanswered. 


4. FRESNEL’S TuEory. Fixep AETHER 


We now turn to the alternative hypothesis, that of a fixed 
aether. This was already assumed by Fresnel. The aether 
must then be able to pass freely through the earth. The atoms 
themselves may well be impenetrable, provided, however, that 
they are assumed to occupy but a small fraction of the total 
volume. We may also assume that the atoms are not impene- 
trable, following a line of thought according to which the atoms 
or their constituents are but special modifications [singularities] 
of the aether. 

In order to explain aberration we again apply Huygens’ 
principle. Let us take the simple case in which the position of a 
star is being determined by means of a primitive sighting appar- 
atus (without lenses), and let the presence of the atmosphere be 

disregarded. The explanation is 

2 # Ss then the same as on the emission 
theory. The propagation of the 

2 wave-front which at a certain 
e Fic. 2, D instant reaches the aperture AB 
(Fig. 2) is unaffected by the 

motion of AB. A simple construction gives us that portion 
of the screen CD which receives the aether disturbance, 
and it is manifest that in this way one falls back to the 
older explanation. But this reasoning cannot be applied when 
we are concerned with refracting surfaces. In this direction 
many experiments on aberration were made. As a typical case 
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we can consider the famous experiment of Arago. If ab (Fig. 3) 
is the direction in which, with the aberration, a star would be 
seen directly, then cd, the direction 

in which the star is seen through a a 6 
prism, will coincide with that of the 5 

refracted ray belonging to ba as 
incident ray. Another experiment Fic. 3. 
[Boscovich - Airy] proves that, in 

observing a star, a telescope filled with water has to be set in 
the same direction as an ordinary telescope. To sum up, these 
experiments show that all refraction phenomena are the same as 
if there were no aberration. 


5. DRAGGING COEFFICIENT 


In order to explain Arago’s experiment Fresnel introduced 
the dragging-hypothesis, which amounts to this : 

Let W be the propagation velocity of light in a given medium, 
when this is at rest. Then, if the medium is moving with a 
velocity p in the direction of light propagation, the velocity of 
light relatively to the aether is, according to Fresnel, not W +p 
but W +kp, where k is a fraction. In other words, the ponder- 

able medium behaves as if it 
| dragged the light with a velocity 
which is only a certain fraction 


H ARB í 
of its own velocity. The co- 
efficient & must then be 
cy /p—_—— kele 
n 
EGF where n is the refractive index 
id of the medium when at rest. 
IG. 4. 


For n=1 we have k=0, as it 
should be; for when n tends to unity, the medium becomes 
indiscernible from the aether. 

This expression for the coefficient k can be easily found by 
the following reasoning. Let us suppose that the position of a 
celestial object is being determined by means of a sighting 
apparatus consisting of a screen with an aperture AB (Fig. 4), 
which is fixed in the aether, and of a second screen with an aper- 
ture CD, behind which is placed some ponderable medium, as 
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e.g. glass, and finally of a third screen which receives the light. 
Since the refraction is such as if the apparent direction of the 
ray were the true one, the segment EF of the third screen which 
receives the light will be found by putting 


sin (ECG) =" sin (ACH). 


This is the result of the experiment. Now, let v, be the velocity 
of light in the aether, v, the light velocity in glass (at rest), and 
p the velocity of the earth, and let us consider the case in which 
the incident light is perpendicular to the first screen. Then, 
since the deviation from the normal remains small, 


sin (ACH) =: tan (ACH) = p/p, 


and therefore, LECG=p/nv, and EG =lp/nv,, if l=CG. The 
segment thus determined differs from that which we would find 
if, starting with the wave-front CD, we applied Huygens’ prin- 
ciple and used in it v, as light velocity ; for then we should find 
a segment E’F’ such that 


E’G:C@G=p:v=Pp:Vln, 


and therefore, E’G=Inp/v,. To reconcile this with the experi- 
mental result we have to assume that within the time //v, the 
glass has dragged the light over a distance 


FE'E = E'G — EG =Inp/vy — lp/nrp, 
that is to say, with the velocity 


Pa): l aai -3) 


6. TnEoRY or ABERRATION 


lf this dragging coefficient is assumed, it can be proved that 
all the phenomena of refraction, etc., are such as if there were no 
aberration. For this purpose we use the artifice of imparting to 
the whole system a velocity equal and opposite to that of the 
earth. The earth is then at rest, while the aether has everywhere 
the same velocity. This proof can be given for the more general 
case of any motion of the aether, provided it has a velocity 
potential, Our first business is to find out how the wave-fronts 
and the light-rays are to be determined. The elementary wave- 
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front spreading out from a point O would, after a time ż, in a 
stationary aether, be a sphere of radius v. Now, if p be the 
velocity of the aether at O, and if it varies continuously from 
point to point, then, neglecting infinitesimal terms of the second 
order, it can be shown that this sphere is simply displaced as a 
whole over the distance pt. The successive wave-fronts are the 
envelopes of the spheres thus determined, every time to be con- 
structed around the points of the preceding wave-front as centres, 
and taking account of the aether velocity at those points. 

To find the wave-fronts in a ponderable medium, let us 
first consider the case in which the aether is at rest, while the 
ponderable medium has the velocity -p. The elementary wave- 
front is then, after a time t, a sphere of radius v,¢ which is dis- 
placed as a whole over a distance — pi(1 —1/n”). If we now give 
to the whole system the velocity p, the elementary wave-front 
will be displaced relatively to the 
ponderable medium, thus brought to 
rest, over the distance 


1 

p = kpt, 
where «=1/n?. This coefficient can 
be said to determine to what extent 
the light has been dragged by the acther. Let now O (Fig. 5) 
be a point of a wave-front. The elementary wave - front 
emanating from O is a sphere with O’ as centre, where OO" = x pt, 
and of radius v,t. Let A be the point of contact of this sphere 
with the envelope which is the new wave-front. Then OA will 
be an element of the light-ray through O, and we find for the 
propagation velocity along this ray, relatively to the earth, 
w=OA/t. Thus, if @ be the angle O’OA, we have 


v? =w? — 2rwp cos d + Kp, 





A 


Fia. 5. 


whence, neglecting higher powers of p, 
k?p? 
w =V; + Kp cos 6+ 95 (co®9-1). . . (5) 
1 


_1),_*p Kp 
= 1 a cos O+ oya (cos? s+). . « (6) 


These formulae hold for any continuous space-distribution of the 
velocity p. 
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Having thus found the velocity along the ray, we can readily 
determine the light path relatively to the earth. For this pur- 
pose we make use of a theorem according to which light follows 
that path to which corresponds the shortest time, a theorem 
which in the present case can also be easily deduced from 
Huygens’ principle. In fact, let AB be the light-ray determined 
by a Huygens construction, and let us consider any other path 
between A and B. Both paths cut the successive wave-fronts. 
The intersection points of the latter with the actual light-ray 
may be called corresponding points. Let S and S’ be two 
successive positions of a wave-front separated by a very short 
time-interval, and let P’, Q’, etc., be points of S’ corresponding 
to P, Q, etc., chosen arbitrarily on S. Then the time of passage 
along PP’, QQ’, etc., always with the ray velocity belonging to 
these lines, will be the same. The time required for covering 
any other path drawn from S to S’, again with the ray velocity 
belonging to it, that is to say, a path laid through non-corre- 
sponding points, will be longer. This appears from Fig. 6, 
in which A is a point of the first wave-front and A’D lies in 

the second wave-front, A’ being the point 

A Cj, in which the latter wave-front touches the 

elementary wave A’C circumscribed around 

a’ A, so that A’ and A are corresponding 

Fic. 6. points. Since the light time for AA’ is as 

long as that for AC, that for AD is evi- 

dently longer. Let us now imagine between the points A and 

B the whole series of successive positions of the wave-front. 

Then, remembering that the light-ray is the locus of corre- 

sponding points, while this cannot be said of any other path 

from Á to B, it will: become manifest that the time along 

the light-ray is indeed the shortest among all paths leading 
from Á to B. 

This being proved, the light path can be determined in the 
following way. The time required for covering an arbitrary 


B 
path between A and B is | ds/w, the integral to be extended 
A 


over the path. Now, this expression has to be a minimum for 
the actual light path. 

Let us first take the case of a non-homogeneous (but isotropic) 
ponderable medium, without surfaces of discontinuity, however. 
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Then «, upon which w depends, is variable, since n as well as v, 
are variable. Further, w depends also on p and 6, which in their 
turn depend on the state of motion of the aether, of which we 
have only assumed that it has a velocity-potential ¢. 

Thus, taking account of (5), and neglecting terms of the order 
»2|v?, we find for our integral 


Bds Kp Bis 8] Pag 

JaC- cos )-[ - f za afs 
Rds -fi 1 06 as. 
“Jav, Janw Os 


Now, nv, =% being constant throughout the space, the last 


integral is Z$ B— $a), and therefore independent of the path 


between A and B. The minimum property is thus influenced 
only by the first integral on the right hand, and since this con- 
tains no trace of motion, the light path will be the same as if 
there were no motion of the aether. Notice that this proof is 
based, first, on the assumption of Fresnel’s dragging cocfficient 
and, second, upon the existence of a velocity-potential. 

It follows from the Huygens construction that the minimum 
property holds also for the case in which the light passes from 
one to another medium. Consequently, the light path can in 
this case be found by the same reasoning as before. In fact, 
considering any path from A to B which cuts the refracting 
surface at C, and putting 1/n?v,? =, we have 


[o-fe* +[% 


Cds 
-J — Bde — ba)+] i m(n- $e) 


Vy 
-[ © -Mbn-64). 
A “Ll. 


In order that this should be a minimum, it is enough to make 
the first term a minimum. Consequently, the light-rays in 
reference to the earth obey the ordinary laws of refraction. 
We have sti!l to consider the relation between the ad 
VOL. I 
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and the wave-front. For these are no longer perpendicular to each 
other. Let O (Fig. 7) be a point of a wave-front, let OO’ = pxt, 
and let P be the corresponding point of a successive wave-front, 
so that OP=wt will be an 
Q element of the light-ray, and 
O'P =v,t. If now the wave- 
front, and therefore the 
wave-normal (O’P), at every 
point be given, then the 
direction of the light-ray 
can be found by combining the velocity v, along the normal 
with the velocity «p. Vice versa, if the ray be given, the direction 
of the wave-normal will be found by combining the velocity w 
along the ray with the velocity —«p. In the latter construction, 
and up to terms of the order p?/v?, the velocity w can be replaced 
by v, taken in the same direction. 

Let us still consider a luminous point L in a homo- 
geneous medium, and let us determine the wave-fronts by the 
indicated construction. The light-rays being straight lines 
diverging from J, it is required to find the direction of the 
wave-normal at any point A. Tf x, y, z be the co-ordinates of 
A, with Z as origin, the components of the ray velocity v, at 
the point A are 





x 2 


"Oy, Ton v 

pr Dpr 

and these must be compounded with -«p, the velocity of 
dragging. The components of the latter are —«Ad/ox, — rdpldy, 
—xop/cz. Thus the components of the resultant will be 


0 0 0 
E-r), gur -Ah p Mr -xe 


and since this resultant has the direction of the wave-normal, 
the equation of the wave-front becomes 


vr — Kp =const. 


In the neighbourhood of L the potential & can be considered as 
a linear function of the co-ordinates, so that the last equation 
assumes the form 

vır - «(ax +by +cz) =const. 
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As we already know, the wave-fronts in the immediate 
neighbourhood of J, are spheres whose centres lie at a certain 
distance from L. It is manifest that the equation agrees with 
this, provided terms with p° are neglected, as in fact they were 
in deducing it. Again, waves converging towards a point have 
the same form as would have expanding ones, if the velocities 
in the aether were everywhere reversed. 

In order to apply this reasoning to the determination of the 
direction in which we see a star through a telescope, let us 
consider a wave-front arriving from a star, and let us derive 
from the given wave-normal the direction of the light-ray. It 
is these relative light-rays which we observe in our experiments. 
In the free aether, far away from the carth, the direction of the 
relutive ray is found by combining the velocity v, along the 
wave-normal with «p, where «=1. In this manner we find, 
obviously, the same direction as according to the elementary 
theory. The further progress of these relative rays can now be 
followed up by the ordinary laws of refraction, etc. And we 
have now to orient the telescope so that the rays thus treated 
should converge upon the intersection point of the cross-wires. 
In this way the theory accounts for all the experimental facts. 

Also the diffraction and interference phenomena are such as 
if the earth were at rest and as if we had to do with the velocity 
vı lf we have, for instance, two paths between Æ and B, the 
interference will depend only on the difference of the light times 
along them, and this difference can be determined by using for 
each path element the velocity v,. For this amounts only to 
omitting, for the two paths, the term p(¢ -,), which is the 
same for both. Here again it is only necessary to assume that 
there exists a velocity-potential, but not that the earth is pene- 
trable for the aether. The latter, however, must be assumed for 
all transparent media in order to have a reasonable explanation 
of the dragging coefficient. The simplest, after all, is Fresnel’s 
theory. 


7. MICHELSON’S EXPERIMENT 


We have thus far neglected all terms having p?/v? as factor. 
See formula (5). 

In certain interference experiments, however, the accuracy 
can be pushed so far that these second-order terms have to be 
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taken into account. Such is the famous experiment of Michelson, 
already suggested by Maxwell. Let Æ and B be two points fixed 
on the earth, and let the latter move with the velocity p along 
AB. What is the time taken by light for a complete to-and-fro 
passage ? This will depend on whether the aether shares in the 
earth’s motion or not. In the former case the time in question 
will be independent of the velocity p, but in the latter case the 
velocity of light relative to the earth will be v-p in one, and 
v+p in the opposite direction, and therefore the time required 
for a complete to-and-fro passage [up to fourth-order terms], 

a se + oP 

V-p vtp v `v 

Since Michelson determines this time by means of an inter- 

ference phenomenon, the influence of the last term will still be 
perceptible, though it amounts only to a fraction of the vibration 
period. Let us see how large J must be in order to give an 
observable effect. The shift of the interference fringes due to 
a time lag of ,', of the vibration period T will be just detectible. 
For this purpose we must have 








B! 
oP l 
2 3 = ol: 
hence, if > be the wave-length of the light, 
2 = ae 107A = 600 cm 
B ~ 10 p? ® 


or /=3 metres. 


In Michelson’s experiment a ray travers- 
‘ing to and fro a certain distance in the 
direction of the earth’s motion is made to 
interfere with a ray traversing the same 
distance up and down in a perpendicular direction. Fig. 8 
gives a schematic representation of the apparatus. This consists 
of two fixed mirrors B and B’ and a dividing glass plate A, 
inclined at 45° to the incident light beam. A part of this beam 
passes through the plate, is reflected at B, and after a reflection 
at the lower face of the plate enters the telescope. Another 
part of the incident light is reflected at A towards the mirror B’ 
and thence through the glass plate into the telescope, where it 
interferes with the first partial beam. The experiment consists 


Fia. 8. 
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in turning the whole apparatus, including the light source and 
the telescope, by 90° and comparing the interference fringes in 
the new and the original orientation. The effect of the phase 
difference is thus doubled. Yet no displacement of the inter- 
ference fringes was observed. 

Let us consider the theory of this experiment somewhat more 
thoroughly. We have to compare the times taken by light to 
traverse two different paths be- 1 
tween two points A and B (Fig. 9). “ 

We have already seen that, if only 

first-order terms are retained, the B 
difference of these times is the 

same as if the earth were at rest. : 

Now, to determine the effect of CD 

second-order terms we must take Fia. 9. 

into account that the light path 

between A and B is itself slightly changed by the motion. In 
fact, the light path could be identified with that corresponding 
to a fixed earth only when we confined ourselves to terms of 
the order p/v. Now, however, the slight change of the path 
implying second-order terms need not be negligible. But it 
can be easily shown that even with the inclusion of second- 
order effects the propagation time from A to B can be 
calculated with sufficient accuracy by extending the integral 
/ds/w, with w given by (6), along the light path which would 
correspond to a fixed earth. Let, eg., ACB (Fig. 9) be the 
light path with the earth at rest, and ADB the actual light path 
for a moving earth. Then the integral /ds/w would have to 
be evaluated for the path ADB. But since the value of this 
integral is just a minimum, it will differ from that taken along 
the path ACB only by a quantity of the second order in the 
deviations of the two paths, such as CD. And since, as 
already mentioned, this deviation is itself of the order p?/v?, 
the difference between the integrals extended over ACB and 
ADB will be of the order p?/v?, and can thus be neglected. 

The same holds also for the other light pencil, as indicated 
in the figure by the lines without letters. 

The propagation time from A to B, r=/ds/w, can now be 
calculated as in Art. 6, developing 1/w into a power series of p/v 
and confining ourselves to the first three terms. 
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Thus, r=7,+72.+73, Where, by (6), 


Rn f ds 
1 ano s 
vi 
» cos 8 
nn i; 
a v 
1 


(which has no influence on the phase difference), and 
2,2 
t3=4 | a (1 +cos® O)ds. 
1 


In Michelson’s experiment we can put « =1, and pis, by Fresnel’s 
theory, the velocity of the earth. To the phase change due to 
the turning of the apparatus the integral 7, contributes nothing, 
so that only 7, has to be taken care of, where @=0° or 180° for 
the first, and @=90° or 270° for the second of the interfering 
rays. Ifl be the doubly covered path, then, as we already saw, 


T; = 2pl/v3 


for the first ray, and 7,=p?l/v? for the second ray. Thus also 
the latter is somewhat affected by the earth’s motion. This 
can also be seen in the following 
simple way. A ray impinging 
upon the glass plate at A (Fig. 10) 
returns from the mirror B and 
meets it again in A’, so that 
the distance over which the plate 
was displaced in the meantime 
is AA’=2lp/v. Consequently the 
length of the light path is 
Fia. 10. 2V2 +P p2/o? =21(1 + p2/20). © The 
difference, as compared with a 
fixed earth, is therefore 192/v?, and the difference of -the corre- 
sponding time, 1p?/v?, which is the value just found for 7;. 

As Michelson overlooked the influence of the earth’s 
motion upon the light time along the path ABA’, in his first 
experiment he estimated the theoretical effect twice too high. 
In this case the distance | was 1-2 metre, so that a shift of ‚y of 
a fringe width was expected. The corrected theory gave only 
z% Of a fringe width, and this was below the threshold of reliable 
observability. Michelson, therefore, repeated the experiment 


B 
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with the modification that each of the two beams passed several 
times between the mirrors. In this manner the light path was 
increased to 22 metres and should have given a displacement 
of ‘o of a fringe width (an estimate in which account was also 
taken of the motion of the solar system). But even now the 
result of the experiment was negative in contradiction to Fresnel’s 
theory. It will be kept in mind that the validity of this theory 
is not limited to the case of an aether at rest as a whole, but 
extends also to types of motion of the aether for which there 
exists a velocity-potential. 

This experiment would also clash with any theory which 
attributes to the relative velocity of the aether and the carth a 
value not sensibly smaller than the translation velocity of the 
earth. In what follows, however, we shall have in mind only 
Fresnel’s theory. 


8. CONTRACTION 1N 'THE DIRECTION OF MOTION 


We can explain the negative result of Michelson’s experiment 
by assuming that the length of the arms of the apparatus is 
changed by turning it through a right angle. ‘This change can 
be assumed to be just such as to give to r, the same value for 
both rays. This calls for a contraction of the path in the 
direction of motion, as compared with the perpendicular path, 
such that the corresponding light time should be shortened by 
pj’. The path then has to be shortened by lp?/v?, and there- 


2 
fore | by uP. This dependence of the dimensions upon the 
v 


orientation with respect to the carth’s motion is not as strange 
as it might seem at first. In fact, the dimensions are deter- 
mined by molecular forces, and since these are transmitted 
through the aether, it would rather be surprising if its state of 
motion had no influence upon the dimensions of bodies. The 
nature of the molecular forces is not known to us. Yet, if we 
suppose that they are transmitted through the aether in the 
same way as electric forces, we can develop the theory of this 
contraction, and we then find for its amount just what is required 
for the explanation of the nil-effect of Michelson’s experiment. 
This contraction would amount for the diameter of the earth to 
6-5 cm., and for a metre rod to +} of a micron. 


Il 
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9, MAXWELL S EQUATIONS 


WE will now consider some theories of the nature of the aether. 
Such a theory must, in the first place, explain the electromagnetic 
phenomena. We shall, therefore, begin with an interpretation 
of Maxwell’s equations. Here one can put different require- 
ments. One can content himself with a theory accounting for 
the phenomena in isotropic and homogeneous media, or try to 
include also the anisotropic and heterogeneous media, and so on. 
To begin with, we shall exclude the conductors only, and shall 
thus consider anisotropic as well as non-homogeneous dielectrics. 
The latter will enable us to treat the boundary conditions, and 
therefore also such phenomena as reflection and refraction. 

The magnitudes appearing in Maxwell’s equations are: the 
e'ectric force E, the dielectric displacement D, the magnetic 
force H, and the magnetic induction B. The equations are 


rot H --D, or 
OH, OH, 10D,) 


dy oz c at 
OH, OH, _10Dy 
% s c dat 
aH,_OH,_10D, 
Oz Oy c ot 

Dz , Dy, @De_ 


div D=0, or a t w oz =0,. . . (8) 


and rot E = at 0B or 
c ot 


yo‘ a (D) 
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OH, E, _10B: 


‘Oy Oz c ot 
OH, OH, 10B, 


Do ee we | Og ; . (9 
Oz ör c Ot (9) 
OH, _ OB; _10B, 
a y cù 
OB, OB, OB, | 
div B=0, or - an tO de eo =0. .  . (10) 


Further, we assume linear relations between the components of 


D and E and between those of B and H, 
Dz =€yHe+eyghy +h, Bre=py Hs + pypHy + pish: 
Dy =€2, Hs + €gghiy + €ogh, By = paH s + poH y + Mog, 
D; =egj Er + eggbiy + egal;, Bz = by Hz + Haad y + asl, 


where 


€12 = E21, €23 = E32 €31 "Ei, Mie Mow Hos = H32: H31 = His 
For isotropic bodies these coefficients are all zero, and 


€11 = E22 5E3356, Miu SH SP33 FF. 
We note that the usual kinematics of continuous media gives 
equations of the same form as rot H= t), In fact, if every 
c 


point of a medium is displaced through an infinitesimal distance, 
of which the components £, n, ¢ along three axes can be considered, 
within a small region, as linear functions of the co-ordinates 
7, Y, Zz, then the complete change in the neighbourhood of a 
point can be represented as consisting of a displacement, a 
rotation, and three dilatations or contractions in three mutually 
perpendicular directions. To see this, it 1s enough to remember 
that a sphere is transformed into an ellipsoid, whose conjugated 
diameters correspond to orthogonal pairs of diameters of the 
sphere. This holds also in particular for the principal axes of 
the ellipsoid. We imagine the sphere to be first displaced so as 
to bring its centre into coincidence with that of the ellipsoid, 
then find those diameters of the sphere which coincide with the 
axes of the ellipsoid, bring them by a rotation into their actual 
position, and finally give to these diameters by a dilatation or 
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contraction the actual length of the axes of the ellipsoid. The 
components of the rotation are 


OEM (AE aL) dn PE 
Pl) ar a D 


In a continued displacement we can compare the state at the 
time ¢ with that at £+dt, and the last formulae give us then the 
connection of the components of the angular with those of the 
displacement velocity. 

We may now try to interpret one group of Maxwell’s 
equations by taking the components of the magnetic force to 
be proportional to the displacements in the aether, or assuming 
that the essence of that force consists in such displacements. 
In symbols, let us put 


H,=m£, H,=mn, H,=mt. 
att 0 u OD: 5 
Then ‘ay emp, —," =2cmg, ap zent, 


Thus, if the observed magnetic force is in its essence a displace- 
ment of the aether particles, proportional to it and having the 
same direction, the nature of the dielectric displacement current 
consists in rotations which are the result of the aether displace- 
ment. But we are then checked by the difficulty that there is 
in such a picture no place for a constant dielectric displace- 
ment, such as occurs always in electrostatics. Moreover, in 
a constant electric field without magnetic force the dis- 
placements č, 7, & vanish, so that there is nothing to 
distinguish such a state of the medium from one in which 
an electric field is absent. 


10. Toe Maenetic Force as VELOCITY. NEUMANN’S 
THEORY OF LIGHT 


We can avoid this difficulty by identifying with the aether 
rotations not the displacement current but the dielectric displace- 
ment itself. Then our equations have to be made to agree with 


op (PL a 
E Haga mae) 
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Consequently, we have to put 


H, =mi, H, =m", Hy=ma, >. ay EE) 


D, =2cmp, Dy =2cmq, D,=2cm. . . (13) 


Thus, wherever there is a magnetic force, we must imagine 
an aether velocity in the direction of this force and proportional 
to it, and we have to look for the dielectric displacement 
in the rotation due to or associated with that velocity. In a 
permanent magnetic field, as e.g. around a steel magnet, we have 
thus to imagine the aether streaming along the lines of force. 

We may also notice that in such a case of continual motion 
the displacements E, 7, £ would not remain infinitesimal. 

It is important to keep in mind that the interpretation here 
given implies necessarily that the coefficient m has also in a 
non-homogeneous medium throughout the same value. Other- 
wise the values (12) and (13) would not satisfy the equations (7). 
The expression 

0H, 0H v, 
Ay az’ 


for instance, would no longer have the value 2m?p/dt, since to 
this the terms 

om OG om on 

dy OL Oz Ot 
would have to be added. This difficulty would assert itself 
especially at such places where the properties of the medium, 
and therefore also the value of m, vary rapidly from point to 
point, as at the boundary of two media. This theory, therefore, 
implies that the coefficient of proportionality between the velocity 
of the acther particles and the magnetic force 1s always and 
everywhere the same. 

We have still to explain the second group, (9), of equations. 
Before doing so we have to consider the energy relations. The 
magnetic energy is a quadratic function of the magnetic force, 
and has thus to be interpreted as kinetic energy, while the 
electric energy will become potential energy. The kinetic energy 
per unit volume, with p as density, is 


ee, 
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and the magnetic energy per unit volume of an isotropic medium, 


with y=], 
-mO 


Thus we have to assume p=, so that also the density will be 
everywhere the same. But the case is different if ~ is not equal 
to 1. In such media the magnetic energy per unit volume is 


o 0 NAN 

3BH = Sum? (ey (a) (a) } 
and therefore p=m?, t.e. the density must be proportional to 
the permeability. According to this theory the characteristic 
feature of iron, for instance, is a large aether density within it. 
For anisotropic media the relat ons are somewhat more compli- 
cated, inasmuch as for these the magnetic energy per unit 
volume is 


Mund + bogll,? + Wag 2* + 23H Hz + 2m3,H.H, + uj. H,). 


Now, if we put here H,=m)£/ct, etc., then we do not obtain an 
expression representing the squared velocity multiplied by a 
certain factor. By a proper choice of the co-ordinate system the 
last expression can be transformed into 


3 (py H x? + Ho H i F [433 HL z’) 


and this can be interpreted as kinetic energy per unit volume, if 
it be assumed that the aether behaves as if it had different mass 
densities for motions in different directions. 

Let us still consider the electric energy which has to be 
interpreted as potential energy. This is, per unit volume of an 
isotropic medium, 


2DE =, = D= p, APm?(p? pP +e +r’). 


The aether has thus to be attributed the property that its 
potential energy is proportional to the square of the rotation 
of its particles. The question how we have to imagine such an 
aether will be taken up in the sequel. 

This theory resembles the light theory of Neumann, who also 
assumed that the aether density is the same in all media, and 
thence deduced that in polarised light the oscillations are in the 
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plane of polarisation. This means, with the present interpreta- 
tion, that the magnetic force is contained in the plane of polarisa- 
tion, as in fact has to be assumed in the electromagnetic theory 
of light. 


11. Tue Evecrric Force as VELOCITY. FRESNEL’S 
THEORY OF LIGHT 


As an alternative we might have interpreted kinematically 
the equations 
OH, OF, 10B: 
dy ie TT tt? 
This is entirely analogous to what precedes, the magnetic induc- 
tion being now represented by the rotation of the aether particles 
and the electric force by their velocity. The formulae now 
become 


ete. 


9 
E, =m": Ey =m 7 E,-= mi 

at ot , « (14) 
B,= —2cm'p, B, = -2cm’g, B,= - 2cm’r 


where m’ is again a constant. There is something strange in 
having to consider the electric force as a velocity. For an 
ordinary conductor carrying a constant charge we would now 
have a permanent, outward or inward, stream of aether. In 
the previous theory the continuous streaming of the aether 
near a magnet was at least circuitous, which state of things 
it was easy to imagine to last invariably for any length of time. 
The electric energy must now be correlated with the kinetic 
one. Its amount per unit volume, which is }eH?, has to be 


made equal to 
C ( acy? 
woi(a) +) +) 


so that p=ım’?. Since m’ is constant and e has for different 
dielectrics widely differing values, the density of the aether which 
is proportional to the specific inductive capacity must have very 
different values in different substances. As the previous theory 
approached that of Neumann, the present one resembles Fresnel’s 
theory. In the latter it is assumed that the displacements of 
the aether particles are perpendicular to the plane of polarisation. 
This means, with the present interpretation, that such also is the 
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orientation of the electric force. Thus on either interpretation 
we remain in harmony with what is generally assumed about 
the light vibrations in the electromagnetic theory. With regard 
to this coincidence of the model with Fresnel’s theory, we may 
still notice that to different specific inductive capacities corre- 
spond different propagation velocities of the electromagnetic 
disturbances, while Fresnel looked for the explanation of the 
different propagation velocities of light in the different densities 
of the aether within different media, and made this density 
proportional to the square of the refractive index, which amounts 
exactly to what we have assumed. With the interpretation 
treated in this section the electrically anisotropic media offer the 
same difficulties as did the magnetically anisotropic bodies in the 
previous one. 

For the interpretation of the second set of equations, (7), we 
consider the magnetic energy which has to be identified with the 
potential energy of the acther. It is a quadratic function of the 
rotation components p, q,r. Thus, as in the previous theory, we 
have to imagine a mechanism whose potential energy is just such 
a function. 


12. THEORY oF ELasticıty. MCCULLAGH’ S AETHER 


In an ordinary elastic body the relations are entirely 
different. Let E, n, č be again the displacement components 
(functions of the co-ordinates 7, y, z). Then we have, in addition 
to the displacement and rotation, the dilatations determined by 


0€ On oc 
ae u er (15) 
and the shears expressed by 


_ mh at ae am 
Ep ENT Fe 9) 


The normal tension components [stresses] are determined by 
X,=2Kir, +02, +9Y,+2,)}, Yy=2Kiy,+ Ma +Yyy +22); 
2,=2K 2, + (8, + Yy + 22)}, 

and the tangential components by 
Y,=2,=Ky. Z,=$,=Kz,X,=Y,=Kz, . (17) 
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where A and @ are the coefficients introduced by Kirchhoff. 
The energy per unit volume is 


K(22 + y,7 +27) + KO(te + Yy +22)? +4K (a,? + y7 +22"). (18) 


Now, for the aether this must be a quadratic function of the 
components of rotation, 
E M (OE Ay (a a 
P= ila a) 1-35, 2! [3 Ax Oy!” 

The next task would be (and such in fact is always the posi- 
tion in the older light theories) to make such assumptions as to 
give the energy that form. 

The most radical means for obtaining this result is to assume 
for the aether the validity of a very peculiar elasticity theory, such 
that the potential energy should in every case be a quadratic 
function of the angles of rotation. This is McCullagh’s aether. 
Let us see how this can be done. Elasticity cannot exist if there 
are no forces tending to bring back the displaced particles to 
their original position (state of equilibrium). What we are con- 
cerned with is the potential energy opposed to these forces. If 
that energy is to depend only upon the rotations, then the system 
of those counteracting forces has to consist only of couples. 
The components of the rotation being p, q, r, we have to assume 
for the components of the moment of the restoring couple, 
per unit volume, 


M:= - (ap +a127 + 497), My = — (ag1p + Aga + aor), 
M,= — (agp + 4999 + dagr), . . ° (19) 


where the a’s are constants satisfying the conditions a,.=dp, 
Ayg=Qgy, @z, =a. ‘These formulae hold also for anisotropic 
bodies. The peculiar feature of the case in hand is that this 
couple of forces should not be produced by the neighbouring 
aether particles. For such a couple would depend only on the 
relative rotation. 


13. QUASI-RIGID AETHER 


To explain these couples we may introduce a second medium 
which remains in its place, not sharing in the rotation of the first 
medium, and which exerts upon the latter those restitutive forces. 
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In the first medium we must then have also some inner forces 

which do not all vanish. To see this let us recall how in the 

ordinary cases the relations Y,=Z,, etc., are deduced. Let us 

Z consider a parallelepipedon dædydz and 

find the couple which produces a rotation 

about the X-axis (Fig. 11). The two forces 

parallel to the Z-axis, acting 

alt upon the faces dzdz, give the 

rose couple Z,dxdzdy, and those 

parallel to the Y-axis applied to 

O Y the faces dxdy yield — Y,dadydz ; 

Y Fia. 11. in all, (Z, - Y,)dzdydz. Now, this 

couple together with M,dzdydz 

must vanish, otherwise it would give an infinite angular 

acceleration, since it is of the third order, while the moment 

of inertia of the parallelepipedon is infinitely small of the fifth 
order. Therefore, for our aether, 

Z,- Yı= -M,. 
Thus there must exist tangential stresses. Similarly, of course, 
we find 


y 
Y, = X, SS M,. 
Moreover, we will make our assumptions as simple as possible. 


To this end we make the normal stress components zero, 
X,=Y,=Z,=0, and take for the tangential stresses 


—Z, = Y,=1M,, 
E X,= 2, =3M,, 
~Y,=Xy=4M, 


In order to see that these assumptions are in agreement with 
the adopted energy expression, it is enough to calculate the 
energy of the medium within a given surface for a stationary 
state. Let do be an element of this surface whose outward 
normal makes with the axes the angles \, », v. Let us imagine 
that the displacements £, 7, ¢ are produced by some external 
forces upon the surface elements, and that these forces as well 
as the displacements and stresses gradually mount to their final 
values. 
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The force acting upon the element do in the direction of the 
X-axis is 
(Xz cos À +X, cos p + X: cos v)do, 
and similarly for the Y- and Z-components. Since all the magni- 
tudes increase proportionally to each other, we find for the total 
work upon the element do 


EX, cosA+Ä,cosu +A,cosv)+n(Y, cos A+Y, cos 
+ Y, cos v)+¢(Z, cos A+Z, cos u +7, cos v)}do. 
A simple reasoning will show that the energy contained within 
the surface œ is given by this expression integrated over the 


surface. This integral can be transformed into the volume 
integral 


N o 
faé + Y 20 +Z) Fa lhe a Yan +Z,,b) 
+, (XE +Y j +Z) ldr, 


dr being the element of the enclosed saa Thus, the energy 
per unit volume is 


} Ha (Xb+ Yan + Zak) +5 (XE + Yon +0 


+ pg Hee + Yon+Z0)|, (20) 
or, by (11), (15), and (16), 
His, + Yyyy+ 2.2; +4(Y, +X )ty,+(¥,—-X,)r 
+3(Zy + Yo)yz + (Zu, - Yo)p +2(Az +2): + (Xz -Ze)g}. (21) 


In fact, we have 
an cE _ m, GE 
Yon Re X, +2, (7 ay) tC - x," ay 


etc., and in order to see that the remaining terms disappear, it 
is enough to remember that we are dealing with a stationary 
state. Thus, e.g., the coefficient of £ in (20) becomes 


res aX, XA arg 


de * Dy + de) Maya =D. 


Now, if the energy (21) is to depend on the rotations only, 
then in the first place the terms containing the dilatations must 
VOL. I D 
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vanish, 1.e. X,= Y,=Z,=0, and in order to get rid of the in- 
fluence of the shears, we must have Y,+Y,=0, Z,+Y,=0, 
X,+Z,=0. Thus we fall back upon the previous values of the 
stresses. 

The energy per unit volume now becomes 


= (Mp +M,q + M:r) 

= 3(a,1P° + 999” + agt? + 2a12Pq + 2agggr +2agırp), 
and it remains to be seen how in the two developed theories the 
coefficients a are correlated with the magnitudes appearing in 


the electromagnetic equations. ln the first case (Art. 10) we 
had 


H, PM H,= BL H, am | 


or’ ot’ ot (22) 
D,=2cmp, D, =2cmq, D, =2emr| 


and the aether density was p =m. 

The potential energy must coincide with the electric energy. 
The latter is, per unit volume, 

I(e',,D? +etc. +2e’,D,D, + etc.), 
since E, = €'ııD: + eiD; + €'15D.; etc. 
Thus the required coincidence calls for the following relations : 
Q12 = 4e’ ,.c?m?, Qo3 = 4e’,,c?m?, Qs1 = 4e’,,c*m?. 

For isotropic media we shall simply have E=e’D, a =4e'c?m?. 

The second set of equations (9) now follows from the equa- 
tions of motion of the aether, 


re 
ia t ay oz Pow O 


In fact, putting here 

X,=0, X,=1M, X,= -1M,, 

M, MN 0% 
Ipe 


dy dz / TP 


we have F 5, 


Now, 


M,= T m E 11 D, +e '»D, +e 132: = —2cmE,, etc. 
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Substituting this and replacing 0£/0t by H,/m, we find 
OF, OF, 10H, 


dy s c Ob’ 


which is the first equation of the second set (for the sub-case 
p=l). 

Similarly, in the second case, the set (7) of equations can be 
deduced from the equations of motion of the aether by con- 
sidering the magnetic energy. 

Attention to such an aether was drawn by Kelvin, who 
called it the quasi-rigid aether. With its aid he proposed to 
account for the magnetic phenomena, and considered, there- 
fore, the magnetic induction as the rotation of the aether elements, 
which coincides with our second case. 


14. QUASI-LABILE AETHER 


An explanation of the phenomena can also be arrived at 
without making about the clasticity of the aether such uncommon 
assumptions. Itwill be enough to consider the case of an isotropic 
homogeneous medium. Let us substitute in the equations of 
motion 


0X, OX, aX, P 
as * ay + Ge Page ben 


the stress components as given by the ordinary theory of elasticity, 


to wit, 
[oE OE dn AGN) 
A,=2K [ez ° +60 + Fay * +3) P 


i K($ + N), X, -K( + Aa 


Then the result will be 


[OE CE oE 0F On 0 Org 
Ki (sata re) +(1 +20), Ox ay i oe) Pop ae: 


We have now to make such assumptions that the left- 
hand members of these equations should become the rotation 
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components of a vector. This will be the case if we put 0 = - 
Thus, in fact, the last equation will become 


0 (OE An\ O/0C É 02g 
Kf? (3) -glos Hl- Pop’ 
or, in terms of the rotations, 


‚(or 0g RE 
Ky j 3) = EP apa’ 
which is the required form. 

With the interpretation of Art. 10 this equation is transformed 
into the first of the equations (9) for an isotropic medium, and 
with the interpretation of Art. 11 into the first of (7). 

The medium just considered is Kelvin’s quasi-labile acther. 
The equilibrium of this aether is labile, the potential energy in 
the state of equilibrium not being a minimum. In fact, it is 
not difficult to see that the potential energy can under the stated 
circumstances become negative. For its expression, per unit 
volume, is 


K(22 +y2 +22) + KO(&, +y, +2) +4K(2,?+y2 +22), 


and this, with 9= -1, is obviously negative provided that the 
shears vanish and z, y,, and z, have the same sign. 

This aether is perhaps not so satisfactory as the quasi-rigid 
one, because the corresponding theory is limited to isotropie and 
homogeneous media, so that also the case of boundary surfaces 
has to be left out of account. 


15. GRAETZ’S THEORY 


The last mechanical aether theory to be still considered is 
due to Graetz. In this the second set of equations is obtained 
in the same way as with the quasi-labile aether. In fact, what 
in the last case was achieved by putting 9= -1 in the equation 
of motion of an ordinary elastic medium, 


pas KIA +(1+ 2 


N ; OF On i 02 gals gc 
with P written for Ar + oo ta — |and A for zat Jp taal 


U MECHANICAL AETHER THEORIES 37 


Graetz obtains by adding to the right-hand member a term 
—2K(1 + 6)0P/0x, which converts it into 
oP Or o 
K(AE - a5) = -2E (ay 02) 

His theory then amounts to retaining for the free aether the 
ordinary elastic equations, so that in a vacuum Maxwell’s equa- 
tions do not hold. After all, since the free aether cannot be 
experimented with, we shall never be able to find out whether 
these equations do hold for it or not. For a ponderable medium 
Graetz assumes that the aether particles are acted upon not 
only by the surrounding aether but also by the ponderable 
substance. The force due to the latter gives then the term 


—2K(1 +6)0P/dx. 


Of this Graetz gives also some account, inasmuch as he deduces 
the said force from a pressure exerted by the ponderable matter 
upon the aether and represented by p=2K(1 +6)P. 


16. THE CHARGING OF A CONDUCTOR ACCORDING TO 
THE ELASTIC AETHER THEORIES 


We have just seen how it is possible, by ascribing to the 
aether various properties, to account for Maxwell’s equations. 
We must come back, however, to a certain difficulty which 
appears in all mechanical aether theories and which presents 
itself when we consider a charged conductor. If we look 
upon the electric force as the manifestation of a velocity in the 
aether, then, as was already mentioned, we must imagine a 
continuous aether stream towards or from the conductor, though 
we perceive in the latter no change whatever. 

How have we to picture to ourselves a charged conductor on 
the theory in which the magnetic force is represented by an 
aether stream ? The dielectric displacement is then a rotation 
about the lines of force in the sense corresponding to that of the 
lines. The components of this rotation are p=4(0¢/0y —0n/oz), 
etc., whence we see that the distribution of the rotation is always 
solenoidal. Let through every point of space a vector be drawn, 
indicating the direction of the rotation. Then the system of 
curves to which these vectors are tangential will be the so-called 
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vortex lines, and of these, owing to the solenoidal distribution, 
we can also construct vortex tuhes. 

Let us now consider a spherical conductor. During the 
process of charging, electricity is being communicated to it 
through a wire. This means that rotations are being produced 
in the aether, and therefore angular velocities exist while the 
conductor is being charged. Let us, then, consider the vortex 
lines in the whole space surrounding the sphere. In virtue of 
the solenoidal distribution these lines must either be re-entrant 
or extend from infinity to infinity. To each of these vortex lines 
must correspond another within the wire. One would now 
have to imagine that the aether fibres inside the wire are twisted 
and that this twist is propagated through the conductor and 
thence along the lines of force in the medium. The rotation in 
the latter is opposed by the elasticity of the medium giving rise 
to couples which, in absence of the rotating force, would at once 
untwist the system, that is to say, discharge the conductor. 
Now, the difficulty consists in finding out what happens when 
the wire is removed. For the state of affairs which actually 
takes place cannot exist in the picture just given. 

In order to see this, we consider an arbitrary surface cœ, 
whose normal at a point x, y, z has the direction angles a, ß, y, 
and which is bounded by a line s. Then, by Stokes’ theorem, 


je cos a +q cos 8 +7 cos y)do = A | (Ede + ndy + dz). (23) 


The field around the charged sphere after the removal of the 
wire is perfectly symmetrical. Now, if we take for o a portion 
of a concentric sphere, the surface integral represents, apart 
from a constant factor, the quantity of electricity which passed 
through this surface and is thus equal to the charge of the portion 
of the spherical conductor which is cut out by the cone subtended 
by s and having its vertex at the centre of the sphere. The same 
charge must also be represented by the line integral in (23). 
We begin with taking for o a small portion of the concentric 
sphere and we let this slowly increase. (The successive boundary 
lines s may, e.g., consist of a system of parallel circles.) Then 
the surface integral will continually increase, while this is 
impossible for the line integral after « has become greater than 
a hemisphere. Ultimately when the boundary s dwindles to a 
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point, the surface integral represents the charge of the whole 
sphere, while the line integral becomes nil. The point is that, as 
is shown by (23), it is for purely kinematical reasons impossible 
to have at each place on the surface of the sphere a rotation 
around the radius which, seen from outside, has everywhere the 
same sense. The difficulty arises from the fact that the rotation 
is throughout solenoidal, while the dielectric displacement has 
not this property wherever there are charges. 

Larmor, Reiff, and others tried to save the theory by giving 
up the symmetry around the sphere. They assumed that at the 
place where the wire originally was the conditions are somewhat 
different from the remainder of the sphere, viz. that there is a 
canal K at that place within which the aether is loosened from 
the surrounding aether, so that while the aether in K is fixed, 
the surrounding one acquires rotations corresponding to the 
electric force, with the result that the line integral of the aether 
displacement along a path embracing the canal is equal to the 
whole charge of the sphere. It should be possible to keep up 
such a state by applying to the aether outside K, all along 
the surface of the canal, appropriate tangential external forces. 
Instead of this, one might think of attaching the twisted aether 
to that within X, which would prevent a complete untwisting 
of the medium outside of K. No objection can now be derived 
from (23), since in this equation &, n, € are assumed to be 
throughout continuous, while this condition does not hold at the 
surface of K. 


III 
KELVIN’S MODEL OF THE AETHER 


17. KELVIN’S MODEL OF THE QUASI-RIGID AETHER 


KELVIN conceived a model of a quasi-rigid aether built up of 
gyrostats. This is a complicated problem. It amounts to 
finding a system which permits all deformations but resists 
such as are associated with rotations and no others. The 
idea occurred to him to meet this requirement by means of 
gyrostats, for these oppose themselves to any change of the 
direction of their axes. Thus his task was to find a system con- 
taining a number of lines which remain parallel to their original 
direction at every deformation devoid of rotation, and which 
change their direction as soon as rotations are produced in the 
system. Along these lines one had then to lay bars bearing 
gyrostats. 

Let us consider a homogeneous deformation in which, that is, 
the components of the displacement of any of the points of the 
system, &, 7, ¢, are linear functions of its co-ordinates z, y, z, 


n = Qs + AoT + (AgoY + (iga2, 


If the coefficients a,,, etc., are chosen arbitrarily, the deformation 
thus expressed is in general associated with a rotation. This 
can be found geometrically by recalling that a sphere is trans- 
formed into an ellipsoid and by determining those mutually 
perpendicular diameters of the sphere which correspond to the 
axes of the ellipsoid. The algebraic representation is found by 
rewriting the formulae thus : 
40 
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É =a; +011 +$ (a1 — Qg1)Y + 3(Q13 — G51)? 

+ 3(aj2 + 091)Y + 3(Q13 + Qg1)2, 
N =Gg + Aggl + (dog — A32)Z + F(ayı - 412)% 

+ }(dgg + Q39)2 + 4(Qyı + 012%, 


C =Gg + Ugg? + 4(Qgı — 413)% +3 (02 - aa3)Y ; 
+ (Ga +G13)% + 3 (Ago + Gog)Y¥- 


Thus the deformation is split into a shift, a dilatation, a 
rotation, and a shear, and we see that the conditions for an 
irrotational deformation are dj. =o), Ge3 =A39, 2313 =@3. In such 
a deformation, therefore, six coefficients, apart from the shifts, 
still remain undetermined. 

Now, to arrive at Kelvin’s model, we construct in a plane 
a system of congruent equilateral triangles fitting to each 
other and erect upon these 
triangles as bases, omitting 
every second (as shown in 
the figure), regular tetrahedra. 
The corners of these tetra- 
hedra lie again in a plane and 
form a system of points such 
as those of the ground-plane. 
On these points, therefore, 
such a system of tetrahedra 
can again be constructed. Fia. 12. 

This we do so that the bases 

of these tetrahedra should have the same position as those 
in the first layer of tetrahedra, so that the former can be 
obtained from the latter by a mere shift. That is to say, if 
the non-shaded triangles of the figure were first chosen as bases, 
we now take as bases the triangles PQR, RST, ete. In this 
way we can proceed, and similarly the system can be in- 
definitely continued on the other side of the ground-plane. 
Thus every corner point of the system will be the common vertex 
of four tetrahedra, at which also twelve edges will meet, two by 
two being prolongations of each other. Kelvin imagines now 
placed at every corner point a ball from which issue six bars and 
as many tubes, all of these being free to assume any direction 
whatever. The bars of one ball are now put into the tubes of 
other balls wherein they can be freely shifted back and forth. 
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This system can now be built up in the form of tetrahedra as 
just described. It will not oppose itself to any pure [irrotational] 
deformation. 

In the next place Kelvin introduces into each of our tetrahedra 
a system of three rigidly connected and mutually perpendicular 
bars of variable length (which may again be accomplished by the 
bar-and-tube method), and whose ends must remain in the grooves 
of the bars of the first system, so that they join each time two 
opposite edges of the tetrahedron. In a regular tetrahedron 
these bars coincide with the lines joining the mid-points of the 
opposite edges, but also in any tetrahedron whatever a set of 
mutually orthogonal intersecting lines, joining pairs of opposite 


A’ 





edges, can always be assigned.* The introduction of these 
systems of bars does not prevent any deformation. We will 
now prove that with an irrotational deformation of the original 
system the new bars are always shifted parallel to their initial 
direction. 

We first consider the case in which a regular tetrahedron 
undergoes an infinitesimal dilatation in the direction of one of 
the edges, say AC, which leaves the plane passing through BD 
and HE, the mid-point of AC, in its place. It is enough to 
prove that a set of mutually perpendicular joins of opposite 


* To see this, notice that through a given point one and only one line can 
always be drawn which joins two skew lines. If, therefore, P be an arbitrary 
point within a tetrahedron, there are through P three determined lines, each 
of which cuts a pair of opposite edges. By requiring these lines to be mutually 
perpendicular we have three equations for the co-ordinates of P. Their 
solution for the case of a tetrahedron which differs but infinitesimally from a 
regular one is implied in the following considerations of the text. 
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edges of the new tetrahedron can be found which are parallel to 
the joins of the mid-points of the original tetrahedron. For this 
will then be the only existing set. Now, leaving the line ZH 
unchanged, let us consider FG. If this is to remain parallel to 
its original direction, and still cut the line EH, its intersection 
point with the plane ADC must remain on the line EF and thus 
be shifted along it, say to F’. It is thus enough to show that 
the line through F’ parallel to FG cuts BC’ or, if G’ be the point 
where that line cuts ZG, that FF’ =GG'. Now, this follows from 
the congruence of the triangles DFF’ and BGG’. Similarly for 
the line joining the mid-points of AB and CD; in view of 
the symmetry, it is obvious that this line will again be shifted 
along EH, so that we find, in fact, three orthogonal joins of the 
opposite edges, parallel tothe originalones. What was just proved 
for a dilatation in the direction of the edge AC holds, of course, 
for the remaining edges, and our proposition will be completely 
established by showing that every irrotational deformation can 
be obtained by the superposition of six dilatations along the edges 
and of a displacement of the wholesystem. This, however, follows 
directly from what precedes, since an irrotational deformation is 
just determined by six independently prescribed ones. More- 
over, the relation between the values of the dilatations and the 
coefficients a,,, etc., in the general formulae can readily be 
established. In fact, if œ, 8,, yı be the direction cosines of 
the first edge, d, the dilatation in this direction, and similarly 
with changed suffixes for the remaining five edges, we have 
(apart from the shift) 


&=d,a,(xa, +yß, +2y1) +dyazlzaz +yBy +2yo)+ . - - 

+ dgag(Tag + YBe + 2Yo); 
n=d,ß,(xa, +yßı+2yı) + daßa( Laa + YBa +zya) +... 

+ deße(£ag + Ye + 2Yo): 
C= dy; (ra, + yBy + 2y;) + deye(tag + yBo+zye)+.. - 


+ dgye(tag + Bg + Ze); 


whence 


dB +da ++... +d U 
diy? tday +... +dgy_? Gy 
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daß: + daß, + e o oà + deaeh = Q19, 
Aßıyıtdaßayat+ : - - +deßeye = Gas 
yy 104 + dyVo0.+. . . +dgygtg =a), 


wherewith the d’s are determined, if the a’s be given.* 

If now at the middlemost bars an arrangement is made of 
such a kind that for rotating the bars a couple is necessary 
which, for an infinitesimal rotation, is proportional to the rotation 
and has its axis coinciding with the axis of the rotation, then 
there is still no force opposing an irrotational deformation of the 
first set of points. But if there is a rotation, a resisting couple 
is produced which is proportional to the rotation and whose 
axis will always concide with that of the rotation, provided 
the three middlemost bars have all the explained property. 
To see this we have only to consider one of the tetra- 
hedra. Let OP, OQ, OR be the directions of the inner bars, 
and let us consider an infinitesimal rotation @ whose axis makes 
the angles X, p, v with the bars. This can be resolved into a 
rotation ¢ cos X about OP, and ¢ cos u, & cos v about OQ and 
OR. In the first rotation OP remains in its place, but for the 
rotation of OQ a couple is necessary of moment Cd cos vr 
and with axis along OP. Similarly for the rotation of OR, 
making in all 2Cg¢ cos à. In quite the same way the rotation 
about OQ requires a couple 2C¢ cos u, directed along OQ, and 
that about OR a couple 2C& cosv. Compounding all these 
couples we have, as announced, the couple 2C¢ whose axis 
coincides with that of the rotation. 


18. SOLID GYROSTAT 


As a first device to make a couple necessary for changing the 
direction of a bar Kelvin proposed the ordinary gyrostat. The 
bar AB (Fig. 14) carries a fixed ring in which a second ring is 
mounted, free to rotate around PQ L AB. The diameter RS of 


* In fact, it can be shown that the determinant of the coefficients of 
d,, « « . ad, does not vanish, so that the equations are compatible with each 
other. In proving this the co-ordinate axes can bs chosen arbitrarily. If for 
these the joins of the mid-points of the opposite edges are taken, one of the 
direction cosines for each edge is zero, while the others become 1/4/2 with the 
same or opposite signs. The absolute value of the determinant is then found 
to be equal }. 
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the second ring carries at the centre O a fly-wheel with RS as axis. 
It is clear, first of all, that this arrangement does not oppose 
itself to a rotation of the bar AB around the axis PQ. Kelvin 
mounts then on every bar two such rings in mutually perpen- 
dicular planes. But let us first 
consider a single gyrostat. Its 


P 

inner ring can spin about PQ, 

and the fly-wheel about RS. fe 

For every contemplated rotation 4 w = 
Q 


and similarly for every moment 

with respect to an axis a certain 

sense will be assumed as the Fia. 14. 

positive one, viz. that from which 

the rotation appears to be anticlockwise. Let for the rota- 
tions in question the positive sense be given by OP and OR. 
Let w be the angular velocity of the fly-wheel with respect to 
the inner ring and let ZAOR=8, and, therefore, 9 the angular 
velocity about OP. Further, let Q be the moment of inertia 
of the fly-wheel corresponding to the axis RS, and Q’ that 
corresponding to the axis PQ. Suppose now that, while APB is 
kept fixed, the inner ring and the fly-wheel are spinning. Then 
Q' will be the moment of momentum of the fly-wheel along 
[about] OP, and its time-rate of change Q’9. This change of 
the moment of momentum can be resolved (Fig. 15) along OR 
and OT perpendicular to OR and OP, and therefore in the 
plane AOR. Instead of the angular 
velocity w and the moment of 
momentum Qw about the axis OR, 
the fly-wheel has after a time df 
the angular velocity w+dw and the 
moment of momentum Q(w+dw) 
about the axis OR’ (<ROR' =d@). 
Thus the components of the moment of momentum are, up to 
terms of the second order, Q(w+dw) along OR and Qwd9 along 
OT, and its rate of change Qw about OR and Qwé@ about 
OT. The couples produced by this change of the moment of 
momentum are due to the forces exerted by the inner ring on 
R and S which, however, can give no couple about OR. 
Thus, Qw =0, whence 





w =constant. 
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Let us now consider the system consisting of the inner ring and 
the fly-wheel. This is acted upon by forces applied at P and Q, 
which thus are unable to produce a couple about OP. If now 
q be the moment of inertia of the inner ring with respect to the 
axis OP, the moment of momentum of the whole system relatively 
to this axis is Q'Ò+gqĤ, and since this cannot vary, @ itself is 
constant. 

Let us now see what happens when the bar carrying the gyro- 
stat is turned from its position A,B, (Fig. 16) to AB contained in 
the plane A,P,B, and making with A,B, an infinitesimal angle e. 
We will suppose that the bar is being kept in its new position 
and that the exterior ring is held in the original plane A)P)B,, 
so that OP lies in this plane perpendicularly to AB. The forces 
which through this change of position are brought to act upon the 
inner ring have their points of 
application at P and Q and can 
thus produce no couple about 
PQ. Let there be initially no 
rotation about OP, and let 
RS fall initially into the direc- 
tion of A,B,(@=0). Then we 

Fia. 16. shall have only the moment 
of momentum about OA, 
amounting to Qw, which we resolve into Qw along OA and Qwe 
along OP. Now, suppose that owing to the said change of position 
the inner ring which originally did not spin acquires an angular 
velocity @. (This will turn out presently to have a non-vanishing 
value.) This angular velocity is shared also by the fly-wheel, but 
the moment of its momentum about the axis OR does not 
depend on the angular velocity @ but only on the angular velocity 
around OR, and the latter must thus still retain its original value 
w, because during the considered change of position the fly-wheel 
was acted upon only by such forces, at R and S, whose moment 
with respect to OR is nil. Thus, after the displacement, we have 
the moments of momenta Qw along OA and (Q' +q)@ along OP, 
and since there is no couple about OP, 


Qwe =(Q' +q), 


> Qwe 
Q= 5 
Q +q 





whence 
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Through the change of position of the bar the innermost ring is 
thus set spinning, the direction and the velocity of this spin 
being determined by the last formula. The subsequent course 
of things can be found by applying the considerations of the 
early part of this article (since these hold for every position of 
the bar AB). We can thus conclude that the acquired angular 
velocity 6 will remain without change of size, so that the angle 
0 will gradually mount to considerable values, and that the 
moment of momentum of the fly-wheel will undergo, per unit 
time, a change of which the component along OT, perpendicular 
to the plane of the inner ring, is Qw#. An equal couple is 
necessary to keep APB in the new position, and substituting 
for @ the value just found, this couple turns out to be 


Qw? 
+g 
and is directed along OT. Resolving it into components along 


OA and ON 1 OA (corresponding to 6 = 90°), we have, along OA, 


_ = i £ sin 6, 
and along ON, 
ee cos 6 
rg cos l. 

Tt is the latter component which is required for our theory. 
In fact, this couple whose direction coincides with that of the 
rotation from A,B, to AB is necessary to keep the bar in the 
new position, in other words, the bar resists this rotation 
with an equal couple of the opposite sense. The other com- 
ponent, proportional to sin 0, is due to the circumstance that 
some forces are also necessary to keep the plane ABPQ in its 
original position. To provide for this we can mount upon our axis 
yet another gyrostat in the same plane with and entirely similar 
to the first, with the only difference that its fly-wheel spins 
originally in the opposite sense. Then @ will have for the two gyro- 
stats always opposite values, and since the initial value of @ is 0 
for both, their angles @ themselves will be equal and of opposite 
signs. The couple about ON, necessary to reset the bar, will thus, 
of course, be twice as great. Similarly, the gyrostat which we 
have already introduced in the plane perpendicular to that of the 
first can be replaced by a set of two oppositely spinning gyrostats, 
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so that ultimately each of the bars carries four gyrostats. Then 
the system resists every motion associated with a change of the 
direction of the bars and no other, exactly as was required. 
The serious objection against this arrangement (as an aether 
model) is that the couple required 
to keep the bar in the new posi- 
tion becomes, owing to the factor 
cos 9, smaller and smaller (and 
finally even negative). 


19. LiQUID GYROSTAT 





Yet another device was pro- 

Fia. 17. posed by Kelvin to ensure that 

a couple should be necessary for 

changing the direction of the bars, viz. a liquid gyrostat. This 

consists of a ring-shaped tube (Fig. 17) filled with a circulating 

liquid and free to rotate about one of its diameters as axis. 

This axis lies in one of our previous bars. Such 

a single ring has the same effect as the previous 
solid gyrostat. 

The theory is much the same as before. Let 
the axis about which the ring rotates be chosen 
as X-axis and let the position of the ring be 
determined by the angle 0 contained between the 
normal of its plane, taken in 
the sense appropriate to the 
motion of the liquid, and the 
axis OZ. For 0=0 the ring 
would then have the position Z 
shown in Fig. 17, where the Zz Fia. 18. 
Z-axis is assumed to point 
forward. The actual position of the ring at any instant follows 
therefrom by a rotation 0 about OX. Compare Fig. 18, where 
OZ’ is the normal to the plane of the ring, while OY’ is 
contained in that plane. 

Let the angular velocity of the liquid be w and the moment 
of inertia of the liquid with respect to the axis of the ring Q, 
and that with respect to a diameter Q’ (so that, the tube being 
very slender, Q’ =3Q). 
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Consider first the liquid itself. This has at first the moment 
of momentum Q’6 about OX and Qw about OZ’, and then 
Q'(6+d6) about OX and Q(w+dw) about OZ’. The vectorial 
difference between the moment Q(w+dw) about OZ" and Qw 
about OZ’ can be resolved along OY’ and OZ’, giving as com- 
ponents of the rate of change of the moment of momentum Qe 
along OZ’ and -Qw along OY’. The liquid is subjected to no 
other forces than the pressure of the tube walls, and in virtue of 
a known property of surfaces of revolution that pressure gives 
rise to a system of forces all of which cut the axis OZ’ and thus 
can give no couple about OZ’. Consequently, Qw=0, and 
therefore w=constant. To consider now the system made up 
of the liquid and the tube, let the moment of inertia of the 
latter relatively to the axis about which it can spin be g. 
Then the rate of change of the moment of momentum of 
the system along OX will be (Q’+ q)6, and since the system 
is acted upon only by forces exerted by the axle, this is 
again zero, so that @ is constant. Thus, while the axis OX 
is kept fixed, the ring can spin about it uniformly and 
at the same time the liquid can circulate in the tube with 
constant velocity. But to keep the axis in position we must 
apply to it a couple -Qw about the line OY’, for, as we saw, 
the latter is the rate of change of the moment of momentum of 
the liquid. 

Suppose now that initially 0=0 and Ĥ=0, so that we have 
to do with the case represented in Fig. 17. No work, of course, 
is required to turn around the bar in the plane XOY, but a 
rotation about OY is opposed by a couple. To see this, we 
consider a rotation about OY through an infinitesimal angle e, 
which brings the bar OX (Fig. 19) into the position OX’ in which 
it is again kept fixed. The normal ON of the plane of the ring has 
now moved into the plane YOZ’ and will first coincide with OZ’ 
or at the utmost include with this direction an angle of the order 
of e; it has acquired, however, an angular velocity 6 about 
OX’ which will differ from zero though it can only be small 
of the order of e (since it is the effect of a change of position 
determined by e). In fact, after the readjustment we have the 
moment of momentum (Q’ +q)@ about OX’ and Qw’ about ON, 
and resolving both along OX and OZ and keeping in mind that 
initially 6 and @ vanished, we find for the increment along OX 

VOL. I E 
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the value (Q’ +¢)@+Qwe, and as for obvious reasons this must 
vanish, we have 

Qw 

0= - FT . 

Let us still notice that resolving the moment of momentum 
of the liquid after the rotation e along the axes OX, OY, OZ, 
the value of the last of these is found to be Qw’, since terms of 
the order e? can be neglected. And since the forces exerted on 

y the liquid during the change of the 
direction of the bar give no moment 
with respect to OZ, w’ must be 
equal to w. 

The value found for 6 is the 

, angular velocity acquired by 

the ring through the change 

7! Fic. 19. X of position of the bar, with 

which it continues to spin. 

For the subsequent motion holds all that was just said about 

the possible motions in the original position of the bar. If OY 

and OZ in Fig. 18 are taken to represent the directions denoted 

in Fig. 19 by OY and OZ’, it will become manifest that the bar 
must be acted upon by a couple about OY’ amounting to 









O 


5 O Qw? 
Qub 7 Q’ + q” 


2.2 
which can again be resolved into vn, cos 0 along OY and 
+ 


2,2 
esin 0 along OZ. ‘The rotation of the bar about OY 
requires thus first of all a couple about OY proportional to the 
rotation. But in addition to this another couple, about OZ, is 
required which can again be avoided by mounting upon the axle 
two gyrostats with liquids circulating in opposite senses. Then, 
of course, the required couple will again be doubled, t.e. amount to 


2,2 
20 +a cos 8. 


The difficulty pointed out in connection with the solid gyrostats 
exists also in the present case. These models can therefore be 
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used only if the bar carrying the gyrostat is not being displaced 
too long in the same direction, as e.g. in the case of periodical 
motions of small amplitude. 

We will consider the case in which e is a periodic function 
of the time (and remains very small), Then the same property 
must also hold for the couple determining the position of the 
bar. This vibration about the axis OY is then associated with 
an oscillation of the ring-shaped tube around the bar. In other 
words, also @ will be a periodic function of the time and the 
deviations will remain small. Moreover, since the position 
changes of the bar are very small, these oscillations of the tube 
can be considered as taking place about the axis OX (Fig. 19). 
Thus we have the following moments of momenta: (Q'+q)6 
about OX, (Q’+q)é about OY, and Qw about the axis ON 
normal to the plane of the ring. This axis changes continually 
its position and makes at a given instant the angles e, —@, 0 
with the axes OX, OY, OZ. In considering the moment Qw 
about ON, the latter axis cannot be replaced by OZ, for though 
these directions differ but verv little, yet the moment of momentum 
itself is not small. This, therefore, must still be resolved along 
OX, OY, and OZ. The component along OX is Qwe, that along 
OY, —Qw0, and that along OZ, Qw. Let us now suppose there 
were no other external forces than the couple producing the 
oscillations of the bar in the plane XOZ, which therefore has OY 
for its axis, and let the moment of this couple be K which, of 
course, is a periodie function of the time. We shall then have 
the following equations of motion: first, by considering the 
moment of momentum about OX, 


Que +(Q’ +q)0 =0, 


the zero on the right hand being conditioned by the vanishing 
initial values of e and @, and, second, by considering the moment 
of momentum about OY, 


-Quß +(Q’+g)e=K, 


while the third equation of motion expresses simply that is 
constant. Eliminating 6, 


Gr e+(Q +g)e=K 
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This has the form of the equation of forced vibrations. We can 
also speak of proper or free vib.ations of our system, which may 
exist when X=0. The frequency of these free vibrations is 


Qo 
Q’+q 


2r Q +q, 


w 


and their period 


In presence of an external force whose period is large 
compared with that of the free vibrations, the state at any 
instant will coincide with that in which the system would be in 
equilibrium under the action of the force prevailing at that 
instant. If the period of the external force is small, the phase 
difference between the vibrations of the system and of the force 
is $7. If the liquid is assumed to circulate in the tubes very 
rapidly, we may always limit ourselves to the first case. Mount- 
ing then on each of the inner bars four gyrostats, in a manner 
already explained, we can represent the aether quite well. 
This can also be so arranged (by varying Q or w) as to yield a 
non-homogeneous medium, by means of which also the refraction 
and reflection phenomena can be represented. For the free 
aether the period of the external couple K must then always be 
assumed very long as compared with that of the free vibrations, 
and thus also with the time of revolution of the liquid in the 
tubes. Otherwise the term with ë in our equation would come 
into prominence, and the aether would have to be given different 
densities for vibrations of different frequencies, which in turn 
would influence the velocity of propagation. 


20. Liguip IN TURBULENT MOTION AS AETHER MODEL 


Kelvin tried also to represent the aether by means of an 
incompressible liquid in turbulent motion. The dimensions 
characterising this medium are those of the vortices. Upon 
these a coarser motion can be superposed, as e.g. a pro- 
pagation along the Y-axis of transversal vibrations in the 
AOY-plane. If f(y,t) be the velocity of a particle due to 
this vibration [along the X-axis] and w’, v’, w the velocity 
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components of the turbulent motion at the point in question, 
the resultant velocity is u=f(y,t)+u’, v=v', w=w'’. The 
velocity corresponding to the vibration can be obtained separ- 
ately by averaging over spaces whose dimensions are small 
compared with the wave-length of the oscillatory motion 
but large compared with the vortices. The average velocity 
is then nil for the turbulent motion, while f(y,t) is nearly equal 
throughout the domain over which we have averaged. The 
averages thus defined will be expressed by bars over the letters 
representing the magnitudes in question. We will now deduce 
two equations which will show that transversal vibrations can, 
in fact, propagate themselves and will give us also their 
propagation velocity. The first of these is arrived at by con- 
sidering the momentum along the X-axis carried across a 
plane perpendicular to the Y-axis. In fact, such a plane is 
traversed by the liquid in either direction and even in equal 
amounts, but the liquid flowing in one direction can have 
a different momentum from that streaming in the opposite 
direction. For the excess of momentum transferred towards the 
positive over that in the negative direction we find, per unit 
area, p uw, if p be the density of the liquid.* By considering 
the increase of momentum within a short cylinder having its 
faces perpendicular to the Y-axis, we find as the first of the 
required equations 

Of(y,t)__ (u'v’) 

ee Ze Ze (24) 


Next, we have the usual equations of motion of a non-viscous 
liquid, 
ðu du du ap) 


a - +9 +W.: = - 
ot x öy % 0x 


Og | 
ae tM og tay ta, y, 
M he a un OP 
ot ox oy oz ð 


„=... (25) 


where p is the pressure divided by the density. 


* This can bo compared with the momentum transfer in a gas due to 
the passage of molecules through a plane, as considered in the theory of 
viscosity. 
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Substituting our special values, u=uw’+f(y,t), etc., the first 
two equations become 


of(y,t) (yt) ou’ -d fy ee 


ide 
au ar te a 


0x 
i u E to ‚ou z ‚ou op 
oy” oz tasj 
ov’ f a , v dp 
a = Fy Ne en Js” T az Toy f’ 


whence 
SUD u'r’) I es); nf 
a oT OF I; wil 


+0) +U 


f awe) ge) yl) yy iP| 
Tq“ On 02 Ox oy | 


Averaging both sides, the equation can be considerably simplified. 
Both f(y,t) and Of(y,t)/et can be considered as constant over the 
averaging space. Thus, 


lyst) _ Of (yst) 
at ae al =0, 


since v’=0. Again we have =0 for every magnitude ¢ which 
C 


fluctuates over distances small compared with the dimensions 
of the averaging space. We have still to take care of terms 
such as w’0(u’v')/Or. If this be written 


(wv) Ou (wv) 0w 
DE dL de oO 


2 


it will be seen that in the process of averaging we are left with 
the average of the expression 


wu (4 av’ 7) 
Ox + dy F az)’ 


but this is everywhere zero [the liquid being incompressible]. 
Thus our equation becomes 


ofwv') — —50fly,t) (op , op 
ae =—-V = _ (v dr +U 2). ‘ i (26) 
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Now, let pop be the pressurg when f(y,t) is throughout zero, 
and pop+p’p the actual pressure, and let us assume that 
f(y,t) is always very small in comparison with the velocities 
of the turbulent motion. In absence of vibrations the left-hand 
member of (26) and the first term of the right-hand member are 
zero, so that 

vf Po 4 y'Po=0. 

Ox oy 


We will assume that this equation holds also when w’ and v’ are 
taken to stand for the velocities of the turbulent motion accom- 
panying the vibratory motion. Then we can put 
‚0 Op Op ap 

Uy thay u +u iy” . oe (27) 
A weak point of this theory, however, is, among others, that one 
does not know whether the properties of the turbulent motion 
remain unaltered by the vibrations. But this and similar 
difficulties need not detain us, our only purpose being to describe 
Kelvin’s scheme in its main lines. 

In order to transform (27) we return to the equations 


Ou pu" ot woe = _ Op} 
t xr Oy z ĉx 





w ‘ w BR = w _0p 
ab an “dy * dz ~ oy | 
aa ae ay a2 ~~ az J 
from which follows : 
Aa Í o a er z =): _ 2( Ow Ow du du >) 
P= -=| ay 02 02 dy * Ox 02 "Oy Ox" 
If now u, v, w are replaced by the values they have in the presence 


of vibrations, f(y,¢) affects only the term with Qu/dy, and thus 
we find 


4 ew (25) 





1 g f(y, t) ov’ 
Ap = -2 y in" 


The solution of this equation can be written symbolically 


Pag _0f (y,t) Ov 
ER WING) OU 
p = -2A y Os 
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The meaning of the symbol A"! ıs known from the theory of 
Poisson’s equation Ab =x, whose solution is 


1 fx 
hehe ia 
d=A „far. 


In our case, owing to the rapid changes of sign of 0v’/0x, while 
Of(y,t)/ay varies much slower, only the nearest neighbourhood of 
the point for which p’ is to be evaluated contributes to the space 
integral. Consequently, ¢f(y,t)/¢y can be regarded as a constant 


factor, so that 
r fY) e, 
p= oy a: On 
Putting 
Op" OP" 
vg oy =Q, 


and keeping in mind that Of(y,t)/0t varies slowly, we find 


z eu t) u 0 1! 
Q= (oz da” wa JA Ou 


fl 2 “2 \ 
E rw Aw... (28) 


oy dx DxO. 


Now, by averaging, 
Oat Pag P, 

and, though this is not quite correct, we may assume that also 
the third mean, 


v? 
go w, 


has the same value. Thus, 





vA A-W = po AA -I hot 


To find the last term of (28), we put a-w =. Then 
y 


wt gee) 08 


z os Ox" 
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In averaging, the first term contributes nothing, so that 


„t ‚or (ou. Ow" 


Oia og Ya +) = T 
ö- UN ay FR ou a nl. 


and by reductions similar to those hia made, 


y of(y,t) , n, 
gas Aa (2-3 vA | 
Au ya AA, 
oy J 
. a a 
G= -pU 


If R? be the mean square of the velocity in the turbulent motion, 
w’? +v? +w? = 3v? = R?, and (26) becomes 





oluv’) RW, 0200.99) 
a oy 


This is the second of the required equations. Combining it with 
our first-found equation, 


Au’v‘) _ _ fyt) 
dy “Oe (30) 
let us eliminate u’v’. This gives 
a Of 
zk? 
oy” 


which shows that in our liquid transversal vibrations can be 
propagated, with the velocity “? R, 


IV 


ATTRACTION AND REPULSION OF PULSATING 
SPHERES 


21. NATURE OF THE PROBLEM 


CoNSIDER an incompressible frictionless liquid in which spheres 
are moving about. This motion may be either a translation or 
a swelling and shrinking of the spheres, both of which may be 
periodic. 

We will first take up the problem of finding the motion of 
the liquid if that of the spheres be given. 

Particles of a non-viscous incompressible liquid which once 
do not rotate will never rotate. There exists then a velocity 
potential, and the velocity components can be represented by 


oo am 2 


“en o= y’ nn: 


Owing to the incompressibility P will satisfy Laplace’s equation 
A®=0. Whence it follows that the velocity at every point of 
a space-region will be determined, if at every point of its boundary 
the normal velocity is known. For then 0@/en all over the 
boundary is known, and thus also ® throughout that space- 
region is determined, since it is continuous and satisfies Laplace’s 
equation. Let us consider the space limited by a fixed closed 
surface at infinite distance and by the surfaces of all the spheres. 
Then the normal velocity all over these boundaries is known, 
the motion of the spheres being given. If a, b, c be the co- 
ordinates of the centre of one of the spheres and R its radius, 
the motion is determined by å, b, é, and R. It remains to find 
a suitable solution of Ab=0. This will be supplied by 


1 ak 
Seth = TAE cok, ites Page Me 
Peal) an. a 
58 
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where A, hy, ... h, are measured along arbitrarily chosen 
directions [and r is the distance from the centre of the sphere]. 
This is a homogeneous function of the space co-ordinates of 
degree -(k+1). If we put 


V= Yr- t0, 


then (the spherical harmonic function) Y, is a homogencous 
function of the zeroth degree. Next, write 


H; = V;r*! 1 = Yırk; 


then H, will be a homogencous function of the %* degree 
which can be represented by a whole algebraic function of the 
co-ordinates. For at every differentiation the highest exponent 
of r in the denominator is increased by 2, so that the highest 
exponent in V, is 2k +1, and the remaining ones are by an even 
number lower. Thus in H, all denominators fall out, while in 
the numerators only even powers of r occur. H, also satisfies 
the equation of Laplace, but it cannot be used for representing 
the velocity potential in a liquid extending to infinity, since it 
becomes infinite at an infinite distance. On the other hand, V, 
cannot be used if the origin of r lies within the liquid. ' 


22. A SINGLE Movina SPHERE 


Consider the space outside a single sphere. Then the simplest 
solutions are, for k=0, Vy=I1/r, Yo=1, H)=1. If we take 
w=C/r, the velocity will be radial and equal to —C/r?. By an 
appropriate choice of C as a function of the time this solution 
can be adapted to a sphere with fixed centre and variable radius R. 
In fact, since at the surface R = — C/R?, we have 
RR? 

r 
For k=1, the next simple case, we find, with the X-axis along 
the direction of A, 
0 /1\ x-a x-a 
hs -h7 r3? He r’ u. 
If we take for the velocity potential outside the sphere 


p- -4)_ Ce-a 1 
o BPR r Č 
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then, at its surface, 


m U RR BR 
This can be adapted to the case in which the sphere has a 
translation velocity along the X-axis. For then the component 
of this velocity along the surface normal is also proportional to 
cos A. Thus, 


om PT 1 _ _20cosA 


-20 cosA A 


“RS =d cos A, 
whence C= -jaR3, 
and ®--WRTZ. . 7 |. 81) 


Since (x —q)/7r3 arises from the differentiation of 1/r, this motion 
of the liquid can also be obtained by a superposition of the 
motions due to two spheres devoid of translation but pulsating 
in opposite phases and placed at an infinitesimal distance from 
each other. 

Let now the liquid be in a certain state of motion and let 
us ask how this motion is disturbed by the presence of a sphere. 
It will appear that due to the latter a certain state of motion is 
superposed over that already existing. We will call this the 
reflected motion. It will be the weaker, the farther away from 
the sphere. 


23. A SPHERE AT REST IN A LIQUID WITH A GIVEN 
MorTIoN 


Let us consider the case in which the disturbing sphere is at 
rest. We choose its centre as the origin of co-ordinates. The 
[irrotational] motion of the liquid being given, so also is the 
velocity potential & which can be developed around the origin 
into 

ee) 


where the values of the derivatives of ® prevailing at the centre 
have to be taken. The state of motion thus expressed can be 
interpreted as if it were due to a superposition of different 
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motions whose velocity potentials are represented by the succes- 

sive groups of terms. Thus, a velocity potential &, represents 
rest, the next group, 

(298.29) 

ae I iy dz 

expresses a constant flow with the velocity prevailing at the 

centre, and so on. This series development can also be put 


into the form 
®=-LChrrY,. 
k 


The reflected motion, combined with the original one, must give 
all over the surface of the sphere the radial velocity 0. Let us 
put for the reflected motion 


o-xe, % 
a kt? 


and let us suppose that to each of the motions to be compounded 
there corresponds a reflected motion. Then we must have at 
the surface, for every k separately, 
Y 
CAR -1Y,- Cr (k + l) pets =Q. 


Consequently, 


k 


s, RB 2k-4-1 
Ci + ¿CR , 


k 
and the velocity potential of the reflected motion becomes 


k Y 
D =E; OREH par 


24, n SPHERES 


Let us now turn to the problem of finding the state of motion 
due to the arbitrarily prescribed motion of n spheres, of which 
the positions, that is, as well as the sizes and velocities, are given 
for all times. The velocity potential @ must then satisfy 
Laplace’s equation and the condition that at each of the surfaces 
ocp/On should have given values. We shall find ® by super- 
posing n solutions, each corresponding to the motion of a single 
sphere. Each of these solutions will be found by proceeding 
first as if the sphere in question existed alone and then by taking 
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account, in an obvious way, of the repeated reflections of the 
state of motion produced by that sphere. Suppose, for instance, 
there were only two spheres, the first of which has a given 
motion, while the second is at rest. Let, in terms of the 
velocity potential, ®, be the motion due to the presence of 
the first sphere alone, ®,; the motion arising through the 
reflection of , at the second sphere, ®,,;; the motion produced 
by the reflection of ®,, at the first sphere supposed to be at rest 
in its instantaneous position, and soon. Then the actual motion 
will be found to be represented by the infinite series 


© =O; -+ Oir + Di + BD Sans 


which will converge the more rapidly the smaller the spheres in 
comparison with their mutual distance. This series gives the 
solution of the problem. For the equation of Laplace is satisfied, 
since each of the terms does satisfy it and the derivatives of 
the series can be found by differentiating separately its terms, 
and all boundary conditions are satisfied. In fact, at the surface 
of either sphere Oc/cn is expressed by a series whose terms, as 
far as the second sphere is concerned, cancel each other in pairs, 
while for the first sphere the first term only survives and this 
has the value which is prescribed for its surface. 


25. Two SPHERES 


We will now determine the forces which two pulsating spheres 
exert upon each other. The equations of motion 


op (eu Ou Ou Ou 
~ Ox eat tn tun.) Re 


become, in terms of the velocity potential, 


Op _ OD ,[ (00)? (a0)? /a® N] 
Ox -Por at +) + (2) (oe) f pete 
and since p is constant, these give, with V written for the velocity, 
oD 
p= -p(y +402) +0, > a a (82) 


where C is a constant. This, however, can be omitted, since a 
pressure which is constant all over a surface cannot produce 
any motion. 
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Let R, and R, be the radii of the spheres and J the distance 
of their centres. Then the velocity potential due to the motion 
of the first sphere is, by Art. 22, 


where r’ is the distance from the centre of this sphere. 

If now the origin of co-ordinates be placed into the centre of 
the second sphere and if the X-axis be laid along the produced 
join of the centres, the last expression can be developed around 
the second sphere into the series 


D= - Reh, 5-5 + a N: 
of which it will be enough to retain the first two terms, as 
according to our assumption the radii of the spheres are small 
compared with l. Since the first term is constant, the motion is 
determined by the second term. At the second sphere a reflected 
state of motion is produced, for which the velocity potential is, 
by Art. 23, 

IR, Rey ols? oy 
where r denotes the distance from the centre of the second sphere. 
Rigorously we should take account also of the motion arising 
from this through a reflection at the first sphere, and so on. 
But since these reflected motions become rapidly weaker, they 
may be disregarded. The velocity potential due to the motion 
of the second sphere is - R2R,/r which, to the desired degree 
of accuracy, need not be supplemented by any reflected motion 
at all. For our purpose is to determine the force experienced 
by the second sphere. This, however, requires only the know- 
ledge of the motion of the liquid in the neighbourhood of the 
second sphere, and the reflection of the last-mentioned motion 
from the first sphere would contribute only an expression con- 
taining a factor of the order 1/l?. Let us still put 


1R? =a, 42,8 =a, ; 


then the whole velocity potential in the neighbourhood of the 
second sphere will be 
_ ol Q5t\da, 1 da, 
an (1-7 ins 2) rd 
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By (32), to determine the pressure, 0/0t is required. Now, 
since in differentiating with respect to the time z, etc., are kept 
constant, 

Oo l (1 > x ( d*a, da, a 1 d?ag 


PD _ _1/1_2\Pa „a ( Ca da) _1 d'as 
AT TIVTU Æ PPAR t u dr 


To find the force upon the second sphere, we have to integrate 
over its whole surface, but in doing so we may disregard the 
term }V? in (32).* Again, all terms not containing z can be 
omitted, since these represent a pressure uniform all over the 
surface which, therefore, as was just mentioned, gives rise to no 
force upon the sphere as a whole. We are thus left with 

BR „N Pay 3 (a da, da, daa) 

PE PP) Ia de® * QRS"? dee” de dt 

oe zl? Pa, 1 da da), 

= PN) 32 di? "2a, de dtf 
Now, a pressure p = -crp gives for the X-component of the force 


upon an element do 
+cpR, cos? bdo, 


whence, the force upon the whole sphere, 
Jeph, cos* Odo = incpR,? = Ancpag, 
and substituting the value of c, 


p(, a, ‚Ta, da, 
dnb $0, die +? dt = 


We will suppose that the spheres pulsate rapidly, so that the 
perceptible force is the average of this expression over a time 
interval comprising many periods. Moreover, let both spheres 
have the same pulsation period. Total derivatives with respect 
to the time can then be omitted, so that the last expression 
for the force upon the second sphere can be written 


Amp| „d (a _da,da)) .. | , (8) 
2 \2dt\2de/” dt df 

* For, evidently, we can limit ourselves to that part of $V? which 
depends on the motion of the first as well as on that of the second sphere, 
that is to say, to the scalar product of the velocity-vector V,, which would 
correspond to the motion of the second sphere alone, into J’,, due to that of 
the first sphere only. Now, at the surface of the second sphere V, is radial 
and V, tangential, so that their scalar product vanishes. 
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and the required average ultimately becomes 


Whence we see that, if the pulsations of the two spheres are 
in phase, this force is negative, that is to say, we have an 
attraction, and if their phases are opposite, a repulsion. Let us 
consider this in more detail for the case of simple harmonic 
oscillations. If the radius varies periodically, so also does the 
volume of the sphere, and with the same period, and if the 
changes remain small, a simple harmonic variation of the radius 
is associated with a simple harmonic variation of the volume. 
We can write, therefore, 


a, = A, +c, cos (nt +e), 
where c, is small compared with A,, and similarly, 

Ay = Ay + Cy COs (nt +€,). 
This gives for the averaged force along the X-axis 

Ko = - “Poop? sin (F é) sin (nf +6). 
In order to find the mean over a long time, it is enough to 
average over a full period. Now, since 
2 sin (nt +) sin (nf +€,) = — cos (2nt + €, +€,) +08 (€, —€,), 

the required average turns out to be 


K,=- ae CyCyn” cos (€, — €a). 


Thus we have an attraction for e, — ¢, < 7/2, reaching a maximum 


for e; =€z, and a repulsion for e; —¢, > 7/2, attaining a maximum 
for ej -e= T. 


26. TREATMENT OF THE PROBLEM BY MEANS OF 
LAGRANGE’S EQUATIONS 


The results just found might also have been deduced from 
Lagrange’s equations of motion. In using these we will introduce 
VOL. I F 
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as co-ordinates, also for the case of a system of any number of 
spheres, the Cartesian co-ordinates @, b, c of the centres and the 
radii R of the spheres. The kinetic energy T is then a quadratic 
function of å, b, é, R with coefficients which are functions of the 
co-ordinates. If X be the component of the external force upon 
one of the spheres, taken along the X-axis, we have 


aoa) a 


The kinetic energy 7 can be split into two parts, one due to 
the liquid, the other due to the spheres, and in accordance with 
this the inner force keeping equilibrium to the force X can be 
divided into two parts. We then find for the force exerted upon 
the sphere by the liquid along the X-axis 


Ti... an 


if T, be the kinetic energy of the liquid. The latter is (owing 
to A =0) 


rb (C2) + 2) je- -ifa 


the last integral to be extended over the surface of all the spheres 
(it being assumed that we need not reckon with the infinitely 
distant boundary), and the normal to be taken towards the liquid, 
i.e. away from the spheres (hence the negative sign). The last 
integral can, of course, be divided into parts, each extended 
over the surface of one of the spheres, thus : 


[00 
T,= -p2 [0 do. 


In the case of two pulsating spheres we have for the second 
sphere OoP/On=R,, and therefore, 


fonao- [1 al -7- aim], are 


where a,, a, are as explained in the preceding article. On in- 
tegrating, the second and the third terms contribute nothing, 
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while the remaining ones are constant, and we are thus left 
== 1,da, , 1 „da 

-4r Rè( RB, + BR, ,?). 

= Ar 2 dt ` R, ade) 


Similarly the integral over the first sphere is 


lda, 1 da 
es 2( 2 aA 
oR a: +p Evii ) 
The symbol a appearing in (34) stands in this case for our l, and 
the first term of that expression gives 


oT, _tmp da, dus, 
l P dt dt 


This is exactly the force we have already found. With regard 
to the second term in (34), notice that what we require are the 
forces which act on the spheres when these have only a pulsating 
motion, while their centres are fixed. ThusÎ=0. In evaluating 
the first term of (34) we can, for obvious reasons, simplify T, by 
putting /=0, as we did in fact. In the second term, however, 
this simplification can be introduced only after the differentiation 
with respect to l. The evaluation of this term, which would 
lead us to the first term of (33), would thus require a further 
consideration of T7,, into which, however, we need scarcely 
enter, since we are concerned only with the time derivative of 
oT joa, and the average of such a derivative over a full period 
is nil. 


27. PEARSON’S THEORY 


Consider the case of an infinitely long period [of pulsation], 
that is to say, a sphere which goes on expanding for ever. Then a 


. , d /4rR,3 
surface enclosing the sphere is traversed by the amount am 


or Anda,/dt of fluid. Let this be denoted by e, Then the 
attraction between two such spheres will be pe,e,/Arl?. 

Now, Pearson abolishes the spheres and supposes that there is 
towards certain points an incessant stream of fluid (aether sources 
or sinks). It is true that one cannot well picture to himself such 
a state of things. But as we can imagine points spread over 
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a surface, towards which the aether streams from the sur- 
rounding space, to spread itself then over the surface, we 
should also be able to place such aether sources in three- 
dimensional space by calling to our aid the four-dimensional 
space. Similarly, there would exist points (aether sinks) where 
the fluid is being annihilated.* Two such aether sources would 
then attract each other, and similarly two aether sinks, while a 
source and a sink would repel one another. 

Apart from the strangeness of such a representation there is 
also another objection. In fact, if the expanding sphere be 
entirely omitted, there would still be a force upon the place 
whence the aether emanates. If a small sphere be described 
around the aether source, the force exerted by the aether will 
be a pressure upon the sphere, and one would have to imagine 
that the source is being displaced together with the sphere. 

There is a striking difference between the present case and 
that of electrical actions. For here we have attraction between 
points [sources] of equal, and repulsion between points of opposite 
signs. 

A similar theory was proposed by Korn to account for 
molecular forces and gravitation. He imagines a number of 
pulsating spheres, all in phase with each other; the latter 
coincidence is secured by enclosing the whole space in a limiting 
surface which is acted upon by a periodical external force. This 
is propagated instantaneously through the incompressible fluid 
and makes the volumes of all the spheres alternately increase 
and diminish in the same phase. 


In what precedes a description was given of some of the 
attempts which were made in order to account for various 
phenomena, and especially the electromagnetic ones, by means 
of speculations about the structure and the properties of the 
aether. Toa certain extent these theories were successful, but 
it must be admitted that they give but little satisfaction. For 
they become more and more artificial the more cases are required 
to be explained in detail. Of late the mechanical explanations 
of what is going on in the aether were, in fact, driven more and 
more to the background. For many physicists the essential 


* [For the three-dimensional beings, that is.] 
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part of a theory consists in an exact, quantitative description of 
phenomena, such e.g. as is given us by Maxwell’s equations. 

But even if one adheres to this point of view, the mechanical 
analogies retain some of their value. They can aid us in thinking 
about the phenomena, and may suggest some ideas for new 
investigations. 
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KINETICAL PROBLEMS 
(1911-1912) 


INTRODUCTION 


In these lectures some self-contained questions concerning kinetic 
theories are treated. They belong partly to the domain of the 
kinetic theory of gases and partly to that of the electron theory. 
Their subject was suggested by Knudsen’s investigations on very 
rarefied gases and by Richardson’s researches on thermionic 
currents. 

The gases offer two extreme cases which can be treated with 
comparative ease, one, in which the dimensions of the containing 
vessel are very large, and another in which these are very small 
compared with the mean free path of the molecules. 

In the former case there will be no sliding of the gas along a 
solid wall and this will have the same temperature as the con- 
tiguous gas layer, while in the latter case the gas is so rarefied 
that the molecular collisions can be disregarded. A volume- 
element does then no more contain the same matter during 
a certain time, as can be assumed to be the case for gases of 
large density. The investigation of cases falling between these 
two extremes offers considerable difficulties. In dealing with 
Knudsen’s investigations one can start from the second 
extreme case, and this is the plan which will be here adopted. 
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CHAPTER I 


INNER FRICTION AND SLIDING, TREATED HYDRODYNAMICALLY 


1. HypRODYNAMICAL EQUATIONS OF AN INCOMPRESSIBLE 
Viscous Liquip 


To begin with, we consider an incompressible viscous liquid and 
write down its hydrodynamical equations. In addition to inner 
friction the sliding along a fixed wall will also be taken into 
account.* 

Let u, v, w be the velocity components of the liquid, p the 
pressure, X,, X,, etc., the inner stress, including the pressure, 
further, u the viscosity coefficient and, finally, p the density. 

For the sake of clearness it may be mentioned that X, is the 
X-component of the tension exerted upon a surface-element 
normal to the X-axis, Y, the Y-component of the stress upon a 
surface-element perpendicular to the Z-axis, and so on. 

In absence of external forces the equations of motion are 





r +% m0 (equation of continuity), . . (1) 
(tut 4 ot = _0X,, 0X, 0X, 
Pay) Oa” Oy * oz 

Gv ov ov. dav, OY, OY, OY, ; 
(ag tae + a +) = a tyta fp © 





OW es Oe we) Ze 02, OL; 
(a Ox Oy Gz) Ox ty toz 


* Helmholtz and Piotrowski investigated whether there is sliding of a liquid 
or not. They observed the oscillations of a hollow metallic sphere filled with 
liquid and suspended on a twisted wire, and they found that the sliding at the 
metal wall was not nil. Wiener Sitzungsber. xl, Abt. I., 1860, p. 607 
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X, -p+ ui 


Y= -p +243; |» 2020.08) 


Ow 
Z: = -p + 2pm 
Ou av 
Xy=Ya=p 5 +n) 


eae Hon tay) nO) 


nee 


* In deriving these formulae one considers a parallelepipedon of edges dz, 
dy, dz within the liquid. Multiplying both sides of each of the equations (2) 
by dadydz, we have on the left hand the product of an element of mass into 
its acceleration, and on the right hand the force acting upon this liquid mass. 

With regard to the equations (3) and (4) we may notice that in absence of 
friction X,= 1',=2Z,= —p, while the tangential stress components are all nil. 

In the presence of friction all stress components will be determined by ex- 
pressions depending on the manner how the velocity varies from point to point. 
The equations representing these connections will be linear ; to a first approxi- 
mation the stress components will be determined by the derivatives of the 
velocities with respect to the co-ordinates. 

We can thus write, in general, 





cu cu Cw cw 
cu E wW 
and so on. 


To determine further these coefficients we take account of the symmetry 
relations. Imagine the liquid reflected at the plane x=0 and notice that for 
the image the same equations must hold. Tt will then become plain that terms 
having an odd number of references to the X-axis change their sign through 
the reflection, while those with an even number of such references retain their 
sign, The reference to the X-axis can occur in a fourfold way, to wit, as in 
ox, Cu, X, and Z,. 

Thus, using a reflection at the plane x=0, we shall find that X, remains 
unchanged, and this gives Qiz =&äy -a13 =la = 0. Again, Y, remains what it 
was, giving 5,. = be, = 043 = 53, =0. 

Next, a reflection at the plane y = 0 leaves X,, unchanged, whence a3 = 43, = 0, 
while Y, changes its sign, and this gives 6,, =6,,=6,,=0. 

From the fact that the exchange of the Y- with the Z-direction can have 
no effect upon X,, it follows that a, =a,,. 

The formula for X, thus becomes 


X=- - Pt Oye ta (+), 


oy 
which, in view of the equation of donitinaiiy, can be written 


= — p + (,; — Gy) a 


Taking into account that the couple resulting from the tensions upon an 
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With regard to the boundary conditions we may notice that, 
if there be no sliding, the velocity of the liquid in contact with a 
solid body is the same as the velocity of the latter, but if there is- 
sliding, this equality holds only for the normal velocity com- 
ponents. In the latter case, therefore, the relative velocity of 
the outermost layer of the liquid and the solid wall has a tangential 
direction. 

Let us now consider within the liquid near the solid wall a 
short cylinder whose dimensions along the normal of the wall 
are infinitely small compared with those in tangential directions, 
and let us express the condition that the forces exerted by the 
liquid upon this volume-element are in equilibrium with those 
due to the solid body. 

Let A be an arbitrary direction in the tangential plane. Then 
the force exerted by the liquid, per unit area, in the direction h 
may be represented by H,,, and that exerted by the solid body by 
AU, Where v, is the relative velocity of the liquid and the solid 
wall and X a proportionality factor. Thus the boundary condition 
will become H, =No. . (5) 


If there be no sliding, then v, =0 and we must put A\=». 

Let the state of motion be stationary, so that du/at, cv/ot, 
dw/dt all vanish, and let us further assume that the velocities are 
so small that their products and their derivatives can be neglected. 
Then, in virtue of (1), the equations of motion (2) become 





op \ 

ae +pAu=0 
0 

a +pAv =0}. ee. (6) 
Op 

=, rn Aw=0 


ab us Eh nat yes E use u a en 
element of liquid must be zero, if we limit ourselves to magnitudes of the third 
order, we find Y,=Z,, whence b,,=b,,. Thus, 
v 
Ya=bas toy) 

Again, since the liquid is isotropic, the expressions for Y, and Z, must 
follow from that for X, by a cyclic permutation, and similarly those for Z, and 
X, from that for Y, Whence it follows that the coefficients in the formulae 
for X,, Y,, and Z, must be equal, and similarly for Y, Zy and X, 

The isotropy of the liquid implies also that if the equations be transformed 
to a new system of axes obtained by a rotation of the original one, the co- 
efficients retain their values. This gives a1, — G33 = 2bgs = 2u. 
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2. EFFECT OF SLIDING UPON A LIQUID FLOWING IN A TUBE 


As a first application we consider a liquid flowing through a 
narrow tube. We put the X-axis along the tube. 

The equations of motion can be satisfied by putting v =w =0. 
This amounts to disregarding small lateral motions near the ends. 
The equation of continuity (1) then calls for du/dx=0, so that 
u and therefore also Aw become independent of x, and, by the first 
of (6), the same is true of Op/dt. By the remaining two equa- 
tions (6) the pressure must also be the same all over the cross- 
section. This will also be the case of dp/dz. Consequently the 
pressure varies uniformly along the tube and depends but on a 
single co-ordinate, x. If p, and p, be the values of p at the 
beginning and the end of the tube, and if I be its length, we have 


0p _Pa = Pı 
Ox l 
and, by (6), 
Ou + Ou _ Pa - M, 
oo ul 


Limiting ourselves to a tube of circular section, of radius R, 
transforming to polar co-ordinates, and noticing that, since the 
motion can be assumed to be axially symmetrical, u depends 
only on r, we can write the last equation 


du l du _p- Pr, 
dr? > dr pl” 
whence, by integration, 


„du _1P- Pip 
"dr 2 pl 


an additive integration constant being omitted, since the velocity 
of flow is a maximum at the middle of the tube, so that du/dr = 
forr=0. Integrating once more, we find 





u= ore +C), 


where C is to be determined from the boundary condition. In 
absence of sliding u =0 for r= R, whence C = — R?, and therefore, 


u Dem, 
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The amount (volume) of liquid streaming through a cross- 


R 
section is, per second, 27 | urdr. Thus, in our case, 
0 


Daph e 


which is proportional to #*, in agreement with Poiseuille’s law. 

In the presence of sliding the friction force -— udu/dr is equal 
and opposite to the force exerted upon the liquid by the tube 
walls. The boundary condition then becomes 


du 
-Pg =Au, for r=R, 
whence 


9 
—__ p_e 
C= -R yh 


and 
_ Pi — Px pe 2p y2 


The corresponding amount of liquid streaming across the section 
of the tube is, per second, 


Tal Rp: - Pa) . « « « (8) 


The term 4u/AR gives the correction to Poiseuille’s law for the 
sliding. For ) =o, or no sliding, (8) reduces to (7). 


3. Toe Draaaina oF A Liguip By A Movine PLATE 


For a second illustration of the hydrodynamical equations let 
us take the case of a liquid contained between two flat plates 
of infinite extension, perpendicular to the Y-axis. Let the 
lower plate (y=0) be at rest, while the upper plate (y = A) moves 
uniformly with the velocity a along the X-axis. Since the liquid 
is dragged by the upper plate, so that its velocity increases along 
the Y-axis, we will write for the velocity components 


u=C0,+Cy, v=0, w=0. 
The friction upon any plane parallel to the plates is, per unit 
area, 


d 
Kay =uQ,. 


VOL. I G 
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For y=0 we have u=C,, so that the relative velocity of the 
liquid and the lower plate is C,, and therefore the force exerted 
by the latter upon the liquid AC,. Thus the first boundary 
condition becomes 


du 
p dy = AC}. 
This gives 
Cz = C. 


For y=A, u=C,+C,A, so that the relative velocity of the liquid 
and the upper plate is a—(C,+C,A), and the force exerted by 
this plate upon the liquid \(a—-C,—C,A). The friction of the 
liquid against this plate is —„C,, and the second boundary 
condition thus becomes 

pC, = Ala = Ci = CA) 
or, substituting the value of C,, 


Aa 
Ga 2u +AA 

whence the friction i 
pO ,=5" — 
+ +A 


For )} =o we have (C,=a/A, so that for finite à the friction will 
be somewhat smaller. 

Notice that u/à has the dimensions of a length. (This follows 
from the equation pdu/dy=XC;,, since C, is a velocity, to wit, 
C,=u for y=0.) The physical meaning of this length can be 
seen by imagining that either plate is moved away from the 
other over the distance w/\, while the liquid expands so as still to 
extend from plate to plate and its state of motion remains un- 
changed. The velocity of the liquid relatively to either plate 
is then nil at both boundaries. In fact, u =0 for y= - w/X, and 
u=a for y„=-A+u/XA. Thus the solution of the problem with 
sliding can be reduced to that of the problem without sliding, 
provided the liquid is given the said expansion. 

Also the result of the problem of a liquid flowing through a 
tube of circular section is in harmony with this property. For, 
if R in the expression (7) is replaced by R+ /d, the expression 
(8) follows, provided u/A is small in comparison with R so that 
the square of u/AR can be neglected. 
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4. EFFECT or SLIDING IN THE CASE OF TRANSLATIONAL 
MOTION OF A SPHERE IN A LIQUID 


We will now determine the resistance offered by a liquid to 
a sphere endowed with uniform rectilinear motion. This will 
lead us to the famous formula of Stokes which, among other 
things, comes into play in the lately developed theory of Brownian 
movement. It is still an open question how far the validity of 
this formula can be upheld for very small particles and irregular 
motions. This will still be discussed in the sequel (Art. 7). 

In dealing with the problem in hand we will assume that not 
the liquid as a whole is at rest and a sphere moves through it 
but, inverting the relations, we will imagine that the sphere is 
at rest and the liquid moves past it, having at infinity a uniform 
rectilinear motion. We put the Z-axis along the direction of 
this motion and take the centre of the sphere as the origin of 
our co-ordinate system. Thus, if u, v, w be the velocity com- 
ponents of the liquid, we have at infinity u =0, v =0, w =a. 

From our previous equations (6) 

pAu = 2, uAv = ne pAw =o, 
and from the equation of continuity (1) we derive 
Ap=0. 


To solve these equations we will follow Kirchhoff and intro- 
duce an auxiliary function ®, such that 


pA® =p. 
To begin with, we could try as a solution 
o> =o om 
U=——, = =. 


Ox oy” U= j 
This satisfies the equations of motion, but not the equation of 


continuity. In fact, since 


ðu ov dw 
en 


the equation of continuity could only be satisfied if A® =0, 
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t.e. if the pressure were everywhere nil. We try, therefore, to 
help matters by adding new terms and write 
om , a, 


od , 
a A R ee »- «+ (9) 


uy = 
Ox 


If the new terms can be so chosen that 


Aw’ =0, Av'=0,Aw’=0 . . . (10) 
and 
Ou’ Ov’ aw’ p 
a dy a ee a : ; ‘ (11) 


our solution will be ready. Now, noticing that we must have 
Ap=0 and that this is satisfied by spherical harmonics of which 
the simplest are 1/r and all partial derivatives of 1/r, we try to put 


p=) 


Account is here taken of the circumstance that the pressure 
p must be an odd function of z, since it must have different 
signs at the opposite poles of the sphere. Further, the function 
of z chosen is the simplest which satisfies Ap=0. Finally, a 
constant factor is inserted which will presently be given a 
suitable value. An additive constant would for the problem in 
hand be without significance, and is therefore omitted. 

Equations (10) and (11) are now satisfied by 


9 
w =0, v =0, w = -”. 
r 
Next, we put 
0 
@ =az +b, O om. 


The first two terms do not contribute to A®, and since 
Or 1 
A(ox) =e ox (Ar ‘= 25 ) 


the last term gives 
Ab =f; 


Notice that all three terms of ®, though odd functions of z, 
are even functions of x and y, which harmonises with the 
symmetry of the liquid motion. 
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The first term gives at infinity, where the derivatives of the 
remaining two vanish, 
w=, 
oz 
Thus is a introduced into the formula. The middle term enables 
us to satisfy the boundary conditions at the surface of the 
sphere. By means of these conditions the constants b and c can 
be expressed in terms of the velocity a and the radius R. 
Thus the solution becomes 


Hig P ay 


al 
en) Ee 


Thus far we have followed Kirchhoff. In the boundary conditions 
at the surface of the sphere we will take account of the sliding 
and in that deviate from Kirchhoff. Since at the surface of the 
sphere the velocity of the liquid must be tangential, the first 
boundary condition is 





ux +vy +wz =O, for r=R, 
whence 9 2 


a+” p7O. . å ‘ ‘ (13) 

Equations (12) show that there is symmetry around the 
Z-axis. Let us then consider a point P at the surface of the 
sphere in the plane XZ, for which <POZ = @, and let us introduce 
a new orthogonal system of axes, of which OZ’ passes through 
P and OX’ lies in the plane XZ. The equations (12), when 
transformed to the new axes, become 


u=- (a—",-2)sind + (3B -S)ar(e cos 8-2’ sin 6) 


v = _ Sy cos #-a’sin 6)}. (14) 


w = (a -2 -°) cosd + © -Se cos 0 -x' sin 0) 
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The velocity of sliding at P follows from these equations for 
x' =y'=0, z =R. Thus, and by (13), its components are 


w= -(a 4-9) sin 0, v’ =0, w’ =0. 


Whence, the tangential tension at P exerted by the sphere upon 
the liquid, 
N a2 -2) sin 0 
R3 R i 


and the friction component, which must be in equilibrium with 


this force, w Ow" 


X mulaz tir) 


Thus the boundary condition becomes 
Ou’ Ow" b e\. 
a = -A0 -77 p) n 8, 


where z’ =R and 2’ =0. 
Evaluating the left-hand member by means of (14), one finds 


„a- r sin 0, 
and the second boundary condition assumes the form 

6ub b c 

He =Manp-h) (1B) 
Since @ has disappeared, this boundary condition can at once be 


satisfied all over the surface of the sphere. 
Formulae (13) and (15) give 


2u 
b__0_,0_,,."R 
p Buy R *% gu 
41475) 1+yp 


Thus b and c are expressed in terms of known data, and sub- 
stituting these values into (12), we find also u, v, w. 

The total force exerted on the sphere by the liquid will fall, 
by reasons of symmetry, into the Z-axis and can thus be found 
by integrating Z, over the surface of the sphere. Now, Z, can 
be determined, in two ways, either from 


Zy =Z, sin 0 +Z, cos 6, 
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or from the tangential force X, and the normal force Z, 
leading to 


Z,=Z, cos 0 - X, sin 0. 


To choose the latter way, Zy must be evaluated from the 
equation 
, Ow’ ; ; 
z= -p +2p o,s for z’ =0, 2’ =R, 
which gives 
+ Guc 12ub 
Zr = ga 08 0-24 cos 0, 


and since we have already found 


oe 1 ain, 
we have 
6ub . 6uc 
Zr = T (1-3 cos? 0) +a cos? 6, 


whence, by integration over the surface of the sphere, the required 
resistance, 


W =27R? "zZ sin 6d0 =8rcp, 
0 


1.6. 9 


For A =œ this expression reduces to 
W =6rpaR, 


which is Stokes’ well-known formula. 
Such, then, is the resistance experienced by a sphere moving 
with velocity a through a stagnant liquid. 


5. Liqgurp Motion DUE TO AN IMMERSED VIBRATING PLATE 


We shall next investigate how far this resistance formula 
can be assumed to hold for a sphere in non-uniform motion, 
such as the Brownian movement. The problem offers con- 
siderable difficulties, and we shall, therefore, confine ourselves 
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to the simple case in which the body has a vibratory motion of 
translation. 

As an introduction we will consider a flat plate in the 
YZ-plane maintained in non-damped vibration in the direction 
of the: Z-axis. Let the plate be unlimited, so that the state 
will be the same all along the Y- and the Z-axes and thus 
depend on x and ¢ only. The velocities being again assumed to 
be infinitesimal, and all relevant magnitudes being functions of 
x and ¢ only, the equation of motion 


Ge ðw dw e _ Op A 
reduces to 
ow Ow 
p dt = p a Pe è i è 5 (1 6) 


As ıts solution we take 
w = geinttBßz, 


Substituting this into (16), we have 
inp =b ‘ať a‘ (1 


whence o 
B=- J PP +i). 
bul ) 
Thus, 


int—4/Np 
w =ae V sel 0, 


of which the real part is 


„Ve, ip 
w=ae — ™ cos (ni v 2); 
2u 


representing waves which issue from the plate along the positive 
X-axis, this result being obtained by taking for 8 as solution of (17) 
the negative root. The wave-length is 2rV/2u/np. That this 
has the dimensions of a length is manifest by (16). The amplitude 
of the oscillations of the liquid is thus decreasing considerably 
when these are propagated over a wave-length, namely, to e`?" 
times its original value. 
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6. Errect oF FREQUENCY IN THE CASE OF A SPHERE 
VIBRATING IN A LIQUID 


Passing to the case of a sphere vibrating in the direction of 
the Z-axis, we may represent its velocity component w by ae", 
where 7 is real, corresponding to non-damped vibrations. The 
motion of the liquid must satisfy the equations of motion 


ou _ Op 


+uAu, etc., 


Par Ox 
and the equation of continuity 


DL JE 

Oz oy z ’? 
while at the surface of the sphere, if sliding be disregarded, 
u=v=0, w=ae”, 

Of these equations two particular solutions can be found, 
whose superposition gives a solution which satisfies the boundary 
conditions at the surface of the sphere. 

For the first solution we put p=0. This will then exhibit 
some similarity with the problem of the vibrating plate, since in 
the case of the latter the equation of motion (16) was free of terms 
containing p. 

Let us, therefore, introduce an auxiliary function ® satisfying 
the equation 
om 


Pog = BAY 


and depending on ¢ and r only. The equations of motion will be 
satisfied not only by ® itself but by its derivatives as well, and 
we can put 
u 8b _ PË we (a +0) 
Tae OT iye T  \ dat * ay?" 
Of these the last is so chosen as to satisfy the equation of 
continuity. 
The differential equation for ® becomes, in polar co-ordinates, 
rd) (rð) 
Pron oR” 


this being of exactly the same form as the equation of motion 
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(16) of the liquid containing a vibrating plate. The solution is 


thus b 
Ö =- eint+ Br 
T 


where B= — Se +?) and b is a constant which for the present 


may be left undetermined. Thus far the first solution. 

The peculiarity of the second solution is that it annihilates 
the last terms in the right-hand members of the equations of 
motion. As such we introduce a function y of ¢ and r which 
satisfies the equation 

Ax =0. 


If then u, v, w are equalled to y or to any of its derivatives, 
we shall have Au = Av = Aw =0. 

In order to satisfy the equations of motion and the equation 
of continuity as well we put 


By Oy 


pu=2X, y= 2X p=- 
0202’ gyoz POT ag? 0208 


and we take for x the simplest harmonic function ( 3 multiplied 
r 
by the factor e'*, to express non-damped vibrations, and by the 
density p, to do justice to the equations of motion ; in fine, 
CD : 
x= = en 


where the constant c is again left undetermined for the present. 
Thus the general solution of our equations becomes 
oD +1 1 0% 
~ dxdz * p 0x02" 
OD. 1 0% 
“iya p 5 yae 
_ (PD  o@) 1 0% 
E - (Da 7 sya) ea 022 
The amplitudes of both functions can be compounded by putting 


f= mberat, 
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Next, if we introduce the amplitudes w’, v’, w’, defined by 
u=w'e™, etc., these will satisfy the equations 


Of ya Of w Of Of 
=a Tayar ” =- (aitai a): 


The boundary conditions at the surface of the sphere are, in 
absence of sliding, u’ =0, v' =0, w’ =a, for r=R. 
The first of these gives 


whence 


(3 -3BR + fB2R2)bF8 +38c=0. . . . (18) 


The second, v’=0, gives the same equation, while the third 
boundary condition leads to 
2 df wry df _ z? +y? df 
rdr r® dr r dre” 


or, since df/dr —rd?f/dr? =0, 


-" Wa, for r=R, 
whence 
(1 — BR)be®* +¢ =10R’, . . (9) 


From (18) and (19) we have, for 6 and c, 


RR. (20) 
o= ($ -4AP + 4ER «©... (21) 


Thus the state of motion is determined. 

In view of the symmetry with respect to the Z-axis the 
resistance opposed by the liquid to the motion of the sphere can 
be represented by 


ƏrR?| Z, sin 640, 
0 


where # is the angle between the normal n (of the sphere) at a 
point of the XZ-plane and the Z-axis 
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If we put Z,=Z,e™, Z,=Z,e™, etc., then 


= TA + zZ, )s 


Z, ul + =) 


„=-p + 2u- 


we 


Evaluating Z,, by means of the condition df/dr —rd?f/dr? =-0, for 
r=R, the equations (20), (21), and the relation (17), we find for 
the required resistance, K = 21k? /"Z, sin 6d6, 


K=4rR(-3+3BR-ıPß2Re)uae. . . (22) 


For very slow vibrations, z.e. for very small values of SR, this 
reduces to Stokes’ formula 


K = —G7pRae™ = —6rpRwephrere. 
For any frequency of vibrations (22) can be written 
K=(g +ih)uael", 
The velocity of the sphere will then be 
w= COs nt, 
and the force opposed to the sphere by the liquid, 
K =pga cos nt -pha sin nt. . . . (23) 


Only the first term is to be considered as a resistance, if “ resist- 
ance ” be so characterised that the work done by it over a 
complete period is negative. Now, the work of K is 


f Kwdt = f (uga? cos? nt - uha? sin nt cos nt)dt, 


and here the first term, whose coefficient g will presently appear 
to be negative, gives on integration over a complete period a 
negative quantity, whereas the second term gives nothing. The 
term — pha sin nt is a force proportional to the acceleration ; the 
‘corresponding effect is thus an apparent increase of the mass of 
the sphere, which is due to the co-vibration of the liquid. 
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Let us still consider the resistance uga cos nt. Here g is the 
real part of 
4rR( -3 +38R -48*R?), 


g= -(1 je 


g= -6rR(1+R w), 


and since 


we have 


and the resistance becomes, with w, written for the velocity 
of the sphere, 

—6rRyw(1+R /™P) 2. 0.0.04 

n Ruu ( + J 7 (24) 


This is again the resistance according to Stokes, but increased 
by a term which for high frequencies can outweigh the Stokes 
resistance. 


7. THE QUESTION OF VALIDITY OF STOKES LAW FOR 
BROWNIAN MOVEMENT 


We can now ask what this resistance is like when the sphere 
is endowed with any variable motion. This is of importance 
in connection with Brownian movement. But we must limit 
ourselves to small velocities, so as to be able to neglect, as above, 
all terms such as udu/az, etc. If the velocity were known as a 
function of the time for the whole duration of the experiment, 
it could be developed into a Fourier series, and the result (24) 
could be applied to each term. Thus the resistance would 
be found. This, however, is not very helpful for a general 
discussion. 

Stokes’ formula cannot be applied to the case of Brownian 
movement, as this is much too quickly variable for such a 


purpose. The term R J A in (24) is to be neglected in presence 


p 
of unity, if the vibration time T is large compared with wpR?/y. 


If we write mpR?/u=6, then it will be possible to apply Stokes’ 
law for slow vibrations, for which, that is, the vibration time 7 
is large in comparison with 6, and this holds also for other 
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motions, provided the velocity does not change much during the 
time 6. 

To illustrate this by an example, let us see whether Stokes’ 
law is applicable to the extinction of the motion of a sphere due 
to the friction of the liquid. If m be the mass and v the velocity 
of the sphere, Stokes’ law would give 


my = —ôrukv, 


whence R 


6nuR, 
= ve” m ” 
so that v would dwindle down to the e-th part of v, after the time 
m/6arpR. If p, be the density of the sphere, this time is 
m _2p,R? 


T= onu 9 p 





If p and p, are comparable with each other, @ is seen to be of 
the same order as the time r. Stokes’ law cannot, therefore, be 
applied to the extinction of motion here considered. 

For particles in Brownian movement 0 becomes quite small. 
If R=5.10-° and y =18. 10-5, then 


-10 
ga” 25:10", 4 10-4 


approximately. Yet the motion of the suspended particles will 
vary during this time considerably, so that the refinements of 
the Brownian movement cannot be mastered by the law of 
Stokes. In many other cases, however, the motion within the 
time @ will change but little and the law will be applicable. 


8. DEDUCTION oF EINSTEIN’s FORMULA FOR THE MEAN SQUARE 
OF THE DEVIATION OF A PARTICLE IN BROWNIAN MOVEMENT 


For the investigation of Brownian movement it is of import- 
ance to correlate the mean square £? of the distance attained by 
a particle within a given time ¢ with the properties of the liquid, 
1.e. with the viscosity coefficient », and with the radius È of the 
particle. 

Without following the particle’s actual crinkly path and its 
rapidly variable motion, we can find the required connection by 
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a roundabout way, namely, by considering the diffusion velocity 
of the particles. On the one hand this can be expressed by u, 
and on the other hand by £. Evaluating thus the diffusion 
coefficient by two different methods and equating the two 
expressions, we shall find Einstein’s formula * for £2. 

First method. Suppose we had a liquid containing suspended 
particles whose concentration varies in some direction or other. 
Along that direction put the X-axis. The concentration of the 
particles is given by their number n per unit volume. Thus n is 
a function of x. The particles will exert an osmotic pressure. 
Now, if it be assumed that the mean kinetic energy of a particle 
is equal to that of a gas molecule at the same temperature T, 
v.€. 3kT, the osmotic pressure is 


p=nkT, 


this being % of the total kinetic energy of the particles per unit 
volume. 

The force driving the particles is equal to the difference of 
their osmotic pressure in two planes perpendicular to the z-axis. 
In a stationary state of diffusion this difference will be balanced 
by the force exerted by the liquid upon the particles. 

Apart from the Brownian movement the particles have a 
common velocity along the X-axis if the concentration decreases 
with increasing x. 

Since the formulae for the stress components and the equations 
of motion of the liquid are linear (products of velocities and 
derivatives of velocities being omitted, as they are small), we 
may say that the force exerted by the liquid upon the particles 
consists of two parts, one corresponding to their common velocity, 
and another to their Brownian movement. The latter will vanish 
for all the particles taken together. 

Thus we have dp 


dx 


= brp Rn, 


* Ann. der Physik, vol. xix., 1906, p. 371. 

t The osmotic pressure is, properly speaking, the momentum along the 
X-axis transferred by the moving particles in the positive less than in the 
negative direction, per unit area perpendicular to this axis and per unit time. 
Consequently, the momentum of particles contained in a layer of thickness dz 
undergoes the change — Pax per unit time, t.e. as if a force of this magnitude 


acted upon the particles. 
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where the left-hand member represents the gradient of the 
osmotic pressure, and the right-hand member the resistance 
opposed by the liquid to n particles ; and since p =nkT, 


dn 
_ kT di = Sru Rn. 
Whence the diffusion current 
ree 


where x, the coefficient of diffusion, has the value 


_ kT 
"6muR 
The diffusion is thus calculated by means of the law of Stokes. 

In the second method we concentrate all our attention upon 
the Brownian movement. We do not consider, however, the 
actual path of a particle, but its total displacement within a 
certain time. Let the distance attained after a time t, reckoned 
from some initial moment, be s. Then s is different for different 
particles. By means of probability considerations it can be 
shown that the mean value of s? is, for all particles, proportional 
to L, say, 


(25) 


K 


s? = Bt. 


The coefficient 8 can be observed, as, among others, was done 
by Perrin.* 

Thereupon will our diffusion theory be based. Suppose that 
the concentration varies from point to point ; then the particles 
from a small volume-element will spread after a time ¢ over a 


sphere of radius vV 8t. Thus the concentration differences will be 
gradually obliterated. 
Now, it can be shown that 


1 
mages e+ e 


where £ is the projection of s upon the X-axis.} 


* Comptes rendus, Paris, vol. cxlix., 1909. p. 477, and vol. clii., 1911, p. 1569. 
t Suppose that during the time t all particles are displaced along the X -axis 
over the same distance l, so that for one half of the particles x is increased, 
and for the remaining half diminished by J. Next, consider a plane V, per- 
pendicular to the X-axis, and two layers, each of thickness /, on both sides of V. 
Let N, be the number of particles contained, at the beginning of the time ¢ 
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From (25) and (26) follows Einstein’s formula 


@ kT RT 1 
t 3nmuR N 3nuR Bo Sali ae 
where È is the gas constant and N the number of molecules per 
gram molecule. 

By means of this formula Perrin determined N from his 
experiments. 


and per unit area of V, in the layer on the positive side, and similarly N, 
on the negative side of V. Then the diffusion per unit time is 4(N, — N,)/t. 
We can assume that within these thin layers 7 is a linear function of x, and shall 


2 
thus find, for the diffusion, “> i and, for tho cooflicient of diffusion, 
lp? 


ar 
which obviously must be replaced by (26), when account is taken of the 


diversity of the č-values for different particles. 
It will be readily scon that f= 48°. 


CHAPTER II 
FRICTION AND SLIDING, TREATED KINETICALLY 


9, FRICTION INDEPENDENT OF THE DENSITY OF THE GAS 


RETURNING to the questions concerning inner friction and sliding, 
we will now treat them on the kinetic theory. Moreover, instead 
of a liquid, we shall now consider a gas. ' If this be strongly 
condensed, there is no sliding. Thus we can at first exclude the 
sliding, to introduce it later on as a correction. 

We will begin by proving that the inner friction is independent 
of the gas density. For this purpose we consider a simple case 
of motion, viz. a gas whose horizontal layers are, as a whole, 
shifted over each other. Let us introduce a co-ordinate system 
whose XZ-plane is parallel to these layers, the X-axis pointing 
in the direction of streaming. If wu be the velocity of a layer, 
the motion of the gas can be expressed by the equation 


u=cy, 
where c is a constant, it being assumed that the layer y=0 is 
at rest. 

The state of a molecule will be determined by £, n, &, the 
velocity components, and x, y, z, the co-ordinates of its centre 
of gravity. The relative co-ordinates and velocities of the parts 
of a molecule with respect to its centre of gravity can here be 
disregarded, as they do not affect our problem. Let us consider 
the molecules which at the instant ¢ are contained in the volume- 
element dS =dadydz at the point P, of co-ordinates zx, y, z, and 
whose velocity components are contained between £ and £+dé, 
n and n+dn, ¢ and &+d£. The number of these molecules can 


be expressed by F(E, 7, ©, y, t)dSda, 
98 
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where dA=d£dnd&£ and F is independent of x and z, since, by 
assumption, the state is the same throughout a gas layer parallel 
to the XZ-plane. 

Further, let bdSd\dt be the number of molecules which are 
thrown into this group by collisions during the time dé, and 
adSdidt the number of those which for the same reason and 
during the same time leave this group. Let external forces be 
absent. At the time ¢+d¢ the molecules of this group will come 
to lie in a volume-element dS’ =dS [by Liouville’s theorem], at 
the point x + dt, y + ndt, z + dt. Thus, 


Fl, n, ©, y, ))dSdA + (b ~a)dSdadt 


will be the number of molecules contained at the instant ¢+d¢ 
in the element dS, constructed at the point x+ &dt, y+ ndt, 
z + dt, and having their velocity components within the domain 
dÀ. Whence, 5 quik OF 

u 
This equation would enable us in general to find the properties 
of the function F. This, however, can be accomplished only in 
the simplest case. If we put F=f for y=0, then 


F(E, N» L, Y, y=0 =f(é, N, L, t). 


If the whole gas is given a translation velocity —cy along the 
z-axis, the state of motion In a plane at the height y becomes 
what it was originally in the plane y=0. Whence we see that, 


generally, FE a, £, y, =f- ey m b t). 

We assume further the state to be stationary, so that 
OF OF 
Ot =(, b -a= ay 


It remains only to determine the state for y=0. From the 
preceding equation we have, for y =0, 


b-a= -0 An, ee (28) 


We will now compare the state of the gas in the said motion with 
the state in which there is no streaming, 1.e. forc=0. For this 
case f can be written Nf, (where N is the number of molecules 
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per unit volume). Then, in the state of motion, f=Nf,+/’, 
where the function f’ contains all the refinement; it may 
assume positive values for some groups, and negative for 
others. The exact determination of f’ is laborious and for 
our purposes unnecessary, since it can be shown by a simple 
reasoning that the friction is independent of the density. We 
assume that the velocity gradient c is small, in other words, 
that there is only a small departure from the state of rest, so 
that f’ is small compared with Nf). Then f in the right-hand 
member of (28) can be replaced by N/, which amounts to 
neglecting only terms of a higher order. This gives 


-ey Io 
b-a=-cN DE Ne 

For the state of rest Maxwell’s velocity distribution * holds, and 
b-a=0. Owing to the change brought about in the function 
f by the motion b-a does no longer vanish. We consider again 
a certain group of molecules and we distribute also the remaining 
molecules into groups according to their velocities and co- 
ordinates. Let the numbers of molecules in these groups be n}, na, 
etc. Owing to the collision of a molecule of the first group with 
a molecule of one of the remaining groups that molecule leaves 
its group. Such a collision then belongs to the type a. Collisions 
of molecules of the remaining groups with each other can increase 
the numerosity of the first group and will thus be of type b. 
Since the number of collisions between two groups of molecules 
is proportional to the numerosity of the first as well as to that 
of the second, we can write 


b ~Q = La NaNe, 


where a; is positive for such collisions as contribute something 
to b, and negative for such as contribute to a. 

For c=0 the numbers of molecules in these groups will be 
denoted by N,9 ao, etc., so that b-a= Fast, N = 0. 

In general we can put 


b -a = Bay +N7')(Nzo +03‘), 


* Viz.: hs — h(£2+4-92+- 
folk» m nel ne nn ”, 


where A is inversely proportional to the absolute temperature. 
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where the a’s remain as before, since we are considering groups 
which are completely determined by their velocities, so that 
when we pass from the state of rest of the gas as a whole to its 
streaming only the numerosity of each of the groups is changed. 

Since we consider only a small departure from the state of 
rest, n,’, ng, etc., are small compared with nio na, etc., and 


therefore 
b - a = Zaya(Niong + Nagri’). 


This leads to the equation 


A, 


Eainig +Nagn, ) = -eN oy, . . (29) 


by means of which the variation of f or the numbers n,', ng, 
etc., can be determined. If the molecules are divided into K 
groups, there are K equations such as (29). 

If we change c, the equation (29) remains satisfied, provided 
n,’, Ng, etc., are changed in the same ratio; thus the whole 
disturbance, and therefore also f’, will be proportional to the 
velocity gradient. Now, if the density be varied, and thus also 
N, the number of molecules in each group for the state of rest, 
Nip Np etc., varies in the same ratio. Then (29) is satisfied, 
provided n,’, n,’, etc., remain the same, so that also f’ does not 
change. 

Whence we see that f’ is proportional to c and independent 
of N, that is to say, independent of the density. 

The friction upon a plane parallel to the XZ-plane is, per 


unit area, 
m fén fdà. 


In fact, this is the momentum along the X-axis, transferred per 
unit time across a unit area parallel to the XZ-plane. Since in 
the state of rest there is no friction, 


mJ én fodd =0, 
and the last expression for the friction can also be written 
mf Enf’da. 
This, as f’ itself, is proportional to c and independent of N, and 


can, therefore, be written ne 


where u is the coefficient of friction. The latter is thus seen to 
be independent of density. 
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10. UNIFORMITY CONSIDERATIONS 


That friction is independent of density can also be proved 
by another method which, though not going so deep into the 
phenomena, is yet exact. We shall in this case also confine our- 
selves to a simple example, but taking at the same time account 
of the sliding. We consider two infinitely extended horizontal 
plates, of which the lower is at rest and the upper moves with a 
constant velocity v horizontally, towards the right, and we pro- 
pose to find the motion of the gas contained between the plates. 

Fig. 1 gives a graphical representation of the motion of the 
gas for the case in which there is no sliding along the walls. The 
y abscissae are the velocities in the different layers, and the 
ordinates the distances from the lower plate. The velocity 
v of the upper plate is represented 
by AB. Manifestly, the line join- 
ing the end-points of the velocity 
vectors is a straight line. 

A sliding along the walls will 
have the effect that the velocity 

O Fia. 1. X of the gas layer in contact with 
the upper plate will be smaller. 
If this be represented by AF (Fig. 2), then AF<AB. Similarly, 
the velocity of the layer touching the lower plate will no longer 
be nil, but will have some small value OE. The distribution 
of velocities over the different layers will no longer be represented 
by a straight line but by the line ZF in Fig. 2, which, when the 
plates are far enough apart, has a straight portion, but is curved 
at the extremities. 

To enter somewhat deeper into the sliding along the walls, 
we must keep in mind that by the velocity of a gas layer is 
meant the average velocity of the particles of which the layer at 
a given instant consists. Among the particles of the lowermost 
layer there will be some which previously belonged to a layer 
with some velocity of streaming and were carried down to that 
‘layer by collisions. For these a ‘certain direction of motion, 
viz. that of the streaming velocity, will be privileged. In this 
layer, however, there will also be present, at the given instant, 
some particles which were there before and which collided with 


B` 
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the solid wall. We have, therefore, to distinguish also different 
cases of collision with the wall. If the wall is a perfect reflector, 
then the particles retain after the collision their tangential velo- 
city component. But in general even the best polished walls will 
have to be considered as rough with regard to impacts of mole- 
cules, so that the reflection will be of a diffuse nature, and there 
will be no privileged direction of motion for the rebounding 
particles. But, as we just saw, there is such a direction for the 
impinging particles, so that, all things being considered, the 
layer in contact with the wall will have a velocity in the direction 
of streaming, very much as in the case in which there is sliding. 
The precise form of the curve representing the velocity distribu- 
tion is hard to deduce from these considerations, at least for 
points near the plate, where the state of motion is very com- 
plicated. Among other things one would have to take into 
account that some of the molecules colliding with the wall 
adhere to it. But the farther away from the wall, the less the 
irregularities, and if we assume the mutual distance of the plates 
large compared with the mean 
free path of the molecules, a con- 
siderable portion of the curve will 
be straight. If the distance were 
very small, the straight-line portion 
would entirely disappear. 

In Fig. 2 the velocity curve is 
represented by HF. (Probably 
its shape should be somewhat 
different, such, e.g., as indicated by 
the dotted line.) If the straight- 
line portion of the curve be pro- 
duced up to the points C and D, 
the former being on the Y-axis and the latter having the 
abscissa AB=v, then C and D will have the following signifi- 
cance. If the gas were replaced by an imaginary gas layer 
extending also beyond the plates and having throughout the 
same constant velocity gradient as that prevailing in the actual 
case only between the plates, then the layers of gas at C and 
D would have the same velocities as the plates, t.e. zero at C 
and vat D. (Cf. Art. 3, p. 82.) 

We will now prove that the segments OC and BD are 





Fia, 2. 


| 
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independent of the distance of the plates. For this purpose let us 
imagine a plate P inserted at the height of G and moved in such 
a way that the state between G and A is not changed, except, 
of course, in the immediate neighbourhood of the plate P itself. 
We make the arc of the velocity curve near P congruent to that 
at E. Now if HK =OE, the velocity distribution represented by 
KF can prevail, provided GH is equal to the velocity of the plate P. 
In fact, if the whole system is given a common velocity equal 
and opposite to GH, the velocity curve retains the same shape, 
but is simply shifted over the distance GH to the left. The state 
at the plate P will then be exactly the same as that actually 
prevailing at the lowermost wall. Whence it follows that OC 
and BD are independent of the mutual distance of the plates. 

In much the same way it can be shown that OC=BD. In 
fact, if the whole system is given a velocity equal and opposite 
to AB, the upper plate is brought to rest, while the lower one 
will have a velocity v towards the left, and the velocity curve 
will be shifted over a distance AB in the same direction. Thus, 
if the figure is turned around by 180°, it must be exactly the same 
as that representing the actual state. The shape of the curve at 
E must, therefore, coincide with that at F, and thus OE = FB 
and also OC = BD. 

Further, if the velocities of the plates are supposed to be very 
small compared with those of the molecules, the state of the 
gas can be considered as an infinitesimal deviation from the state 
of apparent rest. This deviation can then be put proportional 
to the infinitesimal cause, t.e. to the velocity v of the upper plate. 
All horizontal lines of the figure can thus be magnified or reduced 
in the same ratio as v without changing OC and BD. Whence 
we see that OC and BD are independent of the velocity v. 
Again, as we saw before, OC and BD remain also unaffected by 
the change of the distance between the plates. They can thus 
depend only on the nature and the density of the gas. All this, 
of course, holds only if the velocity curve has a rectilinear portion, 
and ceases, therefore, to be true when the distance between the 
plates is of the order of the mean free path of the molecules. If 
we put OC = BD =r, then v will be characteristic for the gas. 
The meaning of v is, that if the velocity gradient were constant 
not only between but also at the walls, the gas at a distance v 
from the plate would attain the velocity of the plate itself. 
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Let the velocities of the plates be zero and v, and their mutual 
distance A. Then, by what precedes, the velocity gradient of 
the gas (i.e. the change of velocity per unit length in the direction 


of the Y-axis) will be en , and, therefore, the friction per 
A+2p 


unit area, parallel to the plates, gl 
A+2p 

We will call v the coefficient of sliding and we will prove that 
it is inversely proportional to the density p, while it will be 
shown once more that u is independent of p. ’ 

For this purpose let us consider a second system which will 
be denoted by II., while the original system will be denoted by 
I. ; the magnitudes relating to the system II. will be distinguished 
by dashes. Let the system II. be so chosen that the correspond- 
ing vertical distances are k times smaller than those in I., while 
the velocities at corresponding points are the same. We will 
now prove that the system II. represents a possible state of 
motion, provided the density in II. is & times that in I. 

Let the function determining the state at a point P of the 
first gas be F,(£, 7, ¢), and that at a corresponding point P’ of 
the second gas F,(£, », £). Then the condition for the density 


ill imply th 
wilimpiythat Fé, 1 O=kF AE, 7, 0). 


Since this holds for each group of molecules, the velocity of 
streaming, t.e. the mean value of £, will, as we assumed, be the 
same at P’ as at P. 

We have now to ascertain whether the condition of a possible 
state of motion is satisfied, and for this purpose we consider the 
equation for a stationary state of the gas, viz. 

b-a nr je pe 
~ a> * dy Oz °° 

If the X-axis is again drawn towards the right, and the Y-axis 
perpendicularly to the plates, this equation becomes 

oF 
b -a= dy T 

It holds for the state I. and, as is easily seen, also for the 
state II. In fact, for the state II. both sides of the equation 
become k? times greater. For the left-hand member expresses 
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the number of collisions (and here the gas may be assumed 
homogeneous around P’), and if the number of molecules in each 
group is k times as large, the number of collisions is increased 
k? times. And with regard to the right-hand member, notice 
that 7 is the same for both systems, and that OF,/oy=k20F Joy, 
since F,=kF, and since all dimensions along the Y-axis in the 
system II. are k times smaller than those in the system I. Thus 
the equation of the stationary state within the gas is satisfied. 

It remains, however, to consider the state at the walls. The 
number of molecules of a given group which, per unit time, strike 
against a wall, will for the system I. be determined by 7F,dA, 
where F, indicates the state at the considered point of the 
boundary layer. Similarly, this number for the system II. will 
be given by 7F,dd. Suppose, further, that the two walls are 
perfectly equal; then the collisions of molecules of the group 
considered will, in the system lI., be & times more numerous 
than, but otherwise entirely the same as, those in the system I. 
Also for the molecules rebounding from the wall will the numbers, 
for groups with the same values of &, n, £. in the system 1I. be k 
times those in the system I. This, moreover, will be the case for 
perfectly reflecting as well as for diffusedly reflecting walls, the 
latter in contrast with the former being such that the molecules 
after reflection retain no trace of the mean motion relative to 
the wall which they had before the collision. If the above is 
assumed to hold also for the case that there is an adhering gas 
layer at the walls, then the system II. will represent also at the 
walls a possible state of motion. 

Since the vertical dimensions in IT. are $ times smaller than 


those in the system I., we have (Fig. 2) B’D' = „BD or 
Ze 
v= Ẹ 
Thus it is proved that the sliding is inversely proportional to 
the density. 


The friction, per unit area, upon a plane parallel to the 
XZ-plane will, as we saw before (p. 101), be represented by 


w =/ménF dd, 
whence w = kw. 
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Since the velocity gradient in the system II. is k times that in I., 
and since w = u times this gradient, we have 


W =h, 
which proves that the friction is independent of the density. 


ll. KUNDT AND WARBURG’S EXPERIMENTAL INVESTIGATIONS 


In 1875 an experimental investigation on friction and heat 
conduction in rarefied gases was published by Kundt and 
Warburg,* confirming the preceding theoretical results. 

A round horizontal disc S, in bifilar suspension, was made 
to oscillate in a gas between two fixed plates S, and S,. The 
coefficient of viscosity was determined by measuring the loga- 
rithmic decrement of the oscillations. This method was already 
applied by Coulomb in the case of a disc oscillating in an un- 
limited gas mass, while Maxwell improved it by introducing the 
two fixed plates, which has enabled him to make the calculation 
somewhat more accurate. 

Kundt and Warburg’s formula is 


MDB = Ap(1 +0), 


where JM is the moment of inertia of the vibrating disc, D the dis- 
tance between the disc and one of the fixed plates, 8 the damping 
(appearing in the angular deviation w=ae~™ cos nt), and u the 
viscosity coefficient, while @ is a number which depends on the 
density of the gas and the distance of the plates and which was 
accurately calculated by Maxwell. The spinning disc generates 
transversal friction waves. If the layer is thin enough, the air 
current is everywhere in phase with the disc, and @=0. Such was 
actually the cage in Kundt and Warburg’s experiments. Again, 
A would be equal to rkt for each plate surface (of radius R) 
exposed to the friction of the gas, provided that the friction 
experienced by a surface-element of the disc having a velocity v 
is taken as equal to the friction upon a surface-element of 
an infinitely extended plate moving with the translation velocity 
v between two fixed unlimited plates placed at the same distance. 
The value of A was more accurately calculated by Maxwell. 


* Ann. Phys. und Chemie, clv., 1875, pp. 337 and 526. 
T Phil. Trans., London, elvi., 1866, p. 249. 
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In these experiments the pressure p of the gas and the 
logarithmic decrement à of the oscillations were measured. The 
latter decreases at smaller pressures, this being due to the 
influence of sliding. In fact, if for a large gas density, for which 
we can put v =0, the logarithmic decrement is Ap, we have for A, 
which is E to the friction, for smaller densities, 


À -1 a\-1 
Ào -45 -(1 +>) =(1 +5) i 
since a_ 2v 
p D 
The factor a will, for a given distance of the plates, be a constant. 
Its value will vary inversely with this distance. 


Here are the results of some of the experiments in which the 
gas was air, and D =0:1104 cm. : 


p A 
320 mm. 0:1318 
20 0:1306 
TT 0:1292 
76 0:1292 
2-4 0:1256 
0:63 0-1109 


For a large number of observations on various gases and at 
different pressures, Kundt and Warburg actually succeeded in 
choosing a so that the values of A calculated by means of a, p, 
and A, agreed with the observed ones. Some of these results 
were : 


p X obs, A calc. 
20 mm. 0-131 0-131 
7-6 0:129 0:129 
2-4 0:125 0-124 
1-53 0-120 0-120 


For different values of D the value of a appeared to be with 
fair accuracy inversely proportional to D. Thus for three 
different distances the values of a, when derived from the 
observations, were : 

=()-149, a, =0-070, a, =0-061. 


eS LS ETT un m name SrtA CD AAD Po ih Ua Md 


* Cf. p. 105. 
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On the other hand, starting from a, and assuming that a, and 
u, are equal to a,D,/D, and a,D,/D,, the result was : 


a = 0:084, az = 0-059. 


For still smaller pressures we are reduced to the case in which 
the straight portion of the velocity-curve (Fig. 2) is absent. 
According to Kundt and Warburg, the formula for A continues 
to hold as long as the distance of the plates is at least 14 times 
the mean free path of the molecules. 

Lastly, they find for air at the pressure of 76 cm. and the 
temperature of 15° C. 

v =0-00001 cm., 


which shows that the sliding at normal pressure is small. But 
at the pressure of 1 mm. one should then still find v =0-0076 cm., 
so that for a distance D =0-2 cm. of the plates the effect of sliding 
would already be very marked. Indeed, the term 2» in the 
M? _ could not be neglected in the presence of D. 

D+2p 

The mean free path of the molecules at normal pressure is 
about 0-0000084 cm., and thus of the same order as v. 


expression 


CHAPTER III 
KNUDSEN’S INVESTIGATIONS ON RAREFIED GASES 


12. FLow oF A RAREFIED GAS THROUGH A 
NARROw TUBE 


WHILE in what precedes we have considered the extreme case in 
which the gas density is large, we will now turn to Knudsen’s 
investigations * which represent the other extreme case, corre- 
sponding, that is, to such a small density that the effect of 
collisions can be disregarded. Knudsen investigates, theoretic- 
ally as well as experimentally, the flow of such a strongly 
rarefied gas through a cylindrical tube whose walls are sup- 
posed to be so rough that the rebounding molecules have no 
trace left of their original streaming motion, being reflected from 
the walls in all possible directions. 

Obviously the problem in hand can be treated on similar lines 
to that of heat radiation. Analogously to the latter we will write 


A cos ddwdo 


for the number of particles rebounding from an element do of the 
wall within a cone of aperture dw, the axis of the cone being 
inclined at an angle @ to the surface normal. This amounts to 
assuming that, with respect to the reflected particles, the normal 
of the surface-element of the wall is privileged. The factor A 
depends on the total number of particles rebounding from the 
element do of the wall, and this is, for a stationary state, equal 
to the number of particles striking this wall-element. Thus A 
is a function of the state of the gas. 


* Ann. der Physik, xxviii., 1909, p. 75. 
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Similarly as for heat radiation, the number of particles which 
two elements do and do’ send to each other can be represented by 


A cos 0dwdo=A cos 0’dw'do’ = 3 A cos @ cos O’dodo’, 


where r=PP' (Fig. 3). 

Thus, in a space occupied by a strongly rarefied gas a state of 
stationary equilibrium, analogous to that of black body radiation, 
is possible, provided the walls are kept 
at a constant temperature. For all ele- 
ments of the walls the value of A is the 
same, and the last expression 
holds also for the number of 
particles which pass through 
two surface-elements placed 
within the gas. 

We will assume for the gas de Fia. 3. 
Maxwell’s law of velocity distri- 
bution, and consider, in the first place, the number of particles 
which traverse, per unit time, a surface-element do at a point P 
within the gas, and whose direction of motion falls within a cone 
of angle dw and axis coinciding with the normal of the surface- 
element. This number can be written 


Adodw, 







and, if f(v)dv be the number of particles per unit volume whose 
velocity falls within the interval v to v + dv, 


vdo aa f(v)dv = Adodw. 
0 Aor 
Thus, 


1 ® 
A= izl, vf(v)dv. 
According to Maxwell’s law, 
f) = Ce-i, 


where 4 is inversely proportional to the temperature and C 
depends on the density of the gas. This gives 
C 


A = Sh 
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Since the gas pressure is equal to two-thirds of the kinetic 
energy per unit volume, we have 


-ImC0 f "Adv = ml Ji D 


A=? fh 


curs 


p= 


whence 


MIT T 


Thus A is a function of the temperature [through A] and of the 
pressure. 

At constant pressure, A diminishes with increasing tempera- 
ture. This can be seen directly. In fact, at a higher temperature 
the molecules strike more vigorously ; if, therefore, the pressure 
is still to remain the same, there must be a smaller number of 
collisions. 

Let us now consider a gas flowing in a vertical cylindrical 
tube downward along the axis, as in Fig. 4, where the X-axis 
is along the axis of the tube. 

We will assume the temperature to 
be throughout the same. If such were 
also the case of the pressure, the gas 
would remain in the state of equili- 
brium described above, and A would be 
throughout the same. But in our case 
the pressure increases up the tube, and 
the same will also be true of 4. We 
suppose p to be slowly variable, so 
that it can be represented by a linear 
function of x, which amounts to re- 
taining of the Taylor-series development 
of p the first term only. In accordance 
with this assumption of small pressure 

Fia. 4. variations, the values of A will be the 

same in all points of a cross-section of 

the tube. In fact, every surface-element can be assumed to be 
hit by as many molecules as if the pressure were everywhere 
the same as in the neighbourhood of that element. For the 
number of particles arriving from above is as much increased 
as that of particles moving upwards is diminished, so that the 
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total number crossing the surface-element is not changed. A is 
thus a function of z alone, and this holds for every tube, no 
matter what the shape of its cross-section. 

Let us now consider an element do of the section v =0 placed 
at the point P (Fig. 4), and let us find the number of particles 
crossing this element pcr unit time in a bundle of directions 
contained within the cone of angular aperture dw whose axis 
makes an angle @ with the normal. This cone cuts the wall in 
an element placed at a point Q. Thus, the particles in question 
issue from an element do’ of the wall within the cone aperture 
dw’, the axis of this cone making an angle 6’ with the normal, 
so that their number can be written 


A, cos H’dw’do’. 
By what was said before this can be replaced by 
A, cos ddwdo. 


Here A, is the value of .1 at the point Q. Let Q, be the 
intersection point of the generatrix of the cylinder [passing 
through Q] with the contemplated section, and let PQ,=r. Then 


QQ, =r cot 6. 

Again, 

A,=A —-r cot 0 st 
dc 


where A and d1/dz are the values of these magnitudes at the 
section x =0. 
Thus the required number of molecules will be 


A cos Odwdo ae r cos l cot 6dwde. 


Integrating this expression with respect to dw over a whole 
sphere, we shall find the current through the element do. Now, 
the first part contributes nothing to this integral, while the 


second gives 


dA 
-2 Tz I" cos 8 cot @dw, 


where the integral is to be extended over the upper half of the 
sphere. 
VOL. I I 
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Thus we can write for the current across the whole section 


dA 
da’ 


where (=2/do/r cos 6 cot Odw. 


-G 


Now, if & be the angle contained between r and a fixed direction 


in the section, du =sin 6d6dd, 


and 


6=2[de [ “rp 008 d0 =}r | do J “nl. (80) 


This integral is related to the self-potential of a substance 
uniformly distributed over the section. 
In fact, let p be the distance between 
two elements do and do’ placed at the 
points P and P’ (Fig. 5). Then the 
potential at P will be 


Fra. 5. Í = = |  [abep = f "i 


p 
and if this be integrated over do, we have 


[do f "rdo. 


For a tube of circular section the integral is easily evaluated. 
Let O (Fig. 6) be the centre of the 





circle and OA a fixed direction, so that Q 

ZAPQ,=d. Put OP=I and OA=a. 2 

Then Sn, A 
QQ =2V a7- le sin? ¢. Qo _/N\, 


In (30) r can be replaced by Q,.Q); 
provided we integrate over & from 0 
to 7. Thus we find 

Gen | do}. ya? — [2 sin? 6 dd. Fic. 6. 
Now, taking for do a ring of area 27ldl, we find 


a 7/2 oan oe 
G=4m| zii Vali sind dd. 
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Inverting the order of integrations we have 
1/2 ze nenn 
G=4m| dé | Va -Esimga 
0 0 
_ intas f" Leos i$- T Ta COS $ 4.008 @ sin? Pas 


o sın?d sin? $ sin? ġ 


1/2 ` n/2 
= sta? | f - ay +cos 6 |dd = sta? tan +sin N = A720, 
° j2 cos? ® 


Thus the total current across the circular section becomes 


dA 
— §ar?a8 dz’ 
where 4P jh 
MTA T 
If u? be the mean squared velocity, then 
2 
Qu? 


and since the temperature is assumed to be throughout the same, 
we find for the current 

8 [Im aè dp 

~3N 2 mu de 


Here the current is expressed by the number of molecules. In 
order to express it in mass units, the last expression has to be 
multiplied by the mass m of a molecule. This gives 
_8 /3n a dp 
3V 2 udr 
The volume of the gas flowing across the section per unit time, 


at the given temperature and pressure, will be equal to'the last 
expression divided by the density p, t.e. 


3m aè dp 
-3 /® ia = Se (31) 


This expression can be compared with that for the current which 
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would follow by applying Poiseuille’s law to these small densities.* 
This would be, in terms of the volume, 


_ma® dp 
8u dz’ 


where the friction coefficient w is approximately 0-31lpua,7 if p be 
the density, u the root of the mean squared velocity, and A the 
mean free path of the molecules. This gives for the current 
according to Poiseuille’s law 


na! dp 
De en (32) 


Apart from the numerical coefficients, the currents (31) and (32) 
are to each other as X to a, that is to say, as the mean free path 
of the molecules to the radius of the circular section of the tube. 

Since in our case X is very large compared with a, the 
amount of gas flowing through the section per unit time will be 
much greater than according to Poiseuille’s law. In fact, it 
appears from Knudsen’s experiments that it may well be 50,000 
times as large. 

Knudsen gives the current in ternıs of volume, at a pressure 
of 1 dyne per cm.? and a temperature r° of the gas. If p, be 
the density of the gas in these conditions, the current is expressed 
by 

8 /3m aè dp 
BN 2 pu dr 


The relation p=4pu? gives in the present case 4p,u?=1 or 
u=~V3/p,. Thus, if pọ be the BR at the pressure 1 and the 


temperature 0°, and therefore p, = „(@= 1/273), we have 


l - 
ultimately the current 
‚dp 
E 
where 
=; m adı/1l +ar l+ar, 


2 VPo 


* Cf. p. 81, equation (7). 
fT Cf. Maxwell, Phil. Mag. (4), vol. xix., 1860, p. 31. 
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13. KNuDSEN’s EXPERIMENTAL INVESTIGATION ON THE FLOW 
OF A RAREFIED GAS THROUGH A Narrow TUBE 


In Knudsen’s investigation this formula was tested experi- 
mentally. In his experiments the stream passed through a tube 
joining two vessels filled with a rarefied gas. Each vessel was 
connected with a MacLeod manometer by means of which it was 
possible to measure very small pressures by compressing the gas 
to a known small volume, measuring its pressure and thence 
calculating the original pressure by Boyle’s law. For pressures 
smaller than 5 mm. the pressure in each vessel is in this way 
measured separately, while for greater pressures the pressure 
difference between the two vessels is measured directly. 

Let the volumes of the vessels be v,, v, and the pressures 
within them p; and p, and let p, > Pa. For a cylindrical tube of 
length I 

dp _ Pe- Pi, 
de l 
whence the current 
-C dp -0 P Pa 
dx I 


The volume of gas flowing through a cross-section of the tube 
per second is, when measured at the pressure I, 


dp y dp, 
dt? dt’ 


“ki 
these being the quantities leaving, per second, the first and 
streaming into the second vessel, so that 


d -P d — 
- = -oPı en, Uy nO A 


Whence 
d(pı a = — (DP ~ Pe) ra 5 = 
or 
log (pı - P2) = - T m Di 
and 


Pi — Pz =Const. e` -Iita v: 
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g; 
ince g- 1%% Alog(p -p:) 
Ti +V At i 


the constant C can be determined from the change of the pressure 
difference during a given time. 
Knudsen determined from his observations the quotient 


C’=C/l 
whose theoretical value is 
8 /ma?Vl+ar 
y a (33) 

In his experiments the value of A/a mounted up to 6000. He 
worked with hydrogen, oxygen, and carbon dioxide. 

In order to find the effect of the length of the tube, different 
tubes were used, of which the first was about 6 cm. long and 
about 0-02 cm. in diameter, while the second was about twice as 
long and had the same diameter. The gas in this case was 
hydrogen. 

The ratio of the observed values of C’ thus found was 1-95, 
while according to formula (33) it should be 2-05. 

Next, to find the effect of the cross-section, carbon dioxide was 
passed through the first and through a third tube whose length 
was about twice and diameter about 1-4 times that of the first 
tube. The ratio of the experimental values of C’ was 1-20, 
while the formula gave 1-15. Whence it appears that the 
current is actually inversely proportional to Z and directly pro- 
portional to a’. 

In order to investigate the effect of the density of the gas, a 
set of 24 parallel tubes were used, each 2 cm. long and about 
0-006 cm. in diameter. 

According to the formula C’ should be inversely proportional 
to V/p, and, therefore, C’Y/p, should be independent of the kind 
of the gas. That such is actually the case, will be seen from the 
following results : 


Gas C” (obs.) V Po CV 
Hydrogen . . . 0168 l 0-168 
Oxygen . .  . 00409 4° 0-164 


Carbon dioxide . . 0-0348 469 0-163 
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To find the influence of the temperature of the gas, Knudsen 
made use of the first tube and filled the vessels with hydrogen 
at a mean pressure of 0-3 mm. In the results of the experiments 
the gas volume flowing through the tube is reduced to the 
temperature 0°, so that 


EN E 
2 ly = +ar) 
Consequently, 
Ya/1 tar 
should, according to the theory, be independent of the tempera- 
ture. Now, Knudsen found : 


m y v 1 + ar 

Temperature C C N cone 
22° 0:0713 0:0713 
100° 0:0641 0-0721 
196° 0-0588 0-0741 


This table shows that also the dependence of C’ on the tempera- 
ture is correctly represented by the formula. 

Finally Knudsen proved the correctness of the constant factor 
in the formula by actually determining in different cases the 
absolute amount of the gas streaming through and by comparing 
the results with those calculated by the formula in question. 
His results were 


Tubes Gas C’ (obs.) C’ (cale.) 
1 Hydrogen . . 0073 0-080 
1 Oxygen . . 0-0187 0-0202 
1 Carbon dioxide . 0-0166 0:0172 
2 Hydrogen . . 0-0375 0-0392 
3 Carbon dioxide . 0:0199 0:0198 
4 Hydrogen .  . 0-168 0-161 
4 Oxygen . .  . 00409 0-0404 
4 Carbon dioxide . 0-0348 0:0344 


With reference to this very rarefied state of the: gas Knudsen 
speaks of a ‘“‘ pure molecular streaming.” He extends also his 
researches to the case of larger densities for which the collisions 
of the molecules can no longer be disregarded and which offers 
considerable theoretical difficulties. In the latter case Knudsen 
speaks of ‘‘ mixed molecular streaming.” 
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14. FLow OF A RAREFIED GAS THROUGH A Narrow ORIFICE; 
KNUDSEN’S EXPERIMENTS 


A second investigation due to Knudsen* concerns the flow 
of a gas through a narrow orifice. 

We assume again that the gas is very strongly rarefied, so 
that the mutual collisions of the molecules can be disregarded, 
and that the dimensions of the orifice are small compared 
with the free path of the molecules. We suppose, further, that 
the temperature of the gas is throughout the same and that the 
diaphragm with the said orifice divides the space occupied by the 
gas in two portions, so that the pressure above the diaphragm is 
throughout p,, and below the diaphragm pe. 

By what was said before the number of particles passing, 
per unit time, through an element do of the orifice from the 
upper to the lower chamber in a direction contained within a 
cone of angular aperture dw can be written 


A, Cos 0 dwdo, 


where @ is the angle between the axis of the cone and the normal 
of the surface element, and A, a factor depending on the tem- 
perature of the walls of the upper chamber and on the pressure 
of the gas therein contained, but independent of the cone chosen. 
Thus the number of particles flying downwards through do per 
unit time will be 


[2 
Í 27A, cos 0 sin Oddo =rA,do, 
0 


and therefore the total number through the whole orifice ¥, 
mAd. 

Similarly the number of particles flying, per unit time, upwards is 
nA,2. 


Whence, the resultant flux downwards, measured by the number 


of molecules, - 
(A, - A,)2. 


Since p y 3 
A= -+ - h=- 
MT T 2u 





* Ann. der Phys., xxviii., 1909, p 999. 
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the last expression becomes 


1 a 
mu 


on Pı = Pe) 2. 


As before, Kundsen expresses ‚the flux through the orifice per unit 
time in terms of the volume the gas would occupy at a tempera- 
ture 7° and a pressure of one dyne per cm.?. If p, be in this case 


the density, then u=v/3/rı (cf. p. 116), and the expression for 


the flux becomes 
1 
a (py -9,)2. 


This formula was again tested experimentally by Knudsen who 
used the same apparatus as in the preceding investigations, but 
having replaced the original tube by a glass tube containing a 
plate of platinum in which an orifice was pierced by means of 
a fine needle. Iie worked with orifices of an area of 0-0005 mm.?. 
According to the last formula the flux is proportional to p- pz 
and can thus be written 
C(v1 ~ Pa); 


the theoretical value of C being %/W2zp,. The experimental 
determination of C proceeds again on the same lines as in the 
previous experiments. In fact, if v, and v, be the volumes of 
the gas at the two sides of the orifice, the volume of the gas 
streaming through per unit time, referred to the pressure of 1 dyne 
per cm.?, is 


d» d 
— v A =V; n= =C(p, - Pr); 


whence 


C= — "W. A log(pı -P Pa) 
ETA At 


The following are some of Knudsen’s results : 


Gas C (caleulated from C (by theoretical 
observations) formula) 
Hydrogen . . . 0225 0:236 
Oxygen . . . . 00565 0-0576 
Carbon dioxide . . 0-0465 0-049] 


These results give, in Knudsen’s words, an experimental 
confirmation of the correctness of the kinetical theory of gases, 
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and especially of the validity of Maxwell’s law of the distribution 
of velocities. If one assumed against Maxwell’s theory that all 
molecules have the same velocity, one would find from the 
theoretical formula for the coefficient C values by 8-6 per cent. 
greater than those of the above table. The observed values of C 
would then differ from the theoretical ones by amounts which 
could not be thrown upon the experimental errors. 

Also in the case of a gas flowing through a narrow orifice 
the phenomena are different at a higher pressure and their theory 
becomes difficult. We have then to distinguish between the 
case in which the mean free path is comparable with the dimen- 
sions of the orifice and that in which it is small compared with 
the latter, which finally leads to the formation of jets. 


15. FLow or A RAREFIED Gas THROUGH A Narrow TUBE 
WHOSE ENDS HAVE DIFFERENT TEMPERATURES 


We will now return to the flow of a gas through a tube 
in order to investigate the effect of temperature differences at 
different places along the tube. 

By what was said before (p. 115), the current is determined by 


2 na? dA 
3 dx 

where ur 
3 

muN Im? 


and it can therefore be written 
_8 a /3m d (P) 
3mN 2 de \w' 
The current is thus dependent on the pressure gradient and the 


temperature gradient. Consequently, the condition for the 
absence of streaming in the tube is 


r- const., 
u 
or p 


Thus we see that, while in a wide tube, in comparison with 
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whose dimensions the mean free path of the molecules is small, 
the condition of equilibrium is that the pressure should be every- 
where the same and, therefore, the density inversely proportional 


with the absolute temperature T (since Zr =const.), in the 


present case the pressure must be proportional to VT and, there- 


fore, the density inversely proportional to vT. 

In a narrow tube whose extremities have different tempera- 
tures these pressure and density differences will arise by them- 
selves. In fact, if originally the pressure is everywhere the same, 
the gas will stream, according to the above formula, towards 
places of highest temperature, so that the pressure will increase 
there. 

This was verified by Knudsen * experimentally. He con- 
nected two MacLeod manometers by a set of tubes of different 
cross-sections. In his first experiments a single capillary tube 
only was used, 9 cm. long and 0-6 mm. in diameter. This capillary 
was placed between two tubes of 14 mm. in diameter. The 
junction of one of the wide tubes with the capillary was brought 
to a temperature of about 350°, which raised the other end of 
the capillary to a temperature of about 100°. The determination 
of these temperatures, however, was uncertain. The two mano- 
meters could also be connected with each other directly by a 
wide tube. The purpose of the observations with a direct 
connection, when the pressure in the two manometers should be 
found equal, was only to form an idea of the accuracy of the 
measurements. The results obtained for hydrogen were : 


Pi Pz 
For direct connection . . . 0-0218 0-0216 
For connection by capillary tube . 0:0223 0-0211 


There is thus, in fact, a pressure difference between the two 
ends of the capillary tube. The whole system shows much 
similarity to a thermo-element. If the two wide tubes are, 
besides the capillary, connected also by a wide tube, the gas will 
continue to flow, and in such a direction that the stream in the 
capillary is from a lower towards a higher temperature. This 
stream is thus analogous to a thermo-electric current. The effect 


* Ann. der Phys., xxxi., 1910, p. 205. 
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will also here be increased by connecting the manometers by a 
large number of wide and of narrow tubes, alternately arranged 
(Fig. 7). At the places marked in the figure by rectangles the 
tubes were heated by means of platinum wires which were brought 
to glowing by an electric current. The temperature difference 
between the heated and the not-heated contact places of the 
tubes was measured by means of a thermo-element and amounted 
to about 500°. The following were the results for hydrogen : 


Pr P: 5 (observed) P sey 
0:00978 0:00419 2°33 2-49 


The pressure difference for different temperatures in different 
parts of a tube filled with a strongly rarefied gas was demon- 
strated by Knudsen in yet another way.* In a glass tube of 


U] | 2 


Fia. 7. 


7-5 mm. in diameter magnesium powder is placed between 
asbestum stoppers. One end of the magnesium column is heated 
electrically up to a temperature of 248°% which gives to the other 
end a temperature of 22°. The pressures at the two ends of 
the tube are measured by means of MacLeod manometers. The 
equilibrium in this case sets in very slowly. Finally we should 


have p _ J T 
P2 T, 
where p, and p, are the measured pressures. 
Knudsen found for p, and p,, after different times : 


Time Pı Pa 
0 0-3000 0-3000 
18h 25’ 0:3308 0:2716 
23h 25’ 0:3347 0-2693 
465 50’ 0-3414 0:2641 
0 0:3428 0-2628 


The last values of p, and p, which would correspond to a perfect 
* Ann. der Phys., xxxi., 1910, p. 633. 
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equilibrium were obtained by extrapolation, it being assumed 
that the relation 


Pe) _ -a| (pı - P2) -4] 


is valid, q being the pressure difference in the state of equilibrium 

From what precedes it is clear that a porous plate should 
show the following effect. If one face of the plate is warmer 
than the other, air should stream from the colder to the warmer 
face, and. this flow should continue so long as there is a tem- 
perature difference between the two faces. 

This phenomenon was demonstrated by Knudsen * in a simple 
way. In a vessel of porous material air is heated electrically. 
Thus the inner surface of the vessel is brought to a higher 
temperature than its outer surface, and air is being drawn 
through the walls into the vessel. If, therefore, the vessel be 
fitted with a tube ending in a bottle containing water, air will 
be seen to bubble up through the water. Knudsen found that 
through a vessel of the capacity of 100 cm.? as many cubic 
centimetres of air could easily be drawn. If the air cannot 
escape, a pressure difference is produced. With a strong heating 
the pressure within the vessel is a few centimetres (of mercury) 
higher than outside. 

Knudsen points out that such phenomena should happen 
often, and that also ın the case of greater gas densities a 
temperature difference will bring about a pressure difference. 
According to Knudsen this phenomenon is undoubtedly of great 
importance in nature. It plays perhaps a röle in the respiration 
of plants and contributes to the refreshing of air in the ‘porous 
soil, air being given out where the surface of the earth is heated 
by the sun, and sucked in where no such heating takes place. 


16. MuTuUAL REPULSION OF Two PLATES AT DIFFERENT 
TEMPERATURES SEPARATED BY A RAREFIED GAS 


Knudsen f has constructed an absolute manometer based upon 
radiometric action. This was suggested to him by an investiga- 
tion on the change of the velocity of gas molecules due to their 
impact against a fixed wall differing in temperature from the 


* Ann. der Phys., xxxi., 1910. p. 207. 
t Ibid., xxxii., 1910, p. 809. 
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gas. This occurs already in Crookes’ radiometer. In this 
apparatus four vanes of mica, covered on one side with lampblack, 
are mounted on the arms of a horizontal cross free to spin round 
on a pivot. If the mutual collisions of the air molecules are 
disregarded, the following explanation of the motion of the 
radiometer seems to be the simplest. Both faces of each vane 
are hit by equally numerous molecules and in like manner, but 
those molecules which impinge against the blackened faces will 
be thrown back with a greater velocity than those hitting the 
clear faces (for, owing to a greater absorption, the former have 
a higher temperature). As a result the vanes will revolve with 
their unblackened faces turned forward. To many physicists, 
however, it seemed doubtful whether this simple theory could 
be here applied. 

Now, Knudsen considers a case for which the theory can be 
adequately developed. He takes two plates S and S’ (Fig. 8) 
immersed in a strongly rarefied gas and placed at a distance 
from each other which is small in comparison with the mean 
free path of the molecules. The plate S has the temperature of 
the surrounding medium, while S’ is brought to a higher tempera- 
ture, so that the plates will exert upon each other a repelling 
force. This force is a function of the pressure. The latter can 
therefore be determined by measuring the former, which gives 
the principle for constructing a manometer. 

We saw from the preceding discussion of the flow of a 
gas through narrow tubes and small orifices that a temperature 
difference gives rise to a pressure difference. Let p be the 
pressure between the plates, po the pressure outside, and T and 
T’ the temperatures of the two plates. Then Knudsen takes for 
the temperature between the plates the mean ¿(T +T’) and 
derives the relation p_vÄT+T) 


Po vT 
If K be the repelling force per unit area, this relation gives 


T+T _ 
Kep-me| | ar -1 |p 


K 
en Er 
rel 


(34) 
"27T 
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The formula used by Knudsen in his experiments with the 
manometer differed from this somewhat ; it was 


2K 
Pom oj + (84) 

ie] 

NE 


The force X is measured by means of a torsion balance to 
which the plate S is attached, while the fixed plate S’ is heated 
electrically. 

Knudsen subjects the formula (34’) to a further scrutiny, 
but we shall not follow him in'this and shall consider the question 
in a somewhat different way. 

We assume that the dimensions of the plates are large com- 
pared with their distance, that the molecules move from one 
plate to the other in all directions according to 
the cos #-rule (cf. p. 110), and that Maxwell’s 
velocity distribution holds. T 7 

Let us now follow a molecule on its path 
through a unit time, during which interval it 
moves many times back and forth, and as many 
times towards the right as towards the left. We S$ S’ 
will calculate the number of flights in a certain Fic. 8. 
direction, say, towards the right, which we will 
denote by v. The number of flights for which the velocity of 
the molecule is contained between v and v+dv and the angle 
of inclination of this velocity to the normal erected upon the 
plates between 0 and 0 +d0, will be found by multiplying the 
number v by the probability that v and @ should fall within 
these limits, that is to say, by the fraction of all molecules 
rebounding, per unit time, from one of the plates, for which the 
velocity and direction of motion satisfy these conditions. 

The number of molecules, per unit volume, whose velocities 
lie.between v and v+dv is proportional to e”*"v2dv, and the 
number of all such molecules leaving per unit time one of the 
plates is equal to Ce ™ vdv, 


where h corresponds to the temperature of that plate. The 
constant C must satisfy the condition 


of eT sdy = 1, 
0 
whence O =2h?. 
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The requirement has still to be satisfied that the direction of 
motion should make with the plate normal an angle included 
between 0 and 8 +d0. The number of particles leaving the plate 
is, with equal cone aperture, assumed to be proportional to 
cos 9. Again, the (angular) space between two cones of apertures 
0 and @+d6 is proportional to sin @d@. Thus, the probability 
that the molecules will have the required direction of motion, 


C” sin 6 cos OdP, 


where C” is determined by the condition 
7/2 
of sin 0 cos 6d9 =1, 
0 
or C’ =2. 


Since these two probabilities are independent of each other, 
the required probability will be equal to their product, t.e. to 


sd bdn, 


Multiplying this by v we shall find the number of passages 
towards the right made by a molecule per unit time under the 
said conditions. The same expression will be found for the 
corresponding number of passages towards the left, provided h is 
replaced by k’, where h’ belongs to the temperature T” of the 
heated plate. 

An expression for the number » itself will now be found by 
noticing that the time taken by a single passage is l/(v cos 0), 
where / is the mutual distance of the plates, and that the time 
taken by all the passages under consideration is the unit of time. 


Thus, 
7/2 %0 -hv w/2 -hr 

Alv | f sin 6d0 f hee ved +. f sin 0d0[ k'e” Pdo} =1. 
0 0 0 0 


Since "het ed 1 m 
| UC v= ie 


this gives 


Ah? sin 6 cos be 


1 
"Wnh +V) 
In order to find the pressure between the plates we have only 


to take the difference of the components along the plate normal 
of the momentum of all the molecules flying away, per unit of 
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time, from the plate S and of all those flying towards the plate S. 
These two components have opposite signs. Consequently, the 
pressure will be given by the sum of their absolute values. 

For the number of times a molecule, within the given limits 
of v and 6, rebounds from the plate S, we have found, per unit 
of time, 

4yh® sin 0 cos Oe "vll. . . . (35) 


With this direction of motion and this velocity the component of 
momentum along the plate normal is mv cos 6. Thus, multiplying 
the expression (35) by mv cos @ and integrating over all values 
of v and @, we shall find the momentum carried by a particle 
which rebounds from S, per unit time. Multiplying this further 
by the total number of particles, we shall have the pressure 
exerted upon the plate by the molecules which move towards 
the right. For the pressure due to the molecules moving towards 
the left a similar expression holds, with A replaced by k’. Divid- 
ing by the area of the plate, we shall have the pressure per unit 
area. Thus, if n be the number of molecules per unit volume 
and | the distance of the plates, the pressure will be 


r[2 = „Ahr! 
p =4vmnl| Í sin 8 cos? @d6 | he" Adv 
0 0 


mie a s -Ke 3 
-} f sin @ cos? 6d0 | h'2e vdy |. 
0 0 . 


Since O hy ; 3 Ir 
e vidv=.. / a.) 
0 8N 


n? 


we have ultimately f 
sal l mn 
= — +: = 
p= ae hiv w) 2Vhh’ 
In order to verify this result we can put h=h’. Since A= fee 
we have 
ana mnu? 
Paap Tg E 
which is the correct value of the pressure when both plates have 
the same temperature. 

Thus the pressure in the space between the plates is deter- 
mined. The pressure p, outside the plates is smaller, and the 
difference p — py gives the repelling force K. 

K 
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The temperature is throughout the greater part of the space 
outside the plates equal 7, and the density is only in a thin 
layer adjacent to the plate S’ smaller than in the remaining 
space. If this thin layer be disregarded, the number of molecules 
per unit volume can be put equal no and the pressure 


MN 
Po =gh ° 
Thus we find 
.„_.m mm _/n h 
K= hie Oh =p( -1). 


nN w (36) 


It remains to establish a relation between n and n,. For this 
purpose suppose first both the plates had the same temperature 7’. 
Then the density of the gas would be everywhere the same. If 
now the plate S’ is heated up to the temperature 7", the first 
consequence will be that at the edge of the plates more particles 
will pass from the space between the plates into the outer space 
than vice versa, until an equilibrium is reached and as many 
particles pass in one as in the other direction. A relation 
between n and no will thus be found by equating to cach other 
the numbers of particles moving in one and in the opposite 
direction. Since, however, the state of motion near the edges is 
complicated, we will simplify the reasoning by an artifice, viz. 
by supposing that there is an orifice in the plate S which is so 
small as not to change the state perceptibly. Then, in the case 
of equilibrium, also the number of molecules passing through the 
orifice into the space between the plates will be equal to the 
number of those passing through the orifice outwards. Let w 
be the area of the orifice. In the first place, the number of 
particles passing through the orifice outwards, per unit time, is 
equal to the number of those which per unit time strike the area 
w of the plate S. By what precedes, this number is 


n 
vnwl Jal Vht VY 

Next, to find the number of particles passing per unit time 

inside, notice that the orifice is made in the non-heated plate, so 

that in the space outside the plates and also at the surface of the 

plate S the temperature can be taken to be T and the number 

of particles per unit volume nọ The number of particles per 
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unit volume having a velocity between v and v+dv and a 
direction of motion between @ and 0 +dd is On, sin de” “ordbde, 
where C is determined by the condition 


© —hv* 
Un, | sin 0d0 [ e wdv=n 


which gives 73 
C = 2 < a . 


The number of particles per unit volume satisfying the said 
condition for v and @ 1s, therefore, 


3 á 
a “ N, sin de" ddr. . . . (37) 


The number of particles striking, per unit time, the area w 
will be the product of (37) into vw cos @ integrated over all v and 
over 6 from 0 to ir, t.e. 


o [k da 9 ahe? Whg 
2 mow} sin 8 cos 6d | e dos -”-. 
T 0 Jo Yyarlı 
Equating this number to that found above for the number of 
particles passing through the orifice towards the left, we find 


PEN won z who 
Vru(vVh+vyh) 2Vah 
or n Vh +V 
My 2h 


Substituting this into (36), we have 


Kap Ge - 1) = ho ol 


or, since h/h' =T'/T, 


K =ip(,/ polyps + (88) 
which is the formula used by Knudsen (cf. (34’), p. 127). 
If 7’ —T is small, this expression, as well as (34), reduces to 


T-T 


K =po -IT (39) 
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Notice that in deducing (38) it was assumed that every 
particle after striking a fixed wall acquires a velocity which 
corresponds to the temperature of the wall. From Knudsen’s 
later investigations it appears that this is not the case, and that 
the mean variation of the kinetic energy is only a fraction of 
what it would be according to our assumption. 


17. KNUDSEN’s MANOMETERS 


_ 2K 
Po < T’ 7 
T 
is based the use of Knudsen’s manometer. Here po is the total 
pressure, including that of the mercury vapour. The formula 


holds for py<zi0 mm. of mercury (4 to 5 dynes per cm.?). For 
greater pressures we can take 


Upon the formula 


9 2K 
Dr E 
NT 


where c>1. The dependence of pọ upon the temperature 
remains as in the original formula, since c is independent of 
temperature. 

The first apparatus used by Knudsen consisted of a platinum 
strip S’ which was fixed and could be heated electrically, 
the temperature being determined by measuring the resistance. 
At a short distance from this was placed a platinum plate 
suspended on the arm of a torsion balance by means of which 
the force K could be measured. The whole apparatus was 
placed under a jar in communication with a Gaede mercury pump. 

The following are the results of some experiments in which 
the plates were placed in air : 


T . . . . . 817 1189 1985 2745 376-2 
To è. . . . . BA BT 247 280 372 
P in dynes per cm.? (cal- 

culated from K). . 2:31 2:45 227 221 2-31 


Since the pressure between the plates was kept constant as far 
as possible, these numbers give a good verification of the formula. 
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The mean of the p,-values is 2:28 dynes per cm.?, while the 
measurement with a MacLeod manometer gives 0-20 dyne per 
cm.?. To the latter must be added the pressure of the saturated 
mercury vapour, for which Knudsen’s previous measurements, * at 
a temperature of about 23°, gave 2-04 dynes per cm.?. Thus the 
total pressure becomes 2-24 dynes per cm.?, which agrees well 
with the pressure found with the new apparatus. Knudsen has 
repeated these measurements with hydrogen. 

The condition that the plate distance should be small 
compared with the mean free path A of the molecules was 
satisfied. In fact, since at a pressure of 1 dyne per cm.? the 
mean free path à amounts to about 10 cm., in the experiments 
under discussion A was about 4 cm., while the distance of the 
plates amounted only to 0-055 cm. 

Other forms of manometers are described in the quoted paper 
by Knudsen. One of these apparatus consisted of a copper 
cylinder with a polished end-surface which served as the heated 
plate. The cylinder is heated by an electric current, and the 
temperature is measured by means of a mercury thermometer 
placed within a cavity of the cylinder. Opposite the polished 
end-surface of the cylinder there is a copper plate attached to 
the bulb of a thermometer which is suspended on the 
arm of a torsion balance. The copper cylinder is 
surrounded by sheets of copper which are not heated, 
and serve as guard rings. By means of a special 
system of pipettes Knudsen was able to introduce into 
the testing space exceedingly small quantities of gas. 

For measurements at high temperatures he used an 
apparatus of platinum. 

Even at very high temperatures the formula (38) 
turned out to agree fairly well with the experimental 
results, while this did not seem to be the case with the 
approximate formula (39). 

Yet another very simple apparatus consisted of a 1 
wide glass tube A (Fig. 9), into which a second glass Fia. 9. 
tube b was sealed. The tube b had on one side, at c, 
an aperture. Opposite this, in the middle of the tube, was 
suspended on two cocoon threads a plate of mica, P. When the 
apparatus is placed in warm water, the plate of mica experiences 

* Ann. der Phys., xxix., 1909, p. 179. 
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a repulsion, due to the difference of its own temperature and of 
that of the outer wall. From the positions assumed by this plate, 
when the apparatus was placed in cold and then in warm water, 
the repulsion, and thence also the pressure of the gas, could be 
measured. After a certain time, when also the innermost tube 
reached the temperature of the bath, the repulsion ceased.* 


18. KNuDSEN’s ACCOMMODATION COEFFICIENT 


In a further paper, entitled ‘“‘ The Molecular Heat Conduction 
of Gases and the Accommodation Coefficient,” Knudsen points 
out that the velocities with which the molecules, impinging upon 
a fixed wall, rebound from it do not correspond to the tempera- 
ture of the wall. He introduces, therefore, an accommodation 
coefficient a, by which the temperature variation which the 
molecules would undergo at the impact, if they assumed the wall 
temperature, has to be multiplied in order to give the true 
temperature variation. The accommodation coefficient depends 
on the nature of the gas and of the wall. The more rough the 
wall, the greater the value of a, tending to 1. 

Let us once more consider the two plates S and S’ (Fig. 8) 
at the temperatures T and 7’. Then the molecules which fly 
away from S towards the right will not have, as previously 
assumed, the temperature T but some other temperature 6, and 
similarly those moving away from S’ towards the left will have a 
temperature 0’ instead of 7’, where @ and @’ lie between T and T”. 

By the definition of the accommodation coefficient we have 


the relations 6’ -0 =a(0' - T) 
0 -0 =a(0 -T'), 
l-a,mı 
whence 0=T toy (Z’ -T), 


r_ rm l-a _ au 
0 =T +y aT T’), 


hea. (I 
0 -0=,2-(T-N). 











tn a 





* For the sake of completeness it may be mentioned that Langmuir (Phys. 
Rev. (2), 1, 1913, p. 337) constructed a manometer based on the properties of 
rarefied gases, by means of which pressures as small as 10-7 mm. of mercury 
could be measured. 

t Ann. der Phys., xxxiv., 1911, p. 593. 
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The accommodation coefficient is here assumed equal for the 
two plates. In calculating the pressure the velocities of the 
molecules can be taken as if a=]1, provided that T and T’ are 
replaced by 6 and 0". 


Whence it follows that in the expressions 


OOM, nd p- 
P ~ Qa/hh! Po= Iho’ 


ħa corresponds to the temperature T, while A and A’ correspond 
to the temperatures 0 and 0’ respectively. 

In order to determine the force K we again assume the 
presence of a small orifice w in the plate S and introduce the 
condition that, per unit time, as many molecules should pass 
through w into the space between the plates as cross it outwards. 
This gives the equation 


—_. On _ ano 
VA Vht+Vh’) 2vVrha 
Thus, 
no Vh+ s/h" 
No j 24 ho 
and 


Vhofı, 1 I ra 
K=p-po= N, + > Ir. (Cf. p. 131.) 


The circumstance that the accommodation coefficient is not equal 
to unity affects only the term - Be : . Since h is inversely 
Vh Vi 


proportional to the absolute temperature, we can put 


C,, C 
h =p: k = F 
Now, T +T'=0+0'. Consequently, in the expression 
l 1 1 
vh + Vh = lv a v0 ) 


0 and @’ can be replaced by T and T’, provided the temperature 
difference is small. Only for greater temperature differences must 
the accommodation coefficient be taken into account. 
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19. Heat ConpucTion IN A RAREFIED GAS CONTAINED 
BETWEEN Two PLATES OF UNEQUAL TEMPERATURE 


We now pass to Knudsen’s investigation on heat conduction * 
and, for this purpose, consider again the plates S and S’ (Fig. 8). 
We have to calculate how much energy (for polyatomic gases 
including also the internal energy of the molecules) is transferred, 
per second, from one plate to the other. We know from previous 
considerations that a molecule flies away from the plate S a 
number » of times per second, where 


] 
=,- 7, 7°: (CF. p. 128. 
ó IVa(vh+vh') Se ) 

The probability that the velocity of such a molecule is con- 
tained between v and v+dv and that the direction of its motion 
makes with the plate normal an angle falling between @ and 
0 +d0 was already deduced and amounts to 


4h? sin 0 cos be” v8d6dv (cf. p. 129), 


where h is determined by the temperature of the plate from which 
the molecule rebounds. 

We now introduce the new condition that the internal energy of 
the molecule should lie between e and e +de. Let the probability 
of this be represented by f (h, e) de, where h indicates that the 
function f depends on the temperature. Since, by Boltzmann’s 
distribution law, the function f is independent of v, the probability 
for all the three conditions to be satisfied is equal to the product 
of the probabilities of each of them, so that the number of times 
a molecule rebounds from S under the said conditions is 


4y sin 6 cos Oh?e~""' vf (h, e)dOdvde. 
Here the function f (h, e) satisfies the condition 


[ r%9=1. 


If e, is the mean internal energy at a temperature corresponding 
to A, oo i 
J ef (h, e)de =€,. 


* Ann. der Phys., xxxiv., 1911, p. 593. 
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Since a molecule leaving the plate carries with it the energy 
Jmv? + e, we find for the total energy transfer, per unit of time, 


— hv? 


Av f “sin 8 cos 6d0 | x nee" v8(4mv® + e) f(h, e)dedv. 


The integration over 6 gives 


vi | | ey hmv? +e)f(h, e)dedv, 
/o -o 


and that over v, 


@ 
mM 
2vhe | (šis +595) f (h, e)de, 
since 


| = ide = a | i ody = Er 


Finally, an integration over e gives 


m &5\_ (Mm — 
avi (Shot) =” (5 te): 
If n is the number of molecules per unit volume and / the distance 
of the plates, the flux of energy per unit area will be 


val +6) = “yada Vuh a): 


The energy gained by the plate, per unit area and unit time, 
from the molecules moving in the opposite direction is given by 
the same expression, only with h replaced by h’ in the last factor. 

Ultimately, therefore, we find for the heat transfer from the 
plate S’ towards the > S 


m 


Va Vh+ Re pter -%). 


Since u? =3/2h, the term m/h can be replaced by $ e,n where eor 
is the mean kinetic energy of translatory motion at a temperature 
corresponding to h. Similarly m/h’ will be replaced by $ Eok 
The heat transfer from S’ to S can thus be written 


TAr En = En enter - 4) 
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The factor 4 shows that in the process of heat transfer the 
energy of translatory motion plays a greater réle than the internal 
energy, which can easily be explained, since with a quicker 
translation the molecules move also more often to and fro between 
the plates. If the temperature difference 7’-T' is small, the 
last expression can be written 


n 4 de, , de\ ,mı 
Tr ZEN sant mr -T. . . (0) 
This holds also for greater temperature differences provided eo 
and e are linear functions of T, which is the case of e, 

The quantities de /dT and de/dT are related to the specific 
heats. If c, be the specific heat at constant volume, and c, that 
at constant pressure, both per unit mass, then 


1 e de ) 


= mdr dT? 
o, -e= 2 o 
pe 38m dT” 
whence = 1 @ de, , de) 
V» m\3 dT dT?” 
ee e A 
pee? m\3 dT drr 


Ultimately, therefore, the expression (40) becomes 


m 


n 
2) al Vh + Vh’) 
which is proportional to the number of molecules per unit volume, 
and thus to the density. 
Knudsen writes for this heat transfer per unit area 


W =(T' -T)pen . . . « (42) 


where p is the pressure (there being, for a small temperature 
difference, no need for distinguishing between the pressures at 
T and T’) and e, is a coefficient introduced by Knudsen which, 
by (41), has the value 


(c,+e)(2’-T), . . (Al) 


€ m aee + Cy) 
pl) 
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Since h=p/2p and mn =p, while Yh + v/h’ can be replaced by 


2Vh, this expression reduces to 
ul: J p 
Eg = / 8 per + Cy). 


If M be the molecular weight of the gas and R the gas constant 
per gram molecule, then 


RT 
P= MP 
R 
and therefore, 
2 I eqgte, 
EON 8r VMT c, -cr 


The heat transfer was thus far expressed in mechanical units. 
To express it in calories, we have to take for e, the last value 
divided by Æ, the mechanical equivalent of a heat unit. With 
E=419.10° and R=83-2. 10° we find 


e l| Spt, 


With a somewhat different deduction of the formula Knudsen 
finds for the numerical factor the value 43-46.10~°. 

We have thus found for the number of calories transferred 
per unit time and unit area from the plate S’ to the plate S 


W=(T’ - T)pe;. 





er =43-44.10- 


The heat transfer appears then to be independent of the 
distance of the plates. 

This formula can also be used to represent the heat transfer 
in other cases, provided the temperature difference 7’ —T is 
small. Thus, for instance, we can imagine a hot body A (Fig. 10) 
kept at the temperature T” and surrounded by an enclosure B 
having a lower temperature T. Since all molecules impinging 
against A arrive from B (mutual collisions of molecules are 
excluded), they will have the energy e corresponding to the 
temperature T, while the molecules coming from A have the 
energy e« corresponding to 7”, the accommodation coefficient 
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being for the present left out of account. The heat supplied by 
A can thus be represented by 


N(e’-e), 


if N be the number of molecules received or sent out by A per 
unit area and unit time. The temperature difference being 

B small, the distinction between the N-values 
corresponding to the temperatures T and 
T’ is of no account, and we can take the 
number which would hold for 7’=T. This 
number will be proportional to the pressure 
and independent of the shape and the 
dimensions of Aand B. The factor e’ - e de- 
pends on the temperature difference T’ -T 
alone and is again independent of the shape 
of A and B. Whence it follows that the 
heat lost by A per unit area should be expressible by the 


sormia W = (T' - T)pex. 


This formula, however, will not hold for the heat gained per 
unit area by B, since some molecules issuing from a certain part 
of the enclosing walls hit another part of these walls and not the 
body A. Only the total loss of heat of the whole surface of A 
will be equal to the gain of the whole surface of B. 


Fia. 10. 


20. THE EFFECT oF ACCOMMODATION UPON HEAT CONDUCTION 


It appears from Knudsen’s experimental findings that the 
amount of heat transferred from one to the other body and the 
coefficient eg calculated therefrom are smaller than their theo- 
retical values. The heat conduction depends on the relative 
dimensions of the surfaces and on their nature. For rough 
surfaces the experimental and the theoretical values of eg differ 
from each other less than for smooth ones. 

To account for this Knudsen considers the effect of accom- 
modation. 

Once more imagine the two plates S and S’ (Fig. 8) at the 
temperatures T and T’. Let @ and 6’ be the temperatures of 
the gas molecules leaving S and S’, and a and a’ the accommoda- 
tion coefficients of these plates respectively. 
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According to the definition of the accommodation coefficient 
(cf. p. 134) we have the relations 


9’ 9 =a(0' -T), 
0 -0' =a'(0 -T'), 


from which follows 


Be nn m 
aaa 
In the expression (42) for the heat transfer 7’ -T has now to be 
replaced by 6’ —@, which can be accomplished by writing 
W=(T’-T)per, 
nel 
“x1 —-(1-a)(1—a’)*™ 


where 


Knudsen points out that what is measured in the experiments 
is this new coefficient eg. Among the special cases he considers 
first that in which both surfaces are of the same kind, so that 
a=a', the coefficient being then denoted by e,,, and then the 
case in which the plate S’ is infinitely rough so that a’ =1 and 
the coefficient is denoted by ej». The coefficients eg, c&u, and 
e1% are related to each other by the equations 


a 
€11 ~ 9 —a Cx, Ely TAER, 


€ x being the value of the coefficient when both plates are infinitely 
rough, so that it might also be written eo». 


The magnitudes À ; ee and | appear to form an arithmetical 
series. EK flo en 

Knudsen determined these different magnitudes by measure- 
ments in which the heat transfer took place between two co- 
axial cylinders. It appeared that the heat transfer depended 
here on the distance of the surfaces and was, for a given pressure, 
greater when the cylinders differed considerably in their radii 
than when they surrounded each other closely. This can be 
explained by noticing that with a large distance between the 
surfaces of the cylinders the molecules which leave the heated 
inner cylinder C’ do not return to C” after a single collision with 
the outer cylinder C, as would be the case for a small distance 
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of the surfaces, but that between two collisions of a molecule 
with C’ several collisions with C can take place, so that the 
molecule has more opportunities to give up heat to the inner 
side of the cylinder C. That the nature of the surface also has 
an effect upon the heat conduction will become plain by assuming 
that a colliding molecule gives up the more of its energy to a 
wall, the more rough the latter. In the case of a perfectly rough 
wall the molecules give up at a single collision their excess of 
heat, and the distance of the surfaces will, therefore, have no 
influence upon the rate of the heat conduction. 

We will consider, then, the two co-axial cylinders C’ and C, 
kept at the temperatures T” and T respectively, and in calculating 
the heat conduction we will follow Knudsen’s not quite rigorous 
method. A perfectly rigorous treatment would offer too great 
difficulties. 

Let us assume that a molecule which flies away from ©’ 
returns to C’ after n collisions with C. Further, let the tempera- 
ture of molecules on leaving C” be @’, after a single collision with 
C, ,, after two collisions with C, @,, and after n such collisions @,. 
If a is the accommodation coefficient for both surfaces, we have 
the relation 

0’ -0,=a(0' -T), 
whence 
9, -T=(1-a)(# - T) =b(0 - T), 


where b=1-a. After n collisions with C we have 
6, -T =b"(0' -T), 

and a collision with C’ gives 
@ —T’ =6(0, -T’). 


From these equations we derive 


, 1-5" ? 
0 um ari (T - T), 
so that 
= 1-5, 
. % E=] _ prti Ex 


This expression gives us the amount of heat carried over from 
C’ to C by molecules which between two collisions with C’ collide 
n times with C. 
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Knudsen splits the heat transfer into parts contributed by 
molecules which between two collisions with C” collide with C 
once, or twice, and so on. For this purpose, however, one has 
first to find the probability that a particle which flies away from 
the outer cylinder C will hit the inner cylinder C’. 

We have always assumed after Knudsen that the number of 
particles which are sent out by a surface-element in directions 
contained within a cone of aperture dw, and with an axis making 
an angle a with the normal to the surface-element, is proportional 
to cosadw. If we now consider a surface-element placed at a 
point P of the outer cylinder, the required probability will be 
represented. by 

wal cos adw 
a: cos adw 





the integral 1 being taken over all cones which have their 
vertex at the chosen point P of the 

outer cylinder and which intersect the r si 
inner cylinder, while the integral 2 X 
has to be extended over all cones, 


with the same vertex P, within a ” 
solid angle 27. We lay through P 
two planes (Fig. 11), PRA, passing 


through the common axis of the Fra. 11. 
cylinders, and PEB, touching the 
inner cylinder. Let QR=¢, 4 QRA =0, ZBRA=9. Then 


do =sin dddde. 


A 


Since a is the spherical distance QA, cos a=sin ¢ cos 0, and 
therefore, 


u I sin? + cos dôd u 
= — - sin d= 
| [ sin? & cos 6dédd 


= 


where r is the radius of the cylinder C’ and R that of C. Thus, 
of all molecules leaving C’ a fraction r/R will again return to C” 
after a single collision with C. The remaining part 1-r/R of 
the molecules will have collided with C more than once. Of 
these again the fraction r/R will return to C’ after a second 
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collision with C. Thus, the part of all molecules which return 
A a and the part 


n-1 
returning after n collisions, if l- 4 ; 


to C’ after two collisions with C is 


The total heat transfer can now be represented by 


W =(T’ —T) per, zs 


where 
ge Si L 1l-(l-a) T \n-1 
ee ae 
while 
C +C l 
‚=43-46 . 10-6 2" ——, . « (48 
=i c,—6, YMT 


This formula was tested by Knudsen in a series of experiments. 
He worked with two glass cylinders, of which the inner one 
consisted of a thin-walled tube around which many windings of 
platinum wire were coiled and fused into the glass. The radius 
of this tube was r=0-340 cm. To the platinum coil two stouter 
platinum wires were attached by means of which the tube could 
be suspended within the wider glass cylinder. In the first 
measurements the radius of the outer cylinder was R, =0-465 cm., 
and in later ones R,=1-6] cm. The outer cylinder is placed in 
melting ice, while the inner one is heated by an electric current 
passing through the platinum coil. The resistance of the latter 
and the intensity of the electric current are measured and from 
these the temperature of the inner cylinder and the amount of 
heat generated are calculated. Since the measurements are made 
after the stationary state is established, this amount of generated 
heat is equal to that of transferred heat. The loss of heat was 
due partly to radiation and partly to molecular conduction. By 
means of a system of pipettes equal small amounts of gas could 
be introduced. This did not change the heat loss due to radiation 
but only that due to conduction, so that from a series of readings 
of the pressure, the temperature, and the amount of heat 
generated the amount of transferred heat per degree of tem- 
perature difference, per unit of pressure and unit of area, i.e. 
the coefficient e, pg, could be computed. The gas chosen was 
hydrogen. 
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Knudsen found for e,z, 1'87.10°° cal. and for e, p, 
2.45.10-78 cal. These results show that, in accordance with 
the theory, the heat conduction increases with the distance of 
the cylindrical surfaces. 

The ratio of e, g, and e,, 7, contains only the accommodation 
coefficient as unknown quantity, so that the latter can be com- 


puted from “nF: One finds a=0-26. 





Cy, R: : 

With this value of a the coefficient ex can be calculated. Its 
value thus found is 11-1 . 10~°, while the, formula (43) based on 
the kinetic theory gives ¢,=11-0.10~°,—a good agreement 
testifying to the correctness of the theory. 

Some further measurements undertaken by Knudsen with 
the purpose of determining the change of the accommodation 
coefficient with the temperature need not detain us here. 

With regard to the accommodation coefficient it may still be 
mentioned that from the theoretical standpoint it is not satis- 
factorily defined. In fact, the accommodation coefficient has a 
meaning only if the state of the gas is completely determined by 
the mean kinetic energy of the molecules, if, e.g., it is assumed 
that Maxwell’s law holds for the velocities of the molecules before 
as well as after their collision with a solid wall. This is for the 
gas between the two cylindrical surfaces not quite the case. 

For a rigorous treatment of the question one would have to 
take into account, for a given temperature of the wall and a 
given velocity of the arriving molecules, the probability of a 
determined state of the molecule after it rebounds from the wall. 
This problem, however, is too intricate. 


21. HEAT CONDUCTION IN A GAS OF GREATER DENSITY 
CONTAINED BETWEEN Two PLATES OF UNEQUAL TEMPERATURE 


After this treatment of the heat conduction in strongly 
rarefied gases we will now give a short account of thermal 
conduction in a gas of greater density. Whereas it is known 
that the coefficient of heat conduction for large gas densities is 
independent of the density, the intermediate domain between 
the very small densities, explored by Knudsen, and great 
densities requires still a detailed investigation. We have 
already touched this domain when treating of a streaming gas 
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layer,* and we have then started from large densities. It 
appeared that with decreasing density next to friction also 
sliding began to assert itself. We will now inquire into the 
behaviour of the heat conduction coefficient in this domain. 

We consider two plates kept at constant temperatures T and 
T’. Let the density of the gas be such that the mean free path 
of the molecules cannot be quite disregarded in comparison with 
the distance of the plates, although there are still very many 
collisions. The temperatures of the gas layers in contact with 
the plates will differ somewhat from T and T’ as was already 
pointed out by Kundt and Warburg. 

Fig. 12 gives a graphical representation of the temperatures 
of different layers. OO’ is the distance of the plates, OA and 

‘A’ represent the temperatures T 
and T’ of the plates, and OB and 
O'B’' those of the adjoining gas 
layers as measured by the mean 
kinetic energy of the molecules. 
For large densities the temperature 
of the gas layers may be represented 
by the straight line 44’. 

At some distance from the fixed 
plates the theory for large densities 
holds, so that the temperature 
gradient is uniform and the tem- 

Fic. 12. perature line has there a straight 

: portion. We produce the latter up 

to the intersection points D and D’ with the vertical lines 

through A and A’. We assume that T’-T is infinitesimal, 

in case the heat conduction coefficient should depend on tem- 

perature. This condition enables us to assume that if the 

temperatures of the plates are increased or diminished by equal 
amounts, the line does not change its shape. 

The whole figure can thus be shifted horizontally without 
changing the length of the lines AD and A’D’. Again, the 
figure can be transformed by increasing or diminishing all the 
distances from AQ in the same ratio. In this case also the 
lengths of AD and A'D’ remain unchanged ; they are thus inde- 
pendent of the temperatures of the plates. lf both plates are of 

* Cf. Chapter II. Art. 10. 
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the same kind, we should find the same curve by interchanging 
their temperatures. This curve will be found by taking the 
mirror image of the original one with respect to AQ. Whence it 
follows that the parts of the curve at B and B’ are congruent, 
and that therefore AD =A’D’. 

Moreover, since the curve retains its shape when the distance 
of the plates is changed,* the length AD=A’D’ for a given gas 
will depend only upon the nature of the plates and will be 
independent of their temperature and their distance. 

Let the distance of the plates be 1 and AD=A’D' =A. 

Similarly to what was said on p. 104 about the significance of 
OC and BD in Fig. 2, a physical meaning can be ascribed to A. 
In fact, if the gas extended also outside of the plates and had 
throughout the same temperature gradient as actually exists 
within the gas, the temperatures at the points D and D’ would 
be equal to those of the plates, 7.e. 7’ and T’ respectively. 

The temperature gradient in the gas at not too small distances 


from the plates is then oe and the heat transfer is k oa 


L+2A 1+2A’ 
where k denotes the conduction coefficient. 

Whence it follows that if the plate distance / is increased the 
amount of heat transferred does not vary inversely proportionally 
to I but slower. 

We will now prove that, for given temperatures of the plates, 
A is inversely proportional to the density while % is independent 
of the latter. 

For this purpose we compare two geometrically similar cases 
[systems] in one of which the gas density is n times greater and 
the dimensions n times smaller than in the other. Let P, and 
P, be two corresponding points, and let the state at these points, 
apart from the density, be exactly the same. Let ds be a volume 
element at P, and Fdsdà the number of molecules contained in 
it, whose velocity components and quantities determining the 
inner state fall within a given domain dà. A similar group in an 


equal element at P, will then contain x Fdsd\. molecules. We 
n 


* That AD=A’D’ is independent of the distance of the plates can be seen 
by inserting between these a third plate P of such a temperature as not to 
change the state of the gas apart from the immediate neighbourhood of P, 
and by applying a similar reasoning to that used on p. 103 to prove that the 
sliding coefficient » is independent of the distance of the plates. 
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assume that the temperatures of the plates are the same in both 
cases; then the temperatures of gas layers at corresponding 
points will also be equal. This harmonises with the assumption 
that at corresponding points the number of the molecules in 
different groups in one gas is the same fraction of their number 
in the other gas; therefore the averages of all quantities are 
equal, and this holds then also for the temperature as the mean 
energy. 

We will first of all prove that the possibility of the existence 
of the first state carries with it also that of the second. We 
consider in the first state a group of particles and follow their 
history during a certain time r. The parameters characterising 
these particles are contained within a given domain, while the 
particles themselves lie in a volume element ds. If there were 
no collisions, we would find these particles at the end of the time 
T in a volume element ds’. Owing to the collisions some particles 
will leave the group and others will join it. The condition that 
the state should be stationary implies that there are in ds’ at the 
beginning and at the end of the time r equal numbers of particles 
whose velocities and parameters determining the inner state fall 
within certain limits. We assume such to be the case for the first 
state and shall prove that this holds then also for the second state. 

The quantities concerning the first state will be distinguished 
by the suffix 1 and those relating to the second state by the 
suffix 2. 

We consider in both cases a group of particles with equal 
intervals for the magnitudes characterising them and contained 
within equal volume elements (ds), = (ds), placed at corresponding 
points. We take the time interval r for the second system n 
times as long. Then also (ds’), and (ds’), will lie at corresponding 
points. If collisions are disregarded, the number of particles 
contained in (ds’), at the beginning of that time interval will be 
n times that contained in (ds’),. This then will also be the case 
at the end of the time intervals 7, and 7,, since all these particles 
arrive from (ds), and (ds),. If there are collisions, then things 
are not so simple. It may be noticed, however, that the number of 
particles which leave the group owing to collisions and the number 
of those that join it for the same reason amount in both systems 
to the same fraction of the total number of particles. In order 
to see this, the collisions may be classified according to circum- 
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stances. Under coinciding circumstances the number of particles 
in the first system is times, and, therefore, the number of 
collisions n? times as large as in the second system. But since 
the time interval for the second system was taken n times as 
large, the number of particles leaving or joining the group in the 
first system will not be n? but only n times as large. The ratio 
of densities will thus remain in both systems the same; also 
with collisions will (ds’), contain » times as many molecules 
as (ds’)s. 

For the interior of the gas it is thus proved that if the state 
of the first system is stationary, so is that of the second. And 
if we further assume that the particles rebound from both plates 
in the same way, this theorem will hold also for the limiting gas 
layers at the plates. 

Lastly, to prove that A is inversely proportional to the 
density, we note that the graphical representation of the tem- 
perature distribution for the second case can be simply deduced 
from that of the first. In fact, since we have assumed that the 
respective temperatures of the plates are equal in both cases and 
that this is true also of the temperatures of gas layers at corre- 
sponding points, the horizontal dimensions of the two figures are 
equal, while the vertical ones of the second are n times those 
of the first. Whence it follows that A,=nA,, which proves 
that A is inversely proportional to the density of the gas. It 
remains only to prove that the conduction coefficient k is 
independent of the density. 


22. HEAT CONDUCTION COEFFICIENT INDEPENDENT 
oF DENSITY 


In order to prove this, we consider the number of molecules 
which pass in equal times through unit area of corresponding 
planes V, and V, parallel to the plates, and we compare again 
similar groups of molecules. 

This number is for the first case n times as great as that for 
the second, and this holds then also for the total energy carried 
across the surface-element. Consequently, the heat conduction 
in the first case is n times that of the second, and since the 
temperature gradient is also n times as great, the coefficient k 
will be the same for both cases. 
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Moreover, noticing that A and k do not change when at 
constant density of the gas the plates are moved farther apart 
or brought nearer to each other, we can say, generally, that A is 
inversely proportional to the density and that % is independent 
of density. 


93. LASAREFF’S EXPERIMENTAL INVESTIGATION 


Lasareff * investigated the temperature distribution in the 
immediate neighbourhood of a wall. He took a very thin gas 
layer, about 9 mm. thick, and worked with highly rarefied gases. 
As walls, metal plates were used which were kept at constant 

temperatures by means of 

water jackets, a cold jacket 

being placed under the lower 

and a warm one above the 

upper plate. The tempera- 

tures were measured by means 

of a thermo-electric pile of 

which one junction was placed 

Sen near the cold plate and the 

Fra. 13. other could be moved up and 

down in the space between 

the plates. Fig. 13 gives a graphic representation of Lasareff’s 

experimental results. The abscissae represent the temperature 

difference relatively to the cold plate and the ordinates the 

distance from the heated plate. The curves AB, CD, EF, and 

GH give the observed results at a prossure of 760, 0-087, 0-065, 
and 0-019 mm. mercury respectively. 

Lasareff found for hydrogen : 


B DFH 


p (mm. mercury) y 
4-5 0-022 
2 0:055 


y is a coefficient, proportional to A. Inasmuch as it can 
be assumed that the temperature of the wires of the thermo-pile 
actually coincides with that of the gas, these experimental results 
prove satisfactorily that A is inversely proportional to the 
density. 

* Lasareff, Ann. der Phys., xxxvii., 1912, p. 233. 


CHAPTER IV 
REMARKS ON LESAGE’S THEORY OF GRAVITATION 


24, LESAGE’s THEORY OF GRAVITATION 


KNUDSEN’S investigations on rarcfied gases may be connected 
with the old theory of gravity due to Lesage.* According to this 
theory celestial space is full of small particles or corpuscles 
moving in all directions with great velocities. Material bodies 
are incessantly hit by these corpuscles and throw them back. A 


| 
ot y bbe 


—> Du 
Tg 
/ \ N 
Fia. 14. 


body placed alone in space will be hit in all directions by equal 
numbers of particles and will thus experience no resultant effect. 
Two bodies 4 and B, however, would partly shield each other 
from the impact of the corpuscles and consequently appear to 
exert upon each other an attractive force (Fig. 14). If the 
corpuscles are assumed to be small compared with the atoms, 
this would lead to an attractive force between two atoms. 

Since the number of corpuscles intercepted by one atom is 
proportional to the solid angle under which it is seen from the 


* Journal des savants, 1764. 
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other atom, the force would be inversely proportional to the 
square of distance of the two atoms. 

But it does not follow from the theory that the force is 
proportional to the masses; it would, instead, depend on the 
dimensions of the bodies. To secure the proportionality of the 
force to mass, yet another hypothesis would have to be introduced. 

Meanwhile the theory of Lesage was shown by Maxwell to be 
incorrect. The motion of the corpuscles is much the same as 
that of gas particles, and as against the fact that the body B 
intercepts corpuscles which in its absence would have reached 
the body A there is this other fact, that due to reflection from B 
some corpuscles will reach Æ which otherwise would not do so. 


AN 


we 3 RR 


Fic. 15. Fia. 16. 


It is not possible to keep the space between two bodies free from 
corpuscles, no more than to keep a space free from black body 
radiation, even if the mean free path of the corpuscles is great 
and their mutual collisions may be disregarded. In the case of 
two parallel plane plates, for instance, only those corpuscles 
which hit the plates perpendicularly would not reach the space 
between the plates; but this is only an infinitesimal fraction of 
the total number of corpuscles. 

It would be difficult to prove that a space can never remain 
free from corpuscles. In fact, one can very well imagine a 
particular state of motion for which a certain domain remains 
free from corpuscles, if, e.g., these move originally outside a 
certain sphere B and inside a second sphere B, concentric with B. 
If the surface of B, is perfectly reflecting and if there are no 
collisions between the corpuscles, then the corpuscles will never 
penetrate into the inner sphere (Fig. 16). But such a state of 
motion will never arise. , 
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The theory of Lesage can be saved by assuming that the 
corpuscles are wholly or partially absorbed by matter. But then 
the picture is deprived of its simplicity. 


95. AN ELECTROMAGNETIC ANALOGUE OF LESAGE’S THEORY 


It was once asked whether Lesage’s theory can be given an 
electromagnetic form. One would then have to assume, for 
instance, that space is full of radiation of a wave-length much 
shorter than that of Réntgen-rays and to show that two particles 
would be driven towards each other by the radiation pressure. 
Such would, in fact, be the case if the particles continually 
absorbed the radiation. 

Consider two electrically charged particles P and Q, at a 
mutual distance 7, in a space traversed by electromagnetic 
waves. Let E be the electric, H the magnetic force, and n the 
frequency. Evidently the particle P will be set into vibrations, 
whether it is free or bound to a position of equilibrium ; only 
the type of vibrations will be different in the two cases. Owing 
to the radiation emitted from Q the radiation field in which P 
is placed will be modified. 

Let us introduce a system of co-ordinates with P as origin, 
and let the co-ordinates of Q be -r, 0, 0. 

In a first approximation the displacement components of P 
can be written 

x=aeE,-beE,, 
y=aeE,-beE,, 
z=ageE, -beE,. 


This is a convenient way of expressing the phase difference 
between the incident waves and the oscillation produced by them. 
The "phase difference is represented by the second terms and is 
caused by the resistance experienced by the particle due to 
friction or emission of radiation. The coefficients a and b are 
supposed to be constants which depend on the mass of P, on the 
restituting force, etc. In E, E, E, is included the field which 
the particle Q, being also set vibrating, produces around itself by 
its radiation. 

Noticing that the work done by the electric force against the 
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resistance, 1.¢. eH,¢ +eH,y +eH,2, for a full period, is positive, it 
will be seen that b must be positive. 

As soon as P is set vibrating, the force upon it is no longer 
E as at the origin [position of equilibrium]. The z-component 
of the force can now be written, in a second approximation, 

(U; -2H,) + (a: +95, +2) +eE,, 
where for H,, H, E, and their derivatives are to be taken the 
values belonging to the original, 7.e. the equilibrium position of P. 
The first term is due to the velocity of the particle, and the second 
due to its deviation x, y, z from the position of equilibrium. 
These terms give rise to very weak oscillations with twice the 
frequency of the incident waves. 

Now, to find the z-component of the force acting upon P, the 
last expression has to be averaged over a full period. The 
required, rather lengthy calculations, in which account must also 
be taken of the disturbance of the field by the vibrations of Q, 
may here be omitted. It appears that, if terms with 1/r?, etc., 
be neglected (which is permitted provided r is large compared 


with the wave-length), the term °(yH , ~2H,) alone survives, and 
c 


this also inasmuch only as it depends on the terms -beÉ in 
the preceding equations. 

As a result of these calculations one finds for the required 
force 


e? .. $ b2n?e? 
-b q (2v8 -Ë,H,) A c ° (EH: -E:H,), 


where for E and H is to be taken the field as it is in presence of 
Q, but neglecting the effect of P upon it. 

This result can be associated with an energy flux. In fact, 
the expression c(Z,H, — E,H,) represents the x-component of. the 
energy flux at P. Since we assume that the rays are propagated 
equally in all directions, there will be in all points of a sphere 
with Q as centre and of radius r the same radial energy flux. If 
E be the flux of energy through the whole spherical surface, 
reckoned positive when directed outwards, then 


E 
o(EyH, ~ E,H,) =» 
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The component of the force upon P taken along the X-axis 
will thus be represented by © 
bne E 


È bart 

Since Q is placed on the negative X-axis, this expression will 
represent an attractive force if it is negative, that is to say, if 
E is negative. Such will be the case if more energy streams 
through the sphere inwards than outwards, and therefore, if Q 
absorbs the rays. 

The electromagnetic modification of Lesage’s theory leads 
thus to a similar result as the original corpuscular theory. 


CHAPTER V 
FRICTION AND HEAT CONDUCTION IN THE PROPAGATION OF SOUND 


26. THE EFFECT OF FRICTION 


WE will now consider the effect of inner friction and of heat 
conduction upon the propagation of sound in a gas. 

Let u, v, w be the velocity components of a volume element 
and s the relative condensation, t.e. if p be the actual density 
and p, that in the original state of equilibrium, s =(p — po)/po. 
Further, let oe pa and let us for the present disregard 

oy 02 
the heat conduction. Then the equations of motion for the 
propagation of sound waves will be 


Krane nn. (44) 
08, 1 oK 
a? ( 
on * a „(Ausg = an 
a 1 OK 
. + _ v(Av Tr, 3 by) 7° (45) 
2 > 10K 
a? 3z T = my (Au + S Oe = (0 


Here a would be the velocity of propagation of sound in absence 
of inner friction, and vy =,/p, where wu is the friction coefficient, or 
the viscosity.* 


* Formula (44) is the equation of continuity for a compressible fluid, viz. 
ep „Apu) | „6(pv) eow) -0, 
ttt 0% vy Oz 
0 
in which the terms u an etc., are neglected, the deviations from equilibrium 


1dp 
being assumed infinitesimal, and the term - 5 A is replaced by - A Te oe 
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Differentiating the three equations (45) with respect to x, y, z 
and taking account of (44) we have an equation for s, 


_ Os = 4 As 

a 3” ai T 

This is, for v =0, the well-known equation of propagation of 
sound waves. 


The problem of propagation in an unlimited gas mass can be 
easily solved for the case of plane waves. 


a?As 0. 200200. (46) 


— — . e = m——n [ln a ad 





The equations (45) follow from the equations of motion 


õu ĉu ðu cu OX, 0X, oX 
-, tU V-— Ft Ww.— | =— = n reo . . 
ae ut ey E =) we Cy wz’ etek = (a) 


(ef. p. 77). 
The expressions for the stress components can be found by following the 
hints given in the footnote to page 78, with the only difference that our fluid is 


. cu ov w 
now compressible so that + +,--+ -_ -K does not vanish. 
ce Ty oz 


Thus one finds 
(uU P 
X, = =p + (Ay) Ay) + aK, 


Uy — Ugg (CV CW 

Ye= "9 (at iy) 

and similar formulae with the same coefficients a,, and a, for the remaining 
stress components. In these equations p is the pressure, as it would correspond 
to the density and temperature at the given point if the gas were at rest or if 
it had throughout the samo velocity u, v, w. 

Between a,, and A, holds, morcover, the relation (a; — 439) + 3a., = 0, which 
can bo deduced by going deeper into the manner of arising of the stresses, but 
which can be accepted here without proof. 

Ultimately, putting a, - Gy =24, we find 

cu 2 
X,= -p+2u 3 uK, 


w iw 
rag) 


and substituting these oxpressions for the stress components in (a), 


(t Wet +w) = ~ Pt usu + ; Ta cte. . . (b) 
Assuming Poisson’s law p=cp*, wo have 
1¢@p_ «nes 
| Pott p òx : 
Again, na, while xa can be replaced by , = ‚so that 
1 Op =q? 08 
p Ca ox 


This, substituted in (b), gives the equations (45), if we put a/p=» and 
neglect products of velocities into their derivatives. 
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If we put s=s,e-%, we can find at once for every frequency 
the corresponding value of g and thence the propagation velocity 
and the damping of sound. 

We will consider here the propagation in a gas contained 
in a cylindrical tube, whose axis coincides with the z-axis. The 
propagation is to proceed along the z-axis, so that all terms will 
contain the factor e™-”, 


Since a =ins and = = —n*s, equation (46) becomes 


(a?+4vin)Astn’s=0. . . . = . (47) 
Since the state in a cylindrical tube is symmetrical about its 


axis, we introduce, instead of y and z, the distance r from the 


axis, so that 
20 „0 „10 
dye dat dee tg OF 
2 
Again, ind =q*s, and therefore, instead of (47), 
073 | lds n? 
or? ror (a + ivn 79 


This is the differential equation of the Bessel function Jy, so that 
we have, for the condensation, 


s=ce™-@] (Ar), . . .  . (48) 


2)s=0. 


2 
ae 
(a? + Swn) 
The argument of I, is thus a complex number. 
The value of u is now to be found from the first of equations (45), 


where A= 


where K = -SE = — ins, so that this equation reduces to 
~ —vAu=q(a?+4wn)s.. .. « (49) 


Of this a particular solution can be found by putting ü=£s and 
by suitably choosing the constant ¢. Substituting this in the 


left-hand member of (49), we have 
" -vAu=£(2 -vAs) =£(ins -vAs), 


so that E(ins —vAs)=q(a?+4wn)s. . .  . (50) 


v PROPAGATION OF SOUND 159 


Since, by (47), Aea s 
asin’ 


we see that equation (50) is satisfied by 
2 
E=9( -i + ty). 


To this particular solution for u an arbitrary solution of the 


equation = -vAu=0 can still be added. 


Since wu is a function of x and r only, the last equation can 


be written 
ou lou à an 
ants opt (Po pen 


and this is again satisfied by a Bessel function, 


in 
u=0ge'"t-] (Br), where B? =q? - a 


Thus the solution for u becomes 
u =e" -éc I (Ar) +c (Br)]. . . . (51) 


It remains only to determine the velocity components v and w. 
Owing to the symmetry the velocity in a cross-section of the 
tube will be radial, so that v and w can be represented by v=yh 
and w=zh, where h will, apart from the factor e“-”), depend 
on r alone. If these values be substituted in the second and the 
third of equations (45), all terms in the second equation will have 
the factor y/r and all terms in the third the factor z/r. On being 
divided by y and z respectively, these equations turn out to be 
identical and after a multiplication by r each of them reduces to 


à 2 
(a? + Kim), + (in -vg?)rh - v(3 a + me =0. . (52) 


The third term of (52) can be evaluated by means of the equation 
of continuity (44) which can be written 


a. - ins +qU 
Oy OZ i 
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Substituting here v=yh and w=zh and differentiating with 
respect to r, we have 

Oh Oh Os Ou 

ae tT aga = aE TT By? 


and this introduced into (52) gives 


a 
- (in —vq")rh = (a? + sim) -qv ae 


Since 7 and ie follow from the solutions (48) and (51) for s and 


u, the radial velocity rh is herewith determined. 

To determine the constants c, and c,, we have two boundary 
conditions. For greater densities, when there is no sliding, the 
gas is at rest at the walls, so that for r =R (radius of the cross- 
section of the tube), u=0 and rh=0. This gives the equations 


éc I (AR) + Cel ( BR) =0, 
(a? + $in -qê AI (AR) - qve,BI (BR) =0, 


where I,’ is written for the derivative of I. 
Eliminating from these two equations ¢,/c., we find 


(a? +$inv - qu£) A 22 + qv&B a Bin =0. 


Here q is the only unknown. Since q is contained also in A and 
B, its determination from this equation is very laborious. One 
has to use approximation formulae for the Bessel function and to 
assume that the effect of friction is small and thus also that 
the state differs but little from the propagation in a non-viscous 
gas. 


27. Errect oF HEAT CONDUCTION 


This problem was treated by Kirchhoff, who has taken into 
account also the influence of the heat conduction. 

In deducing the equations of motion (45) use was made of 
Poisson’s law for adiabatic volume changes (cf. footnote on 
p. 157). If, however, heat conduction is taken into account, 
Poisson’s law can no longer be applied, and one has to write 
down the thermical equation which concerns the change of the 
internal energy of a volume element due to compression and 
heat conduction. 
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Kirchhoff * finds for the propagation velocity in a cylindrical 
tube of diameter 2R 
Eu en: 
v=a(1 cP a=) (53) 
where n is the number of oscillations per second and a the pro- 
pagation velocity in an unlimited three-dimensional space, in 
which case the effect of friction is small; y is a constant 
depending on viscosity and heat conduction, viz. 


evi (p-a) VÉ 


where 
_ Viscosity 


density” 





= Corflicient of heat conduction 
density 


while b is the propagation velocity as calculated by Newton, so 


that ” = 1 Cp 
o Ne 


Formula (53) was repeatedly tested experimentally but was 
never found well corroborated. Kayser t found the deviation of 
the propagation velocity from the value a about four times as 
large as that required by the theory. 

In a space of three dimensions everything becomes much 
simpler than in a tube. The effect of viscosity is then the greater 
the smaller the wave-length. 

Neklepajev ł investigated, in connection with Lebedew’s theo- 
retical researches, § the propagation of very short sound waves in 
air. His source of sound was an electric spark produced at the 
focus of a concave mirror S,. The beam of parallel rays (for 
such short waves, as eg. 0-2 cm., one can speak of “sound 
rays ”’) reflected by the mirror S, fell upon a diffraction grating 
consisting of a series of silvered steel rods. Diffracted bundles 
were thus produced, and one such bundle was concentrated by 
a second concave mirror S, upon a sensitive vane which was 


* Pogg. Ann., cxxxiv., 1868, p. 177; Ges. Abh., Leipzig, 1882, p. 540. 
+ Wied. Ann., ii., 1877, p. 218, 

t Ann. der Phys., xxxv., 1911, p. 175. 

§ Ibid., xxxv., 1911, p. 171. 
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displaced by the pressure of the sound rays. The wave-lengths 
were measured by means of the grating, through the diffraction 
angle. Neklepajev worked with wave-lengths of 2-5 down to 
0-85 mm. 

By sending the sound rays through layers of air of different 
thickness the absorption could be measured. This appeared to 
be considerable. The experiments gave a higher value for the 
absorption than was to be expected from Lebedew’s theoretical 
considerations. Lebedew gives for the distance in which the 
intensity is reduced to rø of its original value, for different 
wave-lengths, the following figures : 


Ain mm. Distance in cm. 
0-8 40 
0-4 10 
0-2 2.5 
0-1 0-6 


For polyatomic gases, in addition to v and & yet a third 
coefficient must be introduced. In fact, when the temperature 
rises while a volume element is being compressed, this will increase 
the velocity of the translational motion of the molecules as well 
as the intensity of their inner motion. In the state of equilibrium 
the energy of translational motion bears a determined ratio to 
that of the inner motion, but while the velocity of translation is 
directly affected by the compression of a volume element, the 
effect upon the velocity of the mner motion is not so immediate. 
The internal energy remains thus in its fluctuations, so to speak, 
behind the energy of translation. This gives a coefficient affecting 
the propagation velocity of sound, viz. making it somewhat 
smaller. It remains to be seen, however, whether this coefficient 
can have a perceptible value. 


CHAPTER VI 


KINETIC THEORY OF SYSTEMS OF ELECTRONS 
RICHARDSON’S INVESTIGATIONS 


28. THEORETICAL INTRODUCTION 


RIcHARDSON* has made an important investigation on the 
emission of negative electrons by a hot metal, and found that 
this is due to the heat motion. The electrons, endowed with 
great velocities, will escape in spite of the forces exerted by the 
metal. The escaping electrons were found by Richardson to 
have a kinetic energy agreeing with that of gas molecules. 

We will, first of all, follow here Maxwell’s considerations on 
the velocity distribution in a monatomic gas acted upon by an 
external force, as e.g. the gravity. 

Let £, 7, © be the velocity components and 2, y, z the co- 
ordinates of a gas molecule. We consider the molecules whose 
velocity components and co-ordinates are contained between the 
limits & and + dé, 7 and » +dn, Cand +d, x and x+dz, y and 
y +dy, z and z +dz. Their number will be represented by 


fdédndfdadydz =fds. 


The components of the external force will be denoted by X, Y, 
and Z. 

All this being valid for the instant t, let us now follow the 
history of this particular group of particles. Mutual collisions 
being excluded, these particles will be contained at the instant 
t+dt within a phase element ds’ placed at a point of the six- 
dimensional space whose co-ordinates are 


x Y Z 
E+ dt, n + 7 dt, G+ dt, xv + Edt, y +ndt, z + Čdt. 


* Phil. Mag. (6), xvi., 1908, p. 353. 
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Let the state of the gas be stationary. This implies that the 
number of particles which fall within the phase element ds’ is 
the same at the instant ¢ as at the instant ¢+dt, so that 


x ' 
Merz... artd,. . „)ds'=fas. 
But by Liouville’s theorem ds’=ds, so that 


X, - Y Z 
J (€ +7 dt,n+ dt, C+ dt, a+ Edt, y + dt, 2+ cat) =f(&,n,¢,2,4,2) 
or, since the difference of these functions is zero, 


fX fY fZ es of, x 
Sm opm fim ae oe oe 


This equation can be satisfied by 
f= ae~ hEn +e) 


where a is a function of the co-ordinates and h is independent of 
them. 


In fact, on substituting we find (with v? = £? + 4? + £2) 
-24 © a(€X +nY +L2) e ann a) 


a(g- - av 5) + (2 -av at) = 0, 


an equation which must’ be satisfied for all values of E, n, ¢, 
z, y, 2, 80 that the coefficients of Ẹ, £3, etc., must all vanish. 
This gives 

oh oh oh Oa 2ha,, 0a 2ha,, da 2ha 

da h gy O a Ben ay me Oe ow 
Thus, the external force must have a potential. If x be the 
corresponding potential energy, 


_ x ox 0x 
er on’ er oy’ Z= 


and 


2 
2 = Ae m, 
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Since f=ae~™*+1"+), we find in the well-known way 3/2h for 


\ 
the mean square of the velocity and a ( : % for N, the number 
of particles per unit volume. 
2h 

This can be written N=N,e ™*,if Ny be the number of 
particles at such places at which the potential of the external 
force is nil. This shows that the density of particles is smallest 
where the potential energy is greatest. 

Since in what precedes no assumption was made about the 
peculiarities of the field of force, apart from the existence of a 
potential, the results arrived at may be utilised in our further 
considerations on the motion of electrons. 

We begin with a limiting case in which the external force 
acts only within a thin layer and is directed normally to the 
boundary of the two media. Let x, x 
be the potential energy in the first, 
and y% that in the second medium. 9 O 
Both magnitudes are constant, while 
the potential energy changes discon- 
tinuously across the boundary. lf 
the medium 1 be a metal, and the Fia. 17. 
medium 2 the space above it, then the 
velocity distribution of the particles in both media obeys 
Maxwell’s formula, the value of h being the same for both. 

Thus the densities will be given by 

2 2h 
a,=ae ™™, a,=ae m*, 

Owing to the force in the boundary layer the velocity of the 
particles will be diminished in traversing this layer. In spite 
of this the mean velocity of the particles will be the same in 
both media, the reason being that only those particles leave the 
first medium which have the greatest velocity. 

To test this result by an explicit treatment, let us intro- 
duce a co-ordinate system of which the yz-plane coincides with 
the boundary of the two media and the positive z-axis extends 
into the second medium. 

We consider a particle [electron] which at the instant ¢ is 
contained in the medium 1 and whose co-ordinates, 2, y, 2, E, m, $ 
fall within the phase element ds. Let 2’, y’, z’, &, 7’, &’ be its 


166 KINETICAL PROBLEMS CHAP. 


co-ordinates at an instant ¢+7, where r is a finite time interval, 
long enough for the particle to have crossed the boundary at the 
instant ¢+7. It will be assumed that it does not collide with 
atoms or with other electrons. 

The condition that the sum of the potential and the kinetic 
energy should be constant gives the equations 


2 
=f-_ (X -x) 
m*a Xa | 0.200. (54 
1 =%); | 
Ç =%. 


Again, noticing that the time required to reach the boundary is 
-z/E, we find the relations 


ı_gr g 

z’ =é (7 + z) | 

y’ =y + m; 

!=2+rl. 
With the aid of (54) and (55) it can be easily verified that the 
state is stationary. In fact, at the instant ¢ there are fds particles 
within the phase element ds and these will have passed at the 
instant ¢ +7 into the phase element ds’, while the latter contained 
f'ds’ particles at the instant ¢. The condition for a stationary 
state will thus be 


(55) 


fds =f'ds’. 
Now it can be readily shown that ds=ds’ and f=f’. In fact, 
| ae! 
a 0 0 0 0 o 
0g’ 
ds a oe 2 
= ee. g 
Ox 

oy 
oy = 
Oz’ 
| a e . ° a az 


En DE EEE, 
DE y OC Ou Oy dz FE" 
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Again, fae Ett) fr ge METH 


and by (54), taking account of the values of a, and ap, these two 
expressions will be seen to be equal to each other. 

We now ask how many electrons pass, in the state under 
consideration, from 1 to 2. We determine first the number of 
particles which, per unit time, pass through a surface-element do 
normal to the z-axis in the first medium, and whose velocity 
components are contained between E and E+d&, n and n+dı, 
¢and ¢+d¢. For this number wefind Efd£dndtdo. Similarly, 
E’f’dE’dn’dt’do’ will be the numbe. of particles which pass, per 
unit time, through a surface-element do’ normal to the z-axis in 
the second medium, and whose velocity components fall within 
the limits & and &’ + dé’, n’ and 7’ +dn’, ¢' and ¢'+d¢’. 

Let us now place do and do’ in the immediate neighbourhood 
of the boundary and assume that £’, »’, ¢’ belong to &, 7, &, t.e. 
that £’, 7’, ¢’ are the velocities, after traversing the boundary, of 
a particle which in the first medium had the velocities £, 7, €. 

It follows then from what precedes that 


f'dg’ =Ed&, n’ =n, C=C, and f’ =f. 
If, therefore, do =do’, then also 
` EfdEdyndlde = &'f'dé'dn'dl'do". 

Herewith is also given the number of electrons which escape 
from the metal. If the state of the first medium is stationary, 
the number of electrons leaving it, per unit time, will be constant, 
even if the state of the second medium is not stationary, which 
is e.g. the case when this medium is unlimited so that the electrons 
are not thrown back. 

At the beginning of the heating the metal plate sends out also 
positive ions, which emission ceases, however, after a certain time. 
These ions come probably from a layer of gas condensed at the 
surface of the plate. 


29. VALIDITY OF MAXWELL’S DISTRIBUTION LAW FOR 
THE FREE ELECTRONS IN A METAL 


The ratio e/m for the negative particles emitted by the plate 
was determined by J. J. Thomson,* who was able to ascertain 


* Phil. Mag. (5), xlviii., 1899, p. 547. 
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that its value was the same as that known to belong to the 
electrons. Thus we have here the same corpuscles. Richardson 
measured their velocity, the arrangement of his experiment being 
as follows. 

Into a quadrangular aperture cut out in a platinum plate 
P (Fig. 18) was fitted the protruding part of the bent platinum 
strip S, the latter being insulated from the plate P by mica. 

p The strip S was heated by an 
? electric current. Another metal 
plate P,, connected with an electro- 
p_ meter, formed with the plate P a 
—___ Lgs condenser. The plate P, gathered 
Fra. 18. the electrons escaping from the 

heated strip S. 

Richardson measured the potential of the plate P, during the 
process of charging, taking at the same time the utmost care that 
the centre of the platinum strip, the glowing spot from which the 
electrons were emitted, should 
remain at the potential zero. 

The velocity with which the 
potential of the upper plate 
mounts in its negative value 
will gradually decrease, because 
owing to the electric field thus 
produced not all the electrons 
leaving the strip S can reach / 
the upper plate (Fig. 19). 

Let us first assume that all Fie. 19. 
electrons leave the strip with 
the same velocity u. There is then a maximum value ®,, to 
the potential of the upper plate. This will be attained when 
the electrons which move normally to the plates are just able 
to reach the upper plate. This gives the condition 


3 2 
ymu* =e®,,, 
whence A mu? 
m 2e s 


We can also find how the potential of the upper plate gradually 
increases. We will assume that the electrons which at the 
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potential ® just reach the plate have left the heated strip at an 
inclination 6 to the normal. These electrons describe then a 
parabola whose apex lies at the upper plate (Fig. 19), and we 
have, for ®, 

mu? co d=eb. . . . = . (56) 


Let now n be the total number of electrons which, per unit time, 
leave the heated strip, and therefore 2” sin 4 cos 0d9 the number 
of those among them which move in a direction contained between 
0 and 0 +d8 (cf. p. 128). Thus, the number of electrons for which 


6 is smaller than @ is . 
2n [sin 6 cos 6d8, 
Jo 


or, in virtue of (56), 
n(1 ae ). 


mu? 
For ®=0 this becomes equal to n, which means that then all 
electrons reach the upper plate. 
If C be the capacity of the condenser, the charge acquired by 
the plate P,, per unit time, is 


OB) 


® (1 emo), 
e 


whence, 


For t=w this gives ®,,, while for ¢=0 we have assumed ® =0, 
The intensity of the current is 


i=0 4P =ne mou, 

Thus, from the measured value of ®,, we can determine mu?/2e, 
and from the rate of increase of the potential the value of 
2ne?/mCu?. From these two magnitudes ne can be calculated 
and, since e/m is known, also the velocity u can be determined. 

We have thus far assumed that all the corpuscles leave the 
hot strip with the same velocity. Let us now see how the above 
result is modified if the velocities are distributed according to 
Maxwell’s formula. It will appear that under these circum- 
stances the potential of the upper plate does not attain a 
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maximum. In fact, when the potential of the plate has risen so 
high that the corpuscles of mean velocity do not reach it any 
more, corpuscles endowed with a much higher velocity will reach 
the plate and charge it to a yet higher potential. 

Let now the upper plate be a circular disc, and let us assume 
the field between the plates to be homogeneous, the potential 
difference of the plates being P. We have to calculate how many 
of the corpuscles leaving an element of the heated strip placed 
at a point O will reach an arbitrary element of the upper plate. 
In doing so we shall assume the heated element to be very small 
compared with the receiving plate, so that all electrons which 
can overcome the potential difference will also be picked up by 
that plate. Let a be the distance of the two plates. If E» n, ¢ 
are the components of the initial velocity of an electron, the 
number of electrons which, per unit time, leave the lower 
plate with velocity components contained between &, and 
Eo +dÉo, n and n +dn, € and ¢+d£, will be represented by 


ag oe ea ta? +9") dé dy dt. 


The coefficient a will be determined by integrating over £, from 
0 to o and over 7 and ¢ from —o to +o and equating the 
integral to n, the total number of corpuscles, which, per unit 
time, leave the element of the hot strip. This gives 


9 
a=" hn. 
T 


The integrated equation of motion of a corpuscle which left the 
glowing element (x=0) at the instant ¢=0 with a velocity &,) in 
the direction of the X-axis is 


umf - 


This gives, for =a, the time ¢ required by the corpuscle to 
reach the upper plate. If the corpuscle actually reaches the 
plate, this equation has two real roots, the smaller of which is 
the required value of ż. 
We solve for 1/t, and have therefore to take the greater of 

the two roots. Thus, 

léo, JEE eð 

t 2a N 4a ma? 


This expression will be denoted by g. 


ed 2 


2ma i ° 
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Let now a rectangular element of the upper plate be hit 
whose sides are dy, dz and which is placed. at a point y, z. Since 
y=nt, z=Ct or n=gy and ¢=gz, the element will be reached 
provided that 7 lies between gy and g(y+dy), and ¢ between 
gz and g(z +dz). 

The number of corpuscles which leave the heated strip with 
the X-component of velocity contained between £, and £)+d£> 
and reach the given surface-element of the upper plate will thus be 


af o€” hto? - hg*(y? te) g2dé gdydz, 


and the total number » of corpuscles hitting that eiement will 
be found by integrating this expression over the values of £, for 
which g is real. Thus, 


v =adydz | E,e~ Mer" +2 gage, o. (57) 
V/Zebim 
This, integrated for y and z over the upper plate, gives the 
required number of electrons as a function of P. Multiplying this 
number by the charge e of an electron we shall have the current 
charging the upper plate, that is, Cd‘b/dt. Richardson calls this 
the thermionic current. 

In Richardson’s experiments the second [upper] plate was so 
large as compared with the mutual distance of the two plates that 
it may be considered as infinitely extended. This facilitates the 
integration of (57). Inverting the order of integrations we find 


i ~h? +0 pgr +0 pnya 
al grée hg, dé, | Š "dy f Š hoz" l 
am 1 j|- am -3> 2, . (58) 
Hl |, ame 7. 





where s = /2e@/m and n is the total number of electrons leaving, 
per unit time, the hot strip. 


The differential equation for the potential of the upper plate 


thus becomes cd m nn... (59) 
dt 
Its solution is Zeh 2 
e™ =] „gehn 2. (60) 
mC 


where we have assumed ® =0 for ¢=0. 
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According to theory ® should thus mount with ¢ continually, 
though exceedingly slowly, the potential tending to become 
logarithmically infinite. Due to unavoidable leakages, however, 
® will practically reach a maximum after some finite time. 

For the thermionic current we find ultimately, by (59) and (60), 


aD en 
1=C — a 
di l+ ie 


so that the current decreases continually, tending to zero for 
t=. 

These theoretical results were tested experimentally by 
Richardson. The value of en is determined from the measured 
intensity of the current at the beginning of the experiment and 


that of - -= from the time rate of change of the current, and 


from these two magnitudes the value of eh/m can be calculated. 
Now, h =3/2u?, where u? is the mean squared velocity of the 
electrons, within as well as outside the metal. Whence, 


eh _ 3e 
m Imu? 


It is now assumed that the mean energy of an electron is 
equal to that of a gas molecule, so that mu? = kT, and therefore, 
eh e 
m 2kT 

Again, if N be the number of molecules per gram molecule, 

eh č eN E 

m NT  2RT 
where R is the gas constant per gram molecule, and Æ (positive) 
the charge of a gram ion of a monovalent electrolyte. The value 
of E is thus known from electrolysis. 

Now, having derived eh/m from his experiments and measured 
the temperature T of the hot metal, Richardson deduced by 
means of the last-written formula the value of the gas constant, 
for a quantity of gas, however, which at the temperature 0° C. 
and the pressure of 76 cm. occupies a volume of 1 cm.?. 

Since the value of R per gram molecule (2 grams of hydrogen) 
is 83-2. 10% and since 1 cm.? of hydrogen at 0° and 76 cm. weighs 
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00000898 gram, Richardson’s theoretical value of the gas con- 
stant, which will be denoted by R, should be 3730. 

The results of the experiments, in which the wire [platinum 
strip] had temperatures from 1473° to 1813° and was heated 
during 16 up to 35 hours, were as follows : 


k=4-1. 10° 
4-2. 10° 
3°5 . 108 
3-6 . 108 
2°9 . 108 
31.109 
32. 10? 
3.4 . 10°? 


In a later experiment 4:04 . 103 was found. 

The satisfactory outcome of these experiments proves that 
the free electrons in a metal have the same mean kinetic energy 
as a gas molecule and that for their velocity distribution Maxwell’s 
law is valid. 


30. VELOCITY DISTRIBUTION OF THERMIONS FOR 
DIFFERENT DIRECTIONS 


In the experiments described in Art. 29 it depends only upon 
the velocity component perpendicular to the heated plate whether 
an electron reaches the plate P, or not. The agreement with 
Maxwell’s law is thus actually proved only for the component £. 
In a second series of experiments * Richardson investigated the 
distribution of the velocity components 7 and & The arrange- 
ment of these experiments was such that the charging of the 
upper plate was mainly due to corpuscles which left the hot plate 
in a slanting direction. 

The electrons are here emitted from a long narrow metal 
strip R (Fig. 20) which fills out almost completely a narrow slit 
in the lower plate P,. This strip is placed perpendicularly to the 
plane of the drawing. The plates are very extended. The upper 
plate P, has a narrow slit BC parallel to that at R. The electrons 
flying through this slit enter into the metal cylinder E which is 
insulated from the plate. The plate P, and the cylinder can be 


* Phil. Mag. (6), xvi., 1908, p. 890; (6), xviii., 1909, p. 681. 
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alternately connected with an electrometer and both can be 
shifted together in a horizontal direction. Thus it can be deter- 
mined in each position what 





E part of the total thermionic 
current passes through the 
A p slit BC. 

BC 2 The measurements are 
i ye made while the electric field 

5 fe is yet very weak. 
7 Let us now assume that 
pa Maxwell’s law holds for all 
Sa eg the velocity components and 
Fra. 20. that for the emission of 


the electrons there is the 
same distribution over different directions as for gas molecules 
flying through a surface-element. In order to calculate the 
number of particles which, starting from an element do placed 
at a point O (Fig. 20) of the hot metal strip, pass through the 
slit BC of the upper plate, we construct around O a sphere with 
OA =z, the distance of O from the upper plate, as radius. Two 
planes laid through O and the edges of the slit BC will cut this 
sphere along two great circles (Fig. 21). Let AB =x, AOB =y, 
and ZYOP=6. The electrons flying across do with a velocity 
u in the direction OP fill 
out a cylinder with do 
as base and the nor- 
mal velocity component 
ucos(AP)asheight. Their 
number is thus propor- 
tional to cos (AP). If we 
now count up the electrons 
flying through do with 
different velocities but all 
in the direction OP, the 
result will still be pro- 
portional to cos (AP), and if this direction falls within a narrow 
cone of aperture sin @d@dy, the number of electrons will 
also be proportional to the latter. Since, in the triangle 
APE, cos (AP) =sin 6 cos y, that number can be written 


c cos ı sin? 0 dOdıb, 
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where c is a constant. Thus the total number of electrons sent 
towards the slit will be 


c cos ob dip f sin? 0d =}rrc cos pdy. 


Now, tan y=2/z, and therefore, paa oe 
cos*v z 


The required number can therefore be written 


‚_.day 22 r 
lro —-( -5-: 4) + 
Bng (z +22 


Richardson’s measurements are in good agreement with this 
result. The number of electrons hitting a narrow strip of the 
upper plate attains a maximum for z=0, while for z = oo it tends 
to zero. 

It is thus proved that for the electrons emitted by a hot metal 
the velocity distribution is the same for all directions ; otherwise 
c would not be constant. 


31. TuE WORK REQUIRED TO DRIVE AN ELECTRON OUT 
OF THE METAL 


Richardson * undertook also some further investigations with 
the object of measuring the work required to drive an electron 
out of the metal plate. For this purpose he investigated how 
the emission of the electrons depends on temperature, having 
assumed in this connection the number of corpuscles per unit 
volume of the metal to be independent of temperature. 

Let us consider a surface-element of unit area placed in the 
second medium near at and parallel to the boundary of 1 and 2. 
In order to determine the number n of corpuscles passing through 
this element per unit time, we introduce again the function f 
(cf. pp. 164 and 165), viz. for the metal 

2h = hv* 
fi=a me 
and for the second medium 


2h, ~ hv? 
Ja =a,0” me 


* Phil. Trans. A, cci., 1903, p. 497. 
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We consider first the corpuscles whose velocities are contained 
between v and »+dv and directions of motion between 9 and 
0 +d0, where 4 is the angle between the direction of motion and 
the normal of the boundary surface. 

The number of these corpuscles passing per unit time across 
the said surface-element is 


2h ~hv? 
ag me Qrv? sin 6 cos Odédv. . . (61) 


Integrating this expression over @ and v we find for the total 
number of corpuscles crossing that surface-element, per unit time, 


~2h, 2/2 ‘%9 A —hvi On 2h, 
n=2ma,e m f sin @ cos ado} ve dv= one " (62) 


If S, be the number of corpuscles per unit volume of the metal, 


we have 
= 00 2h -hv ` _2h 3 
8, =40/ “mye do=a mt, [7 
1= 4r : Ae ve V =A ge ea . « (63) 
whence 
1 Why x1) 
=. . Q@ 7m S,. . « . (64 
A 24/ rh 2 m 
Now, h= so that 
2h 3 1 


m ma ET” a 
We introduce further the quantity 4= -y/e, which will be the 
potential energy for a unit of negative charge, and we write 


e eN E / 
= a- Z ao - = enman 9 e e 65 
kT RT RT en 
where E is the (positive) charge of a gramion of a monovalent 


electrolyte and R the gas constant. 
The formula for n then becomes 


9 E 
n 4/2 ue RIY Yg, 
dr 


Here u is proportional to vT. 

By means of this formula, n being observed as a function of 
temperature, the difference Wa- y, which is positive, can be 
determined. 
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Richardson found from his measurements x, — 1, = 4:1 volts. 
Investigations by Wilson * and by Deininger,} and a second 
determination by Richardson, gave a somewhat greater value. 
The mean of these determinations, 5-54 volts, agrees well with 
the value 5-5 volts found by Richardson § in a later investigation. 


32. Tue DENSITY OF ELECTRONS IN THE METAL 


From Richardson’s experimental findings we may also derive 
the number S, which is here considered as independent of 
temperature. For this purpose we transform the formula for n 
by making use of the relation 


amu? = kT, 
whence follows 


3kT 3 e 3e R 
MEN m N Nmn ESN a p 


The expression under the radical is positive, since e stands for 
the negative charge of the electron. Thus we have 

n=} ze = s —e Bape wg 

From the value of n that of S, can now be calculated, taking 
We — fr, =5-5 volts. Richardson, however, does not give the value - 
of n itself, so that this has to be calculated from other magnitudes 
measured by him. The experimental procedure was as described 
on p. 174 (cf. Fig. 20). The slit and the cylinder E were placed 
at 1 above the middle of the heated strip R. The cylinder £ 
picked then up all the electrons hitting an area of the breadth of 
the slit BC and of the length of 1 cm. The upper plate, the slit, 
and the cylinder could now be shifted in a direction perpendicular 
to the slit BC and thus also perpendicular to the heated 
strip R. In each of the successive positions electrons were 
collected by the cylinder during 30 seconds. The electrons which 
hit the cylinder during a displacement over the full breadth of 
the plate are those which would reach a strip 8, 1 cm. broad, 


* Phil. Trans. A, ccii., 1903, p. 243. 
+ Ann. der Phys., xxv., 1908, p. 285. 
t Phil. Trans. A, ccvii., 1906, p. 1. 
$ Phil. Mag. (6), xx., 1910, p. 205. 
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of the upper plate within 30 seconds; the hot strip R is 
perpendicular to the strip £. 

From Richardson’s measurements the charge carried over, 
per unit of time, from R to 8 is found to be 1-63 . 107° coulombs. 
And since the charge of an electron amounts to 1-5. 10-7 1° coulombs, 
this gives for the number of electrons emitted by R per second 
and hitting 8, 10'°. In order to calculate from this number the 
number of electrons emitted by the metal strip per 1 cm.? and 
per second, we note that the number of electrons flying from R 
towards 8 is equal to that sent out by a rectangle of area 
RA cm.?, if Rand £ be the breadth in cms. of the strips denoted 
by these letters. In Richardson’s experiments, as mentioned 
before, 8 was equal to 1 cm., while the breadth of R amounted 
in some experiments to 0-02 cm., and in some others to 0-04 cm. 
Thus, taking R =0-03 cm., we find n = 10'°/0-03 =33 . 10". 

Again, 

th. -, =5:5 . 108 electromagnetic units, 
E =9650 7 5? 
R =83-2 . 108 5 
—e/m=1:-77 . 10° = 
T =1050 +273 = 1323. 
This gives for the number of electrons per unit volume of the 


metal S, =10%7. 
The result is not very satisfactory, the number found being 


much too large. This is due to the high value of > T (de -Yı) 


which we will denote > by W. In fact, we have found 


n= -W 2 ue — Fate wg. u Zuers, 


Now, in order that the electrons might escape freely we should 

have y.= y,, and therefore, 
nm =const. UN). 

Owing to the factor e7” the number of electrons escaping per 
unit area and unit time will be much smaller than what by the 
last formula would correspond to the number (S,) per unit 
volume of the metal and will decrease considerably with decreas- 
ing temperature. The latter is in accordance with facts. 


In Richardson’s experiments W = 48, so that S, = ~2na/ 3, Im 1, = 
[i.e. 107]. u 
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Such a large number of free electrons per unit volume, however, 
is hard to accept. This is even much greater than the number 
of platinum atoms per cm.3, for which one finds 8. 10°2,* More- 
over, if it be granted that the electrons take part in the heat 
motion and each of them has the mean kinetic energy 3AT, 
such a large number of free electrons would give a much too 
high value for the specific heat of the metal. 

For the number of corpuscles in the second medium we find, 

per unit volume, 


œ _2h —hv* _2hya 73 
> = 4m | Aye m**72e dv = Age m — 
/0 


he’ 
or by (62), 
_ Im la TM y 
IIN a fet 2/5 IET” 


On the other hand, S, can be written 


E(ya-v)) 
S,= 2,/ ype n, . . . (66) 


so that the density of electrons in the metal should be erie í 


times that in the second medium, as we already know. 


33. DIFFERENCE IN POTENTIAL ENERGY BETWEEN THE 
ELECTRONS INSIDE AND OUTSIDE THE METAL ATOMS 


The conclusion of Art. 32 can be avoided by assuming that the 
electrons are partly bound to the atoms and partly free. In addition 
to the force in the boundary layer which hinders the electrons 
from escaping freely from the metal, there is then still another 
force which binds some of the electrons within the metal to the 
atoms. We assume that the potential energy in the metal within 
an atom has a definite value, and we denote by y, the potential 
energy per unit of negative charge for an electron outside, and by 
yr,’ that for an electron within the atom, with the understanding 
that vy,’ is smaller than y,. Further, we denote the part of the 
space occupied by the atoms by a, and therefore the remaining 
part by l-a. 

* In fact, 1 gram atom of Pt weighs 195 gr.; 1 cm.3 of Pt weighs 21-5 gr., 


and thus amounts to } gram atom, while the number of atoms in a gram 
atom can be put at 68. 1022, 
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For the number of electrons per unit volume of the metal 
we thus find 


S,=2 mfa- -aje RT" 4 a0” RIV: ‘lex, . (67) 


JET 


where we have assumed that the medium 1 is in equilibrium with 
the medium 2 and that the density varies with the y’s as indi- 


cated by the expression eRT rep "Our task is now to represent 


the observed facts by means of this formula. Assuming $, 
constant, we have to find how n depends upon T, to compare 
this dependence with the experimental results, and thence to 
calculate yy, and y,'. With these exponents we should then be 
able to deduce again S, from n. 

Two extreme cases may be distinguished. The first, corre- 
sponding to yr,’ =, has just been considered. In this case the 
electrons were subjected only to forces within the boundary layer. 

The second limiting case is that the electrons are acted upon 
by forces only within the atoms and that there are no boundary 
forces. In this case y, will tend to wW. and yy, —y,’ will be much 
greater than W — Yr, so that, even if a be small, the first term can 
be omitted in a first approximation. This gives for S, the same 
formula as in the previous reasoning, but with the factor a, and 
with yy, replaced by y,’. From the experiments it follows then 
again that W.-W,’ =5-5 volts. This, however, gives now for S, 
a value which is only the small fraction a of that found before, 
as if the volume were a times smaller and as if the electrons 
within it were very strongly bound. By choosing a small enough 
the proper value for S, can be obtained. 

One might object to the formula (67) on account of its being 
based upon the assumption of a large number of electrons in 
each atom. That the formula is at any rate correct * can be 
shown by a method of reasoning due to Gibbs.f 

In fact, the number of systems in a canonic assemblage, with 
co-ordinates contained between q, and q, + dg, qa and gq, +dq,,.. . 
qn and q, +daq,, 18 


ed 2. Am 0.0. (68) 


See, however, the remark at the end of this article. 
Tt Fiii Principles in Statistical Mechanics, New York, London, 1902. 


ve KINETIC THEORY OF SYSTEMS OF ELECTRONS 181 


where E, stands for the potential energy which corresponds to 
the given conditions and ® is proportional to the temperature. 
If all mutual action between the parts is disregarded, the potential 
energy of the system is equal to the sum of the potential energies 
of its parts and the expression for the number of systems is split 
into factors, so that (68) assumes the form 


eda dy,dz, e`" dz dy,dza. Kl 


In our case z,, %,, 2 will be the co-ordinates of a corpuscle and Z, 
its potential energy. 

The probability that a corpuscle is situated at a given place 
is thus independent of the remaining corpuscles. The ratio of 
the number of systems in which the first corpuscle is within the 
part a of the space to the number of systems in which it is outside 
that part of the space is then 


ae #2. (1 -a)e ®, >20. (69) 


where Z, is the potential energy of the corpuscle within and E, 
that outside the part a of the space. The same expression gives 
for each system the ratio of the number of corpuscles contained 
within the part a of the space to the number of corpuscles lying 
outside, and thus formula (67) is established. 

In this deduction no assumption whatever was made with 
regard to the number of electrons contained within the part a 
of the space, so that it could even be applied to the case in which 
there are fewer electrons than atoms. What was assumed, how- 
ever, and what may seem objectionable, is that statistical 
mechanics can be applied to the electrons within the atoms. 


CHAPTER VII 
VACUUM CONTACT OF PLATES OF DIFFERENT METALS 


34. POTENTIAL DIFFERENCE FOR VACUUM ÜONTACT 


In connection with Richardson’s investigations the following 
considerations present themselves. 

Suppose we had two metal plates, both at the same constant 
high temperature, and placed in an exhausted space at a certain 
distance from each other (vacuum contact). If the temperature 
is high enough, electrons will be emitted by either plate. If the 
two plates are of different metals, the number of these electrons 
will not be the same for both, and this will produce a potential 
difference of a determined equilibrium value, which must be equal 
to the potential difference between the metals when placed in 
direct contact. For, if such were not the case, the two plates 
could be connected by a wire and, although everything is kept 
at the same constant temperature, we should have an electric 
current which would clash with the second law of thermo- 
dynamics. 

Let the plates be placed horizontally and let the lower plate 
A have a higher potential than the upper, B, the potential 
difference being ®. 

By (58), the number of electrons leaving, per second, the 
lower and reaching the upper plate at a potential difference. 
(which diminishes the velocity of the electrons and thus prevents 
some of them from reaching the upper plate) is 


2evh Ee 
nem =ne" RT * 





a a e EE a nt a n unten menemeni 


* Cf. p. 176, (65) and (65). In (58) $ was the potential of the upper plate 
with respect to the lower one. 
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where n is the total number of electrons emitted, per second, by 
the lower plate. In view of (66), can be written 


aım E 
1 J 2KT ATIS., 


TM 


where r, (the previous ır, -Yr) is the difference in potential 
energy per unit of negative charge within and without the metal 
A, and S, the number of electrons per unit volume of this metal. 

This gives for the number of electrons flying, per second, 


from A to B Eu 
We oer j OS, 


Trin 


while the number of those flying, per second, from B to A is 


9L E 
W wnt TS, 
m 
_ E p è e 
the factor e R? being here left out, since the electric field does 
not counteract the motion of these electrons, so that all of them 
reach the plate 4. 
In the case of equilibrium both numbers must be equal, 
hence the equation 
E 
8, 
Sa 
This gives fer ® 


E E S, 
RT = piP —f) -log s 


The ratio S,/S, is thus in the present case also given by an 
exponential law in which appears the potential energy, split into 
an electrostatic part and one which depends on the difference in 
the attraction between the metals and the electrons. 

The same law for the ratio S,/S, holds also in the case of a 
direct contact between the plates A and B, and since S, and S, 
have fixed values, the potential difference in the latter case is 
the same as for plates separated by a vacuum. 

We can introduce here the same modifications as in the 
preceding investigation by dividing the electrons within the metal 
into two kinds, the intra-atomic and the extra-atomic ones. 

Let the potential energy of an electron in vacuo be by —ep 
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greater than the potential energy of an electron in the metal 
outside the atoms, and by —ep’ greater than that of an electron 
within a metal atom. Since e is always the negative value of 
the charge of an electron, both of these expressions are positive. 
Let, further, Z/RT = u. 

Then, by formula (67), 


2kT 1 
n=} [E am (1 -a)e“? +aer? 2 


For the case of the two plates A and B under consideration 
we find, in view of the potential difference ® (the plate A having 
the higher potential), 


n e~"? =n, 
diog 1, 
uO =log m 


S; (1 = a,)e“Ps + a,e"Ps 
pO log + ao 40,0" 


the suffixes 1 and 2 being attached to all magnitudes belonging 
to the plates A and B respectively. 

In the case of a direct contact between the metals the same 
formula will hold for the equilibrium value of the potential 
difference. 


35. RESISTANCE AT A VACUUM JUNCTION 


Let us now consider a thermo-electric element with two 
vacuum junctions which 


a are kept at the tempera- 
£ tures T and 7” (Fig. 22). 


T T” This will produce. a 
u current, the electromotive 
force being the resultant 

Fra. 22, of two potential differ- 


ences at the junctions 

while the element is open, and of possible potential differences 
in the metals due to the temperature difference. 

Now, we may ask whether this current has to overcome also 

a resistance at the vacuum junctions. To answer this question 
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we assume that the potential difference ® between the plates 
differs very little from the equilibrium potential difference Po. 
From the preceding considerations it follows that the number 

of electrons passing from the plate 1 (of higher potential) to the 
plate 2 ism,e~“*. The number moving in the opposite direction 
is na so that the electric current (reckoned positive from 1 to 2) 
will be 

1 =en, e7? — na). 
Again, 

0 =e(nye-** — na), 
so that 

i =en, (e7 "t —e-H%), 


If u (®-©,) is small, we have approximately 
i= -eun e- "*(® - ®,), 


so that there is, in a certain sense, a resistance at the vacuum 
junctions which is equal to 


This is positive, since e is negative. 


36. PELTIER EFFECT FOR A VACUUM ÜONTACT 


To close this subject, let us still consider the question of the 
Peltier effect for the case of such a vacuum contact. If the 
effect were here the same as in the case of metal plates in direct 
contact, an electric current should give at the place of junction 
(here vacuum) a heat generation when sent around in one, and a 
cooling when sent around in the opposite sense. 

Suppose that within the enclosure containing the plates 
(Fig. 23) everything is kept at a constant temperature T by means 
of a heat reservoir. We have then to calculate how much heat 
must be supplied or absorbed by the latter when an electric 
current passes through the system. Let the current flow from 
the plate 1 to the plate 2. Then the negative electrons move 
from 2 to 1. The heat taken up by the reservoir during a given 
time will be equal to the energy of the electrons which during 
this time enter into the system at P less the energy of those 
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which, during the same time, leave the system at Q, both the 
kinetic and the potential energy being taken into account. 

Since the temperature is the same at P and Q, there will be 
no difference in kinetic energy. Let the potential of the plate 1 
be higher than that of the plate 2, and let their potential difference 
® differ but little from its equilibrium value. Consider a quantity 
of electricity —1 streaming from P towards 
Q. Owing to the potential difference of 
the plates, this negative charge has at 
the entrance P into the system a greater 
potential energy than at the exit Q, the 
difference being just $. In addition to 
this we have also to take account of the 
potential energy of the electrons due to 
the attraction by the metal atoms. The 
state being stationary, the number of electrons contained 
within the atoms will remain unaltered. Thus we can imagine 
that the electrons [involved in the current] pass between the 
atoms. If the potential energy of an electron in the vacuum is 
higher by —ep, than that of an electron in the plate 1 outside 
the atoms, and by -ep, higher than that of an electron in the 
plate 2 outside the atoms, then the negative unit of charge loses 
in passing from P to Q the amount of energy pı -pẹ [This is 
to be added to ®.] 

Thus the total amount of energy lost by the electrons and, 
therefore, gained in the form of heat by the reservoir is, for each 
unit of positive charge passing from 1 to 2, or of negative charge 
passing from 2 to 1, 


Fia. 23. 


O +p, - re. 
We find, therefore, for the Peltier effect r, 
T= D + Py — Po 


37. RICHARDSON’S MEASUREMENTS 


The heat generated in such a system has been measured by 


Richardson.* 
In a vacuum was placed a cold plate P,, kept at a temperature 
T,, and opposite it a hot plate P, whose constant temperature 7 


o 
* Phil, Mag. (8), xx., 1910, p. 173. 
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was so high as to make it emit electrons. In order that all these 
may reach the plate Py, the potential of the latter was kept an 
amount ® above that of P. 

Under these conditions Richardson measured the heat 
generated on the plate P}. 

We consider the plate P, together with the vacuum next to 
P and a part of the wire attached to P, (and consisting of the 
same metal as Py) as a system whose gain in energy has to be 
determined. This will be due, in the first place, to the kinetic 
energy of the electrons entering the system. In order to deter- 
mine this energy we require the mean squared velocity of the 
electrons flying into the vacuum. Now, the number of electrons 
emitted per unit time and unit area by the plate P, with velocities 
contained between v and v + dv, is * 


ar ae~**y8dv, 


where a is a constant. Consequently, the mean squared velocity 
with which the electrons enter the vacuum, 


where u? is the mean squared velocity of the electrons in the 
metal. The mean velocity of the electrons is thus greater at 
their arrival in the vacuum than within the metal, which is to 
be explained by the circumstance that only the electrons endowed 
with the greatest velocities leave the metal.{ 

By what precedes, in conjunction with m/2h =kT, the kinetic 
energy entering the system per unit time can be written 2nk7, 
where n is the number of electrons passing per unit time from 
P to P,. As against this, the kinetic energy of the electrons 
escaping through the wire amounts to 3nkTo. 

Now for the potential energy of the electrons. In the first 
place, owing to the potential difference of the plates the energy 
of the electrons entering the system exceeds by -—en® that of 

* Cf. p. 176, (61), after integration over 6. t Cf. p. 137. 

t This is disregarded by Richardson. The final result, however, remains 
unaffected, since the difference in kinetic energy between the electrons moving 


in and out is, aftor all, negligible when compared with the difference in potential 
energy. 
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the electrons leaving the system (e being always negative). In 
the second place, owing to the attraction by the metal atoms, 
the potential energy of an electron will be greater outside than 
within the system. Denoting, therefore, this excess by —ey we 
shall have for the heat generated on the plate Py 


W =n(2kT ~3kT, —e® — ep). 


[In Richardson’s experiments eb +e} was so much greater than 
2kT —3kTo, that we may as well write 


W = —n(e® +e). 


The plate P, was inserted in one of the branches of a Wheat- 
stone bridge, so that the heat generated could be measured 
through the change of the resistance. Richardson made use of 
the circumstance that W is a linear function of ®. Taking his 
measurements at different values of P and determining the ratio 
of the W’s, he was able to calculate y. 

This quantity appeared to be different according to the manner 
in which the piece of platinum was treated before the experiment, 
and depends probably on the presence or absence of occluded 
gases. The value found for wy by Richardson was again 4 to 5 
volts, and thus agreed with the value found by previous measure- 
ments (cf. Art. 31). 

There remains, however, yet one great difficulty. We saw 
(Art. 33) that the potential difference of 5-5 volts which was 
obtained in Richardson’s experiments described in Art. 31 had 
to be explained mainly by the difference in potential energy of 
the electrons in the vacuum and of the intra-atomic ones, viz. 
by assuming that the electrons are strongly bound to the atoms. 
Now, in connection with the last experiments there is no question 
of electrons in the atoms, so that the potential difference of 4 to 
5 volts must here be explained by the difference in potential 
energy of the electrons contained in the metal but between the 
atoms and the electrons in the vacuum. Thus the agreement of 
the results of the two sets of measurements is rather unintelligible, 
and it must be admitted that not all is clear in these investigations. 


CHAPTER VIII 
PROBLEMS IN WHICH THE MOTION OF ELECTRONS PLAYS A PART 


38. NicHOLS’ EXPERIMENT 


THE presence of free electrons in a metal suggested to Nichols * 
the following experiment. He set a metal disc into rapid spinning 
motion about its axis, expecting the electrons to be driven by 
the centrifugal force towards the rim and thus to produce a 
potential difference between the rim and the centre of the disc. 
Let -e be the charge and m the mass of an electron, w the 
angular velocity of the disc, and ® the potential at a point of the 
disc at a distance r from its centre. Then the condition of 
equilibrium is 

mw?r = —e se 

dr’ 

whence, by integration, 


=- 1 MR + const., 
2e 


or, if P, be the potential at the centre and ®, that at the rim 
(R being the radius of the disc), 


©, - 0,-5.7 aR, 


The question is whether this potential difference can be detected. 
If we take 
w =100 x 27 sec.-!, 


< =1:8 . 10” (electromagnetic units), 
R=10 cm,, $ 
then 
®, —®, =1-1 electromagnetic units =1-+1 . 10-8 volt. 


* Phys. Zeitschrift, vii., 1906, p. 640. 
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Nichols spent much time in trying to detect this potential 
difference, but his apparatus seems not to have been sensitive 
enough. Owing to the rotation of the disc, the measurement of 
the potential difference was connected with great difficulties ; 
sliding contacts had to be applied. He arrived at the conclusion 
that a potential difference 1000 times as great, such as should 
manifest itself if the positive ions were free to move, could be 
detected by his apparatus. The absence of the effect was thus a 
proof that the positive ions are bound.* 


39. Kinetic ENERGY OF THE CONDUCTION-CURRENT ELECTRONS 


The following question suggests itself in connection with the 
preceding considerations. Can the magnetic energy 4L:? of an 
electric current, which is often referred to as the [electro-] 
kinetic energy, be interpreted as the total kinetic energy of the 
electrons, Z(4nv?)? The answer is decidedly in the negative. 
For 3Z2? depends on the self-induction and, therefore, on the 
shape of the wire, and is in common cases much greater than 
S (4mr?). 

The kinetic energy of an electron, whether entirely or partly of 
electromagnetic nature, is at any rate localised in the immediate 
neighbourhood of the electron. The total kinetic energy of all 
electrons is found by addition, since the spaces in which this 
energy is located do not overlap. This gives &(}mv?). On the 
other hand, the quantity 4Z:? represents the total magnetic 
energy to which every electron contributes through its weak 
magnetic field which extends to a considerable distance. In the 
case of an electric current, however, the magnetic fields of different 
electrons are equally directed and will, surely, overlap. Since 
the magnetic energy is a quadratic function of the magnetic 
force, the total magnetic energy will not be obtained by simply, 
adding up the energies due to the separate electrons. Thus, for 
instance, if there are N electrons and if each of them produces 
at a given point of space the same magnetic force H of exactly the 
same direction, the resultant magnetic force will be NH, and the 

* [This experiment has since been carried out successfully, in a modified 
form and with more refined means, by R. C Tolman and others. For its 
recent history and the results obtained see Physical Review, vol. viii., 1916, 


p. 97 and p. 753; vol. ix., 1917, p. 164; vol. xxi., 1923, p. 525; vol. xxii., 
1923, p. 207.] 
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magnetic energy will be proportional to N?, while the quantity 
X(4mv*) is proportional to N itself. 

Since the number of electrons involved in an electric current 
is very large, it will be seen that 174? may be much greater than 
X(danv2). 

For the sake of illustration let us consider a simple case in 
which the conductor consists of a cylindrical wire sheathed by 
a co-axial tube. Let the current i flow upwards in the wire and 
downwards in the sheath. The magnetic force at a distance r 
from the axis is 

1 
Irer’ 





where c is a constant, the propagation velocity of light. If a, 
and a, be the radii of the wire and the tube, the magnetic energy 
in the space between them is, per unit length, 


v? As 
Ane? log oo fe = (70) 


On the other hand, if N, and N, be the numbers of the electrons, 
per unit length, in the wire and the tube, respectively, and v, v> 
their velocities, the kinetic energy of the electrons will be 


tm(N 0 +N). aa 0 


Let us now compare the expressions (70) and (71). 
We have t = N ev; = Nev and, using the electromagnetic mass, 


m= < u (if R be the radius of an electron). Thus (71) becomes 
ôr Re? 
1? 1 1 
Taata N,“ v) 
Now, if a,/a, is a moderate number, (71) will be much smaller 
than (70), provided that 
a ( l, z) 
3R\N, N,/ 
is a very small fraction. 

For that purpose RN, and RN, must be very large, that is 
to say, the number of free electrons in the conductor, contained 
between two parallel planes at a distance equal to the radius of 
the electron, must be very great. 
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Take R=1-5 . 10—*3 cm. and consider a copper wire of 1 cm.? 
cross-section. Then the number of centres of metal atoms 
contained between the said planes will be 14 . 10°. 

This number is so great that the number of free electrons can 
satisfy the requirement of making RN very large, and yet be 
small compared with the number of atoms,—the latter condi- 
tion being indispensable in order to avoid difficulties with regard 
to specific heat. 

Only for extremely thin wires (of a diameter of the order 
of a wave-length of light) would the value of 4/1? be comparable 
with that of 2(4mv*). 
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